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The relative specific ionization along an initially mono-energetic 4-Mev proton beam has been 
measured. The stopping power of aluminum relative to air is found to be 1.48 mg/cm? per cm of 
air at 15°C and 76 cm Hg, and is also found to be independent of the energy in the interval from 
1.5 Mev to 4.0 Mev. That this should be so is adequately explained in the interpretation by 
J. A. Wheeler. The stopping power of a variety of metal foils was also measured and the results 
are given. From the above measurements it is possible to calculate J, the mean excitation energy 
of atoms with many electrons (a constant which enters in Bethe’s formula for the rate of loss of 
energy of heavy particles due to ionization) to be 0.85 Z(me*/2h) or 11.5 Z in electron volts. 


A. INTRODUCTION 


N connection with an experiment on proton- 
proton scattering! in which the small Berkeley 
cyclotron (4-Mev protons) was used, it was 
desirable to obtain a mono-energetic proton beam 
of known energy. The raw beam of the cyclotron 
is inhomogeneous in energy by as much as 
twenty percent.” This difficulty was circumvented 
by appropriately placing three one-mm slits 
about 30 cm apart and in the fringing magnetic 
field of the cyclotron.* Although this decreased 
the intensity of the original beam by a factor of 
about one thousand, there was ample current 
left for scattering and range experiments. To 
determine the energy of this so obtained mono- 
energetic beam, it was decided for reasons of 
expediency and accessibility to measure the 


* The experimental part of this paper was done in 
August, 1940 at the Radiation Laboratory of the Uni- 
versity of California. 

1R. R. Wilson and E. C. Creutz, Phys. Rev. 59, 916 


2R.R Ap Phys. 11, 781 (1940). 
3E. C. Creutz and R. Wilson, Phys. Rev. 59, 916 


ionization in air produced by the protons as a 
function of their range. It is felt that the results 
of these measurements, although terminated 
before their completion, are of sufficient interest 
to warrant publication. 


B. RELATIVE IONIZATION BY HIGH 
ENERGY PROTONS 


Figure 1 illustrates the apparatus and method 
of measuring the relative ionization produced by 
the protons. Because of the difficulty of keeping 
the beam current constant, it was necessary to 
use two ionization collectors as shown; one to 
measure the specific ionization at various points 
along the proton path, and another, which was 
kept in a fixed location, to serve as a monitor ion 
collector. The protons emerge from the vacuum 
system through a one-half mil Cellophane foil 
and first pass through the monitor ion collector. 
Then they pass through air until they reach the 
movable shallow ion collector. This collector 
consists of a thin aluminum foil (2.16 mg/cm?) 
located 0.5 mm in front of an insulated circular 
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brass button of 3” diameter. The foil was raised 
to a potential of a few hundred volts with respect 
to the button, and the ionization current to the 
button was measured with a sensitive galvanome- 
ter. The monitor collector consisted simply of two 
parallel plates (with guard plates), one of which 
was raised to a potential of about a thousand 
volts and the other which led to a galvanometer 
that measured the ion current. The voltages on 
the collectors were adjusted high enough so that 
the ionization currents were very nearly inde- 
pendent of voltage. 

The ratio J of the ionization current in the 
shallow ion chamber to that in the monitor 
collector, normalized to unity at its maximum 
value, is plotted in Fig. 2 for various positions of 
the shallow ionization collector. It is seen that 
the points fall along a typical Bragg curve. To 
determine if the results were influenced by 
heating of the air due to the energy loss of the 
protons (current about 10~'° amp. and of cross 
section 1X2 mm), a run was made with one-half 
the proton current used before. Within the 
accuracy of the measurements, the points were 
unchanged from their previous values. The 
temperature of the air as well as the atmospheric 
pressure was read at each measurement so that 
the results could be reduced to standard values. 
Because of multiple scattering in the air, the 
proton beam, as indicated by the fluorescence it 
produced on a ZnS screen, spread out to nearly 
the full aperture of the movable ion collector at 
its greatest distance from the entrance foil, and it 
is possible that the points near the end of the 
Bragg curve are slightly influenced by this effect. 

The striking characteristics of the curve are the 
pronounced decrease of ionization with energy, 
the very sharp maximum, and the rapid linear 
decrease of the curve after the maximum. The 


curve can be fit adequately, except for the last 
two cm, by an equation of the form 


I=A/(Ro—R), (1) 


where A and R, are constants with values 4.2 and 
27 cm, respectively. The sharp maximum is an 
indication of how mono-energetic the proton 
beam was. The intersection of a straight line 
fitted to the linear part near the end of the Bragg 
curve with the range axis gives the extrapolated 
ionization range. It can be seen that this quantity 
is very precisely determined from such an 
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ionization range curve. However, the relation 
between the extrapolated ionization range and 
the mean range, from which can be obtained the 
energy,‘ is not yet accurately known for protons. 
The difference between the two ranges can be 
calculated if one knows how the ionization varies 
along the path of an individual proton, especially 
near the end of the range. Rado’ has attempted to 
do this with the available data. Inasmuch as 
these data are so widely variable from one 
observer to another, he was unable to reach any 
definite value. His work does suggest, however, 
that 0.6 cm should be subtracted from the 
extrapolated ionization range of 23.2 cm, as read 
from Fig. 2, to give 22.6 cm as the mean range. 
When this range is reduced to 15°C and 76-cm 
Hg pressure, and when the stopping power of the 
foils is added, we obtain 4.00 Mev for the energy 
of the protons as given by the Cornell range 
energy relation.5 The determination of the 
stopping power of the foils will be explained in 
the next section. 

It is felt that obtaining an ionization range 


4M. S. Livingston and H. A. Bethe, Rev. Mod. Phys. 
9, 264 (1937). 
( om Thesis, Massachusetts Institute of Technology 
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curve would be a convenient and accurate method 
of determining the energy of a proton beam if a 
relation between the energy and the extrapolated 
ionization range existed. It is hoped to obtain 
such a relation experimentally by using the 
magnetic analyzer developed at Princeton for 
measuring proton energy. . 


C. Tue StoppinG POWER OF ALUMINUM AND OF 
OTHER METALS 


A rotating foil wheel was located as shown in 
Fig. 1 so the protons could pass through any one 
of the several aluminum foils of different thick- 
nesses. The thickness of each foil was determined 
by measuring its area and weight with an accu- 
racy of about one percent. The area of each foil 
was about one square cm. A chemical analysis of 
the aluminum from which the foils was cut 
showed it to be 99.35 percent aluminum, 0.35 
percent iron, and 0.25 percent silicon.® 

With each foil in place, the same type of 
measurements were made as described in the 
previous section. The results of this series of 
measurements are shown in Fig. 3. All the data 
were taken at one time. The curves all seem to 
have the same shape as that obtained with no 
foil, except that they are shifted to the left. This 
indicates that the straggling does not differ 
appreciably in air and aluminum. In Fig. 4 the 
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stopping power of each foil, determined by sub- 
tracting the extrapolated ionization range ob- 
tained with the foil from that without a foil, is 
plotted. The extrapolated ranges were reduced to 
the same temperature and pressure (300°K and 
73.31-cm Hg) before the subtractions were made. 
Figure 4 shows that the stopping power of each 
foil when plotted against its weight per mg falls 


6 The foils and analysis were kindly supplied by the 
Reynolds Metals Company. 


surprisingly close to a straight line that passes 
through the origin. This is interpreted as meaning 
that the stopping power relative to air is inde- 
pendent of the energy in the interval from 1.5 to 
4.0 Mev. The slope of the straight line of Fig. 4 
gives the specific stopping power to be 1.48 
mg/cm? of aluminum per cm air when reduced to 
15°C and 76 cm Hg. 

In the same way, the stopping powers of 
several other metal foils were measured. In these 
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cases only one foil of each metal was available. 
The ionization range curves are shown in Fig. 5, 
and in Table I is listed the specific stopping 
power for each of the foils measured and the 
quantity g which is explained in the interpretation. 


D. INTERPRETATION OF THE RESULTS’ 


Bethe has given a well-known formula’ for the 
rate of loss of energy per unit of path length by 
heavy particles of non-relativistic velocity. For 
protons his relation reduces to 


—dE/dx = (4rNZe*/mv*) In(2mv?/I). (2) 


Here v is the speed of the particle, and e and m 
represent the electronic charge and mass, N gives 
the number of atoms per unit of volume, Z is the 
atomic number of the stopping material and J isa 
certain mean excitation energy characteristic of 
the atom in question. The value of J varies with 
atomic number in a complicated way for light 
atoms. On the other hand, for atoms containing 
many electrons, Bloch has shown’ on the basis of 
? Discussion contributed by J. A. Wheeler. 


8H. A. Bethe, Ann. d. Physik 5, 325 (1930). 
*F. Bloch, Zeits. f. Physik 81, 363 (1933). 
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the statistical atomic model that the mean 
excitation energy should be proportional to the 
atomic number, Z. No one has as yet computed 
from first principles the constant of proportion- 
ality,'"° which therefore has to be found by 
comparison with experiments. This Bloch did.°® 
He found =0.961(me*/2h?) = 13.1 ev. Unfor- 
tunately the velocities at which measurements 
were available in 1933 were limited to 1.92 X 10° 
cm/sec. (7.6 Mev alpha-particle; same speed as 
K electrons of element Z=9). In the meantime 
further theoretical work of Livingston and Bethe‘ 
has emphasized the importance of the correc- 
tions which must be introduced into the simple 
relation (2) when the velocity of the par- 
ticle is low enough to be comparable to the 
speed of the bound electrons in the stopping 
material. One would arrive at a misleading esti- 
mate of Bloch’s constant if he were to look apart 
from such corrections even in discussing the 
present experiments of Wilson, performed with 
protons of velocity ranging up to v=2.77X10° 
cm/sec. Thus if Bethe’s simple formula applied, 
and I were exactly proportional to Z, one would 
expect there to exist a linear relationship between 
relative stopping power and logarithm of atomic 
member : 


( —dE/NZdx) material 
(—dE/NZdx) sir 


In(2mv?/constant) —InZ 


(3) 


In(2mv?/constant) — In7.22. 
The left-hand side of relation (3) may be repre- 


10 An ap roximate theoretical estimate was, however, 
made by H. Jensen, Zeits. f. Physik 106, 620 (1937). 


sented by the letter q: 
(—dE/NZdx) material 

ai 

(NZdx) sir 


equivalent stopping 


atomic weight of material 


twice its atomic number 


(4) 


( mg/cm? of air 


mg/cm? of equivalent stopping: 


This quantity has been computed from the experi- 
mental data in Table I and is plotted in Fig. 6 as 
a function of InZ. It is seen that q is a practically 
linear function of InZ. Without further con- 
sideration one would be led to the false conclusion 
that he could obtain a fairly accurate value 
of Bloch’s constant from the intercept of the 
straight line with the horizontal coordinate axis: 
In(2mv?/constant) = 6.40, whence I/Z=effective 
mean of initial and final energy of proton, 
divided by 276,000. Actually, however, the mean 
energy of the protons employed in the experi- 
ments differed considerably (2-4 Mev) from 
element to element, so that on the basis of 
Bethe’s simple theory one should not expect the 
points in Fig. 6 to come at all close to a straight 
line, as they do. Moreover, for a definite element, 
such as aluminum, the stopping power relative to 
air is seen from Fig. 4 not to vary with energy as 
it should according to (3). 

Both of the difficulties just mentioned are 
resolved when one notes that the speed of a K 
electron of aluminum is 2.8X10° cm/sec., or 
practically equal to that of the proton. This 


TABLE I. Specific stopping power. 


SPECIFIC 
TOTAL STOPPING 
STOPPING POWERIN ENERGY 
MG BY Fort MG INTERVAL 
Foi or INCM PERCM In MEV 
Al 17.30 11.72 1.48 1.5-4.0 0.864 
Cu 25.35 12.55 2.02 2.6-4.0 0.665 
Fe 7.47 4.07 1.83 3.7-4.0 0.718 
Mo 27.85 11.87 2.34 2.7-4.0 0.597 
Ni 11.48 5.97 1.92 3.4-4.0 0.668 
rm 17.80 5.40 3.30 3.4-4.0 0.465 
Ta 47.20 14.72 3.20 2.2-4.0 0.475 
Zn 8.29 4.13 2.01 3.7-4.0 0.668 
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circumstance requires the introduction of the 
corrections of Livingston and Bethe‘ into the 
right-hand side of (3), which then, quite inde- 
pendent of Bloch’s theory of excitation energies, 
becomes : 

In(2mv?/Ta1) —(Ci/Z) ar 


—(Cp/7-22) air 


(5) 


For J,i, Livingston and Bethe find 80.5 ev or 
5.92(me'/2h*) from the observed absolute stopping 
power.’ Their Fig. 28 gives the corrections Cr. 
Wilson’s Fig. 4 shows that 27.8 mg/cm? of Al are 
equivalent to 20 cm air (300°K, 75.31 cm pres- 
sure). This comparison gives for the ratio (5) the 
value (26.97/26)(20 X1.166/27.8) =0.870, con- 
stant over the range of velocities investigated. 
These data are used in Table II to compute Ja). 

The values in the last column come out to be 
strikingly consistent with one another. Since 
corrections for binding of the LZ electrons will be 
least important at the higher energies, it is 
reasonable to take 8.83 as best value of 
In(2mc?/I,,). This result gives I,;=150 ev 
= 11.0(me'/ Table II suggests why the stop- 
ping power of aluminum is constant relative to 
air over a considerable interval of energy: the 
correction for binding of the K electrons is 
decreasing for one element, increasing for the 
other. The complicated nature of the compensa- 
tion shows how difficult it is to give any simple 
interpretation to the linear relationship in Fig. 6. 

TABLE II. Determination of mean excitation energy of 


aluminum from stopping power relative to air observed by 
Wilson for protons of various energies. 


ENERGY In NUMERATOR In 
MEV (Ck/7.22)air OF (5) (Ce/Z) (2me®/T 


2.0 4.00 0.11 3.38 0.04 8.87 


2.5 4.22 0.10 3.58 0.05 8.86 
3.0 4.40 0.10 3.73 0.05 8.83 
3.5 4.55 0.09 3.87 0.06 8.83 
4.0 4.69 0.08 4.00 0.06 8.82 


As basis for fixing the constant of proportion- 
ality between J and Z there is really available at 
present only the data for air and aluminum: 


Air, I/Z=11.15 ev=0.820(me*/2h*), 
Aluminum, J/Z = 11.54 ev =0.848(me'*/2h?). 


Since Bloch’s statistical theory of the proportion- 
ality is the more nearly correct the higher is the 
atomic number, it seems to be reasonable in the 
light of the evidence to write for the mean excita- 


06} 

+ 

0.4) 

° 2 53 4 > 7 
ln z 
Fic. 6. 


tion energy of atoms with many electrons 
I/Z=11.5 ev=0.85(me'/2h*). 


The constant in (6), about 10 percent lower than 
that given by Bloch, is employed in the following 
paper to determine the range-energy relation for 
fast particles in lead. 

In conclusion the author takes pleasure in 
expressing his gratitude to Professor Ernest O. 
Lawrence for the privilege of using the cyclotron. 
I am especially grateful to Professor John A. 
Wheeler for pointing out the importance of the 
measurements and for contributing the interpre- 
tation. Acknowledgment is made to the Research 
Corporation for financial support in this project. 


| 
| am 
al 
te 
> 

| 
e 
ly 
n 
yf 
e 
t | 
t, 
| 

4 
re 
K 
or 
‘is 

| y 

q 
965 
18 
397 
968 
165 fi 
175 
568 


DECEMBER 1, 1941 


PHYSICAL REVIEW 


VOLUME 60 
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Observations on the curvature of the track of a meson before and after traversal of a known 
amount of matter allow a determination of the mass of the particle. A discussion is presented of 
the accuracy and advantages of this method of measuring mass, the principle of which is well 
known but has so far been little used. In this connection the stopping power for fast particles 
has been computed on the basis of the experiments reported by R. R. Wilson in the preceding 
paper. Discussion of the formula used is contained in the appendix. Curves are presented for the 
dependence of range on energy and momentum for swift hydrogen and helium nuclei as well as 
for mesons. The method of momentum loss is illustrated and values are obtained for the mass of 
the meson from measurements published by Anderson and Neddermeyer in 1934. The evidence 
available from these mass determinations and those published by other authors shows the 
necessity for many more mass determinations before one can decide whether there is a distri- 


bution of masses of the meson or not. 


INTRODUCTION 


EASUREMENT of the range of a meson 

and of the curvature of its track by a 
magnetic field, i.e., its momentum, are well 
known to suffice for a determination of its mass. 
For this purpose use is made of the theoretical 
relation between initial momentum and _ total 
path length. In a variant of this method one 
applies the differential form of the momentum- 
distance relation, and measures the change in 
magnetic rigidity of the particle over a small 
portion of its path. This procedure has been 
employed by several observers, and Corson and 
Brode! have discussed it in some detail. 

The longer is the segment of path length con- 
sidered, and the higher is the atomic number of 
the traversed material, the more surely one can 
exclude the possibility that the particle in ques- 
tion is an electron. But then one can no longer 
directly apply the differential expression for the 
rate of change of momentum with distance. 
Still it is in principle possible also in this case to 
determine the mass of the particle, as Nedder- 
meyer and Anderson? remark in a recent review 
of methods of measuring this quantity. From an 
experimental point of view it is only necessary 
to observe the curvature before and after 
traversal of a plate of known thickness. How- 


1D. R. Corson and R. B. Brode, Phys. Rev. 53, 776 


(1938). 
2S. H. Neddermeyér and C. D. Anderson, Rev. Mod. 


Phys. 11, 200 (1939). 


ever the procedure in question has been but 
slightly used. 

This way of determining mass by the method 
of momentum loss appears to deserve further 
emphasis for three reasons: 

(a) The necessary measurements of curvature 
of track in a magnetic field and of thickness of 
material traversed are among the most direct in 
the whole field of cosmic-ray research. Difficult 
determinations of the density of ionization along 
the track are avoided. 

(b) Accuracy in this method, as in almost all 
other methods of determining the mass, requires 
that the velocity of the meson be not close to 
the speed of light. In the acceptable range of 
velocities, however, practically every meson 
whose curvature can be measured before and 
after penetration of the plate will permit an 
evaluation of its mass. Contrast this attainable 
abundance of data with the rarity of those 
collisions of a meson with an electron in which 
the energy transfer is sufficient to allow a 
determination of mass from the laws of con- 
servation of momentum and energy. 

(c) The law for the stopping power of matter 
for fast particles has an exceedingly reliable 
theoretical basis.* 

The question whether all mesons have the 
same mass is so important that it appears worth 
while to give a qualitative discussion of the 


3 See in this connection E. J. Williams, Phys. Rev. 45, 
729 (1934). 
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resolving power of the “method of momentum 
loss.’’ If the energy of the meson is large in 
comparison with its rest energy, the loss of 
energy, AE, in a given thickness of material will 
be practically independent of energy, according 
to the theory of stopping power. In the rela- 
tivistic relation between losses of energy and 
momentum, AE~vzdp, the velocity v will be 
essentially a constant, the speed of light. We will 
therefore find fast particles losing momentum of 
the same order of magnitude whatever be their 
mass. Consequently the resolving power is poor 
for mesons of relativistic velocity (momentum 
>100 Mev/c; Hp>4X105 gauss cm). 

This insensitivity of the momentum loss at 
relativistic velocities to the mass of the particle 
conversely makes it possible to test the law of 
stopping even for mesons whose precise masses 
are unknown. It is only necessary that their 
energies be sufficiently great : 7p> 10° gauss cm. 
Measurements under this condition of the change 
in curvature on traversal of a plate were made by 
J. G. Wilson* and were found to verify a theo- 
retical expression for the absolute rate of loss of 
energy not very different from (1) below. 

For non-relativistic speeds the rate of energy 
loss per unit of distance may be written very 
roughly in the form 


dE/dx ~constant 


From this it follows that the mass, yp, of the 
particle is given by 


u~ p(dp/const. dx)!. 


The percentage error in yw will be given by the 
sum of (a) the percentage error in the mo- 
mentum, ~p, or corresponding magnetic rigidity, 
Ip, of the particle; (b) one-third the percentage 
error in the measurement of the momentum loss 
(if that loss is small compared to p) ; (c) one-third 
the percentage error in the thickness of the 
stopping material; (d) one-third the percentage 
error in the theoretical constant of the stopping 
power formula. 

These considerations on the whole favor the 
“method of momentum loss” for non-relativistic 
velocities. 


4 J. G. Wilson, Proc. Roy. Soc. A172, 517 (1939). 


RANGE-MOMENTUM RELATION 


A more complete treatment of the method in 
question requires a discussion of the law of 
stopping power. Fortunately, radiative losses can 
be neglected for mesons of much greater than 
electronic mass as for protons and alpha-par- 
ticles, even for the highest velocities (v=0.9c) 
and greatest atomic numbers (Pb; Z=82) which 
will be of interest. The so-called “radiation 
unit” or ‘‘unit length” of the radiative theory 


will have a value for mesons greater than that 


for electrons (0.4 cm in the case of Pb) approxi- 
mately in the ratio of the squares of the two 
masses. Thus the mean energy dE, lost in the 
form of radiation by a meson of 150-Mev energy 
and of mass 200 m traversing 1 cm of lead will 
be given by® 


dE,/E~dE,/150 Mev ~(1 cm) In2/(200)?0.4.cm, 


whence dE,~0.01 Mev. This value is trifling in 
comparison with the loss of ~10 Mev by ioniza- 
tion, which alone therefore need be considered. 

According to the theory of stopping power, 
the rate of dissipation of energy through ioniza- 
tion and excitation by a fast particle of velocity 
v and of not too great charge number, z, will be 
given by the relation,® 


dE/dx = (4xNe'z?/mv*) 
X 2mv?/T(1 —v*?/c?) ]— }. (A) 


Here N represents the number of electrons per 
cm* of stopping material, J is a certain mean 
energy of excitation for these electrons and m is 
the electronic mass. For atomic hydrogen,’ 
I =1.103(me'/2h?) ; for helium,’ J = 3.19(met/2h?*). 
For air (Z=7.22), Livingston and Bethe® deter- 
mined the mean excitation energy from the 
observed stopping power : J =5.92(me*/2h*). The 
experiments of R. R. Wilson reported in the 
preceding paper permit a similar evaluation of 
the mean excitation energy for aluminum :'® 
IT=11.02(me'/2h?). Elements of high atomic num- 


’ For a more detailed treatment of the radiation loss, 
see H. J. Bhabha, Proc. Roy. Soc. Al64, 257 (1938). 

® For literature and discussion see Appendix. 

7H. A. Bethe, Ann. d. Physik 5, 325 (1930). 

8 E. J. Williams, Proc. Camb. Phil. Soc. 33, 179 (1937). 

®M. S. Livingston and H. A. Bethe, Rev. Mod. Phys. 
9, 267 (1937). 

10R. R. Wilson, Phys. Rev. 60, 749 (1941) and discussion 
at end of Wilson’s paper by J. A. Wheeler. 
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(MOMENTUM c) “ LIGHT VELOCITY TIMES MOMENTUM CHANGE (Mev) 
so oc 
VELOCITY OF LIGHT 
TIMES CHANGE IN MOMENTUM (Mev) | 
: Zo2 Zo. 
10% TIMES CHANGE 
in He (Gauss cm) 
= © 
'N 


R Lo 12 Le 18 20 22 24 


Fic. 1. Calculated relation between range and momentum for hydrogen and helium nuclei and for mesons. For the 
stopping materials computations were based on values of the mean ionization energy mentioned in the text. For Pb and 
Al the ionization energy was determined from Bloch’s formula with the constant of proportionality based on the data of 
R. R. Wilson in the preceding paper. To interpolate for the stopping power of a compound substance such as PbCl: note 
that the logarithm of the effective atomic number is (82 In82+34 In17)/(82+34). Given that a meson (s=1) with 
a given magnetic rigidity, Hp, of 528,000 gauss cm passes through a 3.70-cm lead plate, corresponding to a range num- 
ber of 9.92, and emerges with a rigidity of 238,000 gauss cm. We compute the change in Hp and plot the point A 
in the upper diagram. The secant CD must be parallel to OA but the positions of C and D individually remain to be found. 
A large scale transparent model is constructed of the template in the lower right corner of the figure. Lines SC and SD 
are drawn such that D’B’ : C’B’=DB : CB=528,000 gauss cm: 238,000 gauss cm. = 2.22. With its rulings parallel to OA 
and with corresponding points C’D’ or CD, etc., kept on the curve for lead, the template is moved until the prolongation 
B’ of C’D’ or B of CD, etc., lies on the line QR extended. The graphical solution so found is indicated by the points CD 
on the range curve. From the location of D it is seen that 1.49 = (m/)(Hp,/1704) =310 (m/u). This result gives u = 209 m 
for the mass of the meson as determined from the given observations. ; 


ber are most effective in distinguishing mesons For atoms to which the Thomas-Fermi model 
from electrons in the method of momentum loss, applies this formula gives 


but for these no direct determination of the 

2 = om 
mean excitation energy is available. Fortunately In(2me*/I) = 11.391 —InZ. (2) 
(a) the stopping power is not very sensitive to Let us now state the connection between 


the value of I and (b) Bloch’s theory of excitation range and energy or momentum in the form 
energies" permits a satisfactory estimate of this which is convenient for general use and at the 
quantity. According to Bloch, the mean excita- same time suited for determining the mass of 
tion energy for atoms containing many electrons _ the meson. The range of a particle of mass u and 
is proportional to the atomic number Z. For the charge number z may be considered to depend on 
constant of proportionality he gave the value any one of the following quantities: 

I/Z=0.96(me'/2h?) = 13.1 ev on the basis of the 


experiments available in 1933, but the data re- kinetic energy : E=c*(coshé— 1) 


ported in the preceding paper by Wilson lead to momentum: — b=ne sinhd 
a more reliable estimate of Bloch’s constant :!° magnetic rigidity: Hp=y(c?/ze) sinhé (3) 
velocity : v=c tanhé. 


I/Z=0.85(me*/2h?) = 11.5 ev. 
1 F, Bloch, Zeits. f. Physik 81, 363 (1933). Better than any of these quantities for inde- 
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as 


CHARGE OF PRIMARY 


MASS OF PRIMARY 
ELECTRONIC MASS 


RANGE NUMBER = 


TIMES Mic 


TIMES THICKNESS IN G /cm? 


Ht 
it 


mo 
PARTICLE (u/m) | 


H'+ 1836 938 


H7+ 3671 1875 


Hé++ 7295 3727 


) TIMES RANGE NUMBER 


8 


“a2 04 a6 as Lo 12 4 
(KINETIC ENERGY/ pe) = ENERGY/ 0.5109 Mev) 


Fic. 2. Calculated relation between range and energy for hydrogen and helium nuclei and for mesons. The point x 
corresponds to an alpha-particle with the following properties: energy =0.2 X 3727 = 745 Mev; range in aluminum =0,0046 


X (7295/4) X (26.97/13 X0.601 cm*/g) =29 g/cm? or 10.7 cm. 


pendent variables, however, are the dimension- 
less numbers sinh@ and cosh@, used in Figs. 1 
and 2, respectively. If the dependent variable or 
range is expressed in g/cm? it will fluctuate in an 
irregular way from element to element because 
of the lack of any precise relation between atomic 
number Z, and atomic weight, A. For this 
reason let us express the range in terms of 
a dimensionless “range number,’’ defined as 
follows: 


‘range number” = 42-Avogadro’s number 
- (e?/mc*)*-(Z/A)- (range in g/cm?) (4) 
= (0.601 cm?/g) -(Z/A)- (range in g/cm’). 


In the case of CH, the expression Z/A will 
have the value (6+4)/(12+4). In the case of 
dry air at 15°C, at a pressure of 76X13.6980 
dynes/cm?, division of the range by 2714 cm 
gives the range number. For lead of density 
11.3 g/cm’ range number equals thickness multi- 
plied by 2.68 

In terms of the “range number,” relation (1) 


takes the form (z*m/,) - (range number) = 


f d(cosh@+cosh~'@— 2) (s) 


In(2mc?/I) +2 In sinhé—tanh?6 


Numerical integration gives the curves” drawn 
in Figs. 1 and 2. The constant of integration was 
fixed in the case of air so as to give a range of 
238 cm for 15 Mev protons;? it was estimated in 
the case of the other substances listed. Negligible 
on the scale of the figures is the slight dependence 
of the constant of integration in (5) upon the 
charge of the primary, a dependence which 
arises principally from the phenomenon of cap- 
ture and loss of electrons at the very end of the 
range. Of course, the constant of integration has 
no effect on the determination of the mass of the 
meson by the method of momentum-loss. 

Figure 1 illustrates how measurements of 
magnetic rigidity before and after a meson 


® As long as they last, large photostats of these curves 
may be obtained from the Gcaumer. Department of 
Physics, Princeton University. We are indebted to Mr. 
Ralph Thompson for much help with the numerical com- 
putations. 
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TABLE I. Columns 1, 2 and 3 give measurements of Anderson and Neddermeyer on single cosmic-ray particles. On the pre- 
liminary assumption that each particle is a meson, the range-momentum relation allows as described in Fig. 1 a computation 
of the ratio (u/m) of unknown mass to electronic mass. Knowledge of the mass and momentum after traversal of the lead plate 
suffices to determine the factor 1/8? = (velocity of light/velocity of particle)? which is approximately proportional to the ioniza- 
tion. It is seen that this ionization would be exceedingly large for the last three particles if they were mesons. However, the ob- 
served rate of ionization, according to a kind personal communication of Anderson, is certainly not so great for these particles 
as twice the rate for a fast electron. Therefore the last three particles cannot be mesons of the computed masses. In the absence of 
information about the ionization of the other particles, we have calculated the ‘energy decrement number,” I, to test whether it 1s 
reasonable to assume that these particles are electrons which have experienced radiative losses. According to theory, all values of 
I between 0 and 1 are equally probable for electrons. Consistently small values of I are therefore an argument against attribu- 
ting electronic character to a group of particles. The angles @ elec. and @ meson are defined in the text. Arguments based on 
scattering indicate that no one of the last five tracks, even if it represents a meson, will permit a mass determination. 


1 2 3 4 5 6 7 8 9 10 11 12 13 

THICK- MOMENTUM ANALYSIS IF ELECTRON ANLAYSIS IF MESON 

NESS OF Cpe t ELec. cAp RANGE Ho 1 @ MESON 
LEAD (IN MEv) (MEv) t I (DEGREES)|} (MEV) NUMBER Hp: me B (DEGREES) 
1.35 113 86 102 3.1 0.00 33 27 3.6 1.32 150 1.8 44 
1.1 220 160 173 2.5 0.01 16 60 3.0 1.33 590 4.6 35 
‘2 200 125 138 2.5 0.02 21 75 3.0 1.60 560 6.2 52 
1.1 240 220 233 2.5 0.00 12 (20) 3.0 1.09 (310) 1.5 14 
13 38 6 19 2.5 0.07 430 32 3.0 6.3 [45 15.7 1700 
0.7 63 23 31 1.6 0.27 90 40 1.9 2.7 [110 6.9 238 
1.0 140 20 32 2.3 0.34 123 120 2.7 7.0 390 101 1230 
1.0 106 26 38 2.3 0.20 95 80 ar 4.1 250 26 475 
1.5 110 12 30 3.5 0.09 250 98 4.0 9.2 220 88 2350 


penetrates a plate of known thickness allow a_ the more reliable where the fractional loss in 
determination of the mass of the meson. momentum has been the greater. 

From the data alone it is not possible to tell 
DISTINCTION BETWEEN ELECTRONS AND MEsONsS_ which tracks represent mesons and which repre- 


‘ 
exist for the change in momentum in any given 


measurements on the momentum loss of single. 
: . , instance. Either the loss is due to a high rate of 
cosmic-ray particles which had traversed ~l 
ionization, in which case the particle must be 
of lead without undergoing multiplication. Their 
observations are reproduced in the first three moving slowly, and the quotient (momentum/ 
velocity) will be large and will indicate a meson; 
columns of Table I. As the cascade theory and 
or the loss is due to radiation, which is only 
, oye agree possible if the mass is comparable to that of an 
since emphasized, many particles of this char- h 
principle, however, three 
acter at sea level are mesons. On the preliminary 
qualitative means to distinguish between elec- 
i mesons, trons and mesons: (a) ionization, (b) scattering, 
(a) A meson emerging from the lead plate into 
quently the methods outlined in the caption of the cloud chamber will ionize more than an 
will of electron of the same momentum. The relative 
hi rate of ionization will be given approximately 
computed in column 12 of Table I on the assump- 
tion that the particles are mesons. The calcu- 
different from those to be expected for electrons, 


13. C. Anderson and S. H. Neddermeyer, International and - fact Dr. Anderson informs us that an 
Conference on Physics (London, 1934, Table II, p. 179). inspection of the last three tracks shows that 
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this criterion excludes the possibility that these 
tracks represent mesons. It will be noted that 
this test requires only a qualitative knowledge of 
the ionization, in contrast to those methods of 
determining the meson mass which depend on an 
accurate count of numbers of ions. 

(b) Particles of the same momentum will 
suffer the more scattering the greater is their 
mass. Following the treatment of multiple scat- 
tering given by Williams," we may express the 
arithmetic mean projected angle of deflection due 
to multiple scattering in lead approximately in 
the form 


6 in degrees = (2500 degrees Mev/pc). (6) 
(thickness in cm)!(1+[c?/pe 


Here there is some ambiguity as to the proper 
mean value to insert for the momentum, ?, 
which changes in traversing the plate. The 
scattering is greatest when the particle is travel- 
ing most slowly. We have used for p in (6) the 
momentum after traversal of the plate in order 
to obtain a rough estimate of the order of magni- 
tude of the scattering to be expected for each of 
the tracks in Table I. As columns 7 and 13 of 
this table show, there is not sufficient difference 
in the expected scattering for the first four 
particles to distinguish between electrons =m 
in Eq. (6) ] and mesons (u from column 11 of 
table). For the last five tracks the computed 
scattering is so great that the calculation itself 
has no meaning. It is clear that these tracks will 
give no information about the mass of mesons. 

Generalizing the foregoing conclusion, we may 
say that it is important to avoid large scattering 
in determining the mass of the meson by the 
method of momentum loss. This condition will 
be satisfied if the momentum of the particle 
after traversal of the plate is not too small. On 
the other hand, the momentum afterward must 
be considerably less than that before traversal if 
the measurement of the momentum change is to 
be accurate. Both conditions will be abundantly 
fulfilled if the primary momentum is of the 
order of 100 Mev and the momentum decreases 
in the plate by a factor of the order of 2. 

(c) In addition to ionization and scattering as 
criteria in selecting acceptable tracks, we have in 


4 E. J. Williams, Phys. Rev. 55, 303 (1940). 


the theory of radiation a means to analyze the 
possibility that the particles in Table I are fast 
electrons. If this is the case, they will have lost 
energy by ionization at a rate of ~12 Mev per cm 
of lead. We can therefore correct the values, co, 


_ of the final energy to obtain the value, E,, which 


would have been observed if radiative losses 
alone occurred. Column 4 of the table gives FE, 
and column 5 shows the thickness, ¢, of the plate 
as measured in terms of the so-called ‘“unit- 
length”’ of the theory of radiation. There is to 
be expected a statistical connection between E, 
and ¢. According to Bethe and Heitler'® the 
chance, dN, that an electron emerges from the 
plate with a final energy EF, (corrected for 
ionization losses) is given by the expression 


(7) 


Let us measure the loss of energy through radia- 
tion by an “energy decrement number,” J, 
defined as follows: 


In 
I= f y+ exp(—y)dy/(t—1)! 
=I(InE/E,,t-—1). (8) 


This quantity will be zero if there is no energy 
loss, unity if the particle loses all of its energy 
by radiation. It is a purely experimental quantity 
and is easily computed, for extensive tables of 
the function J(z, p) are available.'® In terms of 
the “energy decrement number,” relation (6) 
takes a very simple form: 


dN=dl. (9) 


In words, all values of the energy decrement 
number between zero and unity are equally probable, 
whatever be the original energy of the electron 
or the thickness of the plate which it penetrates. 

The energy decrement number, J, for each 
particle in Table I has been computed according 
to the definition in Eq. (7). The entries in column 
6 show that the values of J are on the whole 
quite small, particularly for the first four tracks. 
In the absence of other information, consistently 
small values of J will be an argument against 
attributing electronic character to a group of 
particles. 


%H. Bethe and W. Heitler, Proc. Roy. Soc. A146, 
83 (1934). 

*K. Pearson, Tables of the Incomplete Y-Function 
(London, 1922). 
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TABLE II. Previously published data on the mass of the meson. The methods used for the determination are: curvature in a 
magnetic field and ionization density (c and 1), curvature and range (c and r), change of curvature in a known amount of matter 
(c and c), curvature and collision with an electron (c and e). The errors given for the momentum and the mass are those cal- 

culated by the authors. The references are listed separately. Notes give additional information from original paper. 


cp (IN MEV) MAss (yu /m) METHOD REF. cp (IN MEv) MASss (yu /m) METHOD REF. 
1 33 220+50 candi 1 | 13 7.7 120+30 cc 9,3 
2 55 >430(< 800) candi 1 Note: Curvature of track changes in the gas by 0.6 
3 ; 44 190+60 candi cm/cm. Observed length of 5 cm sets a lower limit for the 
a 34 160+30 candi range which gives according to (3) an upper limit of mass 


Note: Values of 1, 3 and 4 are considered reliable; as 
upper limit of u/m for particle 2 is set by a minimum range 
equal to the observed length of track in the chamber, the 
lower limit is estimated from the lower limit of ionization. 


5 45 250 candi~ 2,3 
Note: Mass value considered as ‘‘good.”’ 
6 29 ~130 candi 4 


Note: Value only approximate. Corson and Brode, (3), 
calculate u/m = 160 with a different method for evaluating 
the ionization loss. 


as 110 m. 


222+3-+150+3 180+20 cc 10 

Note: Particle loses momentum by penetrating 4.8 cm 
of Pb; calculated by use of Bhabha’s relativistic formula 
[Proc. Roy. Soc. A164, 255 (1938) ]. 


15 11.6+0.2 170+8 candr 10 

Note: Radius of curvature is an average over range of 
6.5 cm, therefore calculated mass is somewhat doubtful, 
in spite of small margin of error given. No allowance made 
for the scattering which is important at such a low velocity. 


7 57 ~200 candi 5 | 16 — 2.4 20? candr 11 
Note: Value “probable,”’ but not very good. r 11 
1 120+30 candr 11 

8 16.5 <200 candr 14,3 | 49 442 55435 candr 11 
Note: Observable range was 18 cm, track out of focus, 20 1246 170+ 100 candr 11 


so that ionization density not measurable; if it was 10 
times normal, mass would be 125. 


9 16.5 ~350 candr 6,3 
Note: Observed by (6), calculated by (3); according to 
private communication from Anderson curvature is aver- 
age value over whole length of track (4 cm), particle be- 
sides probably scattered, computed value has ‘‘not much 
meaning. 


10 =42 < 1000 candr 6,3 

Note: The same critical remark a fortiori, as the particle 

‘does not stop in the chamber (range >5 cm) and the Hp 
value given is only an upper limit. 


11 52 220+35 candr 7,8 
Note: Considered by Anderson as reliable measurement. 


12 180 <65 cande 8 

Note: Value calculated from possible collision with an 
atomic electron which coming out of the plate shows an 
energy of 16 Mev; other interpretation possible. 


Note: Values obtained by use of ‘‘slow’’ cloud chamber; 
masses of particles 17 and 18 considered as better values 
than those of the other particles. No allowance made for 
the scattering which is important at such low velocities. 


21 93 240+20 cande 12 
Note: Value reliable. 
22 221126 250+50 cc 13 
23 165— 87 170+20 cc 13 
Note: 23 considered reliable. 
24 30 180+25 cand?; 
cande 15 


Note: Reliable. 


Average of reliable values: 180 


1E. J. Williams and E. Pickup, Nature oo. 648 (1938). 
2D. R. Corson and R. S Brode, Phys. Rev. 53, 215 (1938). 
3D. R. . Brode, Phys. Rev. 53, 773 (1938). 
4J. C. Street and E. c yee Phys. Rev. 52, 1003 (1937). 
5 P, Ehrenfest, Jr., C. R. Paris, 206, 428 (1938). 
6 C. D. Anderson and S. H. Neddermeyer, Phys. Rev. 50, 263 (1936), 
and 13). 
H. Neddermeyer and C. D. Anderson, Phys. Rev. 54, 88 (1938). 
S. H. Neddermeyer and C. D. Anderson, Rev. Mod. Phys. 11, 191 


(1938). (Fig. 19). 


MAss OF THE MESON 


The criteria of ionization and scattering ex- 
clude the last five tracks in Table I, and the 
change of momentum is too small for the fourth 
track to permit a reliable mass determination. 
Without a qualitative knowledge of the ioniza- 
tion we cannot exclude the otherwise unlikely 
possibility that the first three tracks represent 
electrons. If they represent mesons, however, 


*A. J. Ruhlig and H. R. Crane, Phys. Rev. 53, 266 (1938). 
asap. N. Nishina, M. Takeuchi and T. Ichimiya, Phys. Rev. 55, 585 

"H. Maier-Leibnitz, Zeits. f. Physik 112, 569 (1939). 

12 L. Leprince Ringuet, S. Gorovdetzky, E. Nageotte and R. Richard- 
Foy, Phys. Rev. 59, 460 (1941). 

183 J. G. Wilson, Proc. Roy. Soc. A172, 521 (1939). 

4 R. B. Brode and M. A. Starr, Phys. Rev. 53, 3 (1938). 

% Donald J. Hughes, Phys. Rev. 60, 414 (1941). 


they by no means lead to a unique value for the 
mass of the meson. Nor do the mass determina- 
tions previously published and compiled in Table 
II allow a decision of the very important question 
whether the mass of the meson is unique. In this 
connection we may recall the remark of Anderson 
and Neddermeyer in 1939? that ‘it has become 
increasingly likely that a complete interpretation 
of the experimental data is not to be found in 
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the single assumption of unstable particles with 
unit charge and a unique mass of the order of 
200 electron masses.’’ Also Weisz'’ has pointed 
out that assumption of a suitable distribution of 
meson masses allows a reasonable correlation of 
otherwise discordant estimates of the lifetime of 
the meson. 

Before any statistically reliable conclusions 
about the uniqueness of mesons are possible, it 
will be necessary to obtain many more measure- 
ments of the masses of penetrating particles. 
For this purpose the method of momentum loss 
appears especially satisfactory in view of the 
considerations presented in this paper. 

We are indebted to Dr. Carl D. Anderson for 
correspondence and to Dr. Marcel Schein for a 
number of interesting discussions on the proper- 
ties of the meson. 


APPENDIX. REMARKS ON EQ. (1) FOR STOPPING POWER 


The non-relativistic form of relation (1) was first derived 
by H. Bethe,'* who has also given the above relativistic 
formula. Equation (1) does not, however, agree with the 
relativistic expressions for stopping power as published in 
various general references.'® In the derivation one classifies 
the various energy losses, Q, into two groups according as 
they are less than or greater than a certain quantity, Qi, 
which is small in comparison to mc*, large in comparison 
with atomic binding energies, but which is otherwise 
arbitrary. The rate of loss of energy by collisions in the 
first group is 


(2 Ne*Z?/mv*) {In [2mv*Q,/I?(1 —0*/c*) ]— (v*/c*) } 


as Mgller®® and Bethe*® have both proven. E. J. Williams* 
obtained the same result by the method of impact parame- 


17 P. Weisz, Phys. Rev. 59, 845 (1941). 

18H. Bethe, Ann. d. Physik 5, 325 (1930); Handbuch 
der Physik (1933), second edition, Vol. 24, Part I, p. 523, 
Eq. (56.16). 

18 Mott and Massey, The Theory of Atomic Collisions 
(Oxford, 1933), p. 269; Heitler, The Quantum Theory of 
Radiation (Oxford, 1936), p. 218; Livingston ‘and Bethe, 
Rev. Mod. Phys. 9, 263 (1937), Eq. (750); Neddermeyer 
and Anderson, Rev. Mod. Phys. 11, 199 (1939). 

20 Mgller, Ann. d. Physik 14, 579 (1932); Bethe, Zeits. 
f. eT 76, 293 (1932). 

21 E. J. Williams, Proc. Roy. Soc. A139, 175 (1933). 


ters and was able to show that the relativistic terms in this 
expression arise from the Lorentz contraction of the field of 
the primary particle in distant collisions, and should there- 
fore be considered as especially reliable. In collisions of the 
second group where the energy transfer is large (Q>Q,) and 
lies between Q and Q+dQ, the calculated average rate of 
loss of energy is 


(24 max) +correction 


Here Q max is the maximum 
loss which a primary of momentum fo and total energy Eo 
can experience in a single encounter. The calculated value 
of the correction term is in general very large and depends 
in a striking way upon the spin and magnetic moment 
assumed for the heavy particle.” The calculated correction 
term is however negligible if the energy of the struck 
electron in the center-of-gravity frame of reference is small 
in comparison with mc. This condition is satisfied if the 
energy of the primary in the laboratory frame of reference 
is small in comparison with (u/m)uc*, or consequently 
Q max is considerably less than Eo. Then impacts for which 
Q is between Q, and Q max result in an average rate of loss 
of energy, 


(24 Ne*Z?/mv*) {In(Q max/Q,) — (v*/c*) |. 


If in addition uw is much greater than m, then Q max will be 
given by 2mpe?/u? = 2mzv?/(1—v*/c*), and the total average 
rate of loss of energy will be represented by expression (1). 
Relation (1) is valid only if the primary is considerably 
heavier than an electron; only if its energy is small in 
comparison with (u/m)yc*; only if it is moving faster than 
the bound electrons of the stopping material; only if 
capture and loss of electrons by it can be neglected and its 
charge is not too great ;* only if the energy is not so great 
that the density of the material has an influence upon its 
stopping power ;** only if the average energy loss and most 
probable energy loss are essentially identical.** All these 
conditions are satisfied in the present application. 


2H. J. Bhabha, Proc. Roy. Soc. Al64, 257 (1938); 
H. S. W. Massey and H. C. Corben, Proc. Camb. Phil. 
Soc. 35, 463 (1939); H. C. Corben and J. Schwinger, 
Phys. Rev. 58, 953 (1941); W. Pauli, Rev. Mod. Phys. 13, 
230 (1941), Table II. 

%M. S. Livingston and H. A. Bethe, Rev. Mod. Phys. 
9, 265 (1937), Eq. (758). 

*F. Bloch, Ann. d. Physik 16, 285 (1933); N. Bohr, 
Phys. Rev. 59, 270 (1941). 

% E. Fermi, Phys. Rev. 57, 485 (1940). 

2% N. Bohr, Phil. Mag. 30, 581 (1915); C. Moller, Ann. d. 
Physik 14, 531 (1932). 
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The present work is a study of the behavior of potassium, rubidium, and cesium ions emitted 
from a solid block of Kunsman catalyst, scattered at an angle of 90° in mercury vapor and 
collected on a stationary cylinder. The values of the incident and scattered currents were 
measured simultaneously. The incident beam used was of the order of several times 107% 
ampere, while the scattered beam was of the order of 10~* to 108 part of that value. Scattering 
was found to increase linearly with the pressure of the scattering vapor, to decrease with 
increasing energy of the incident ions for all ions studied, and to decrease with the mass of the 


incident ions. 


HE classic experiments of Rutherford on the 

bombarding of gold foil by alpha-particles 
and his theory of a nuclear atom! provided a 
working hypothesis for the study of the scattering 
of charged particles by atoms. Since those 
beginnings, impacts between particles of small 
mass and gaseous atoms have been studied by 
others using in their theories mechanisms of 
electrostatic repulsion, short range attractive 
forces, and polarization of the particles. These 
investigations also aimed toward securing a 
possible explanation of results from classic or 
wave mechanical laws for a rigid sphere atomic 
model. Heavier particles, such as positive ions of 
the alkali metals, also have been used for the 
study of scattering laws. Scattering at very small 
angles was investigated by C. A. Frische®? and 
A. G. Emslie,* angular distribution and energy 
loss of scattered positive ions by A. G. Rouse.‘ 
Scattering at a single angle obviates some of the 
experimental difficulties by simplification of the 
apparatus and elimination of moving parts. A 
scattering angle of 90° has several advantages. 
For, using a suitable set of slits at right angles to 
the path of the incident ion beam, it is possible to 
register practically all the ions scattered at this 
angle; hence the pressure of the scattering gas 
may be kept lower than would be feasible with a 
smaller movable collector, or alternatively, at the 


* Now at Mary Manse College, Toledo, Ohio. 

1 E. Rutherford, J. Chadwick, and C. D. Ellis, Radiations 
from Radioactive Substances (Cambridge University Press, 
1930), pp. 191-280. 

? C. A. Frische, Phys. Rev. 43, 160-168 (1933). 

3A. G. Emslie, Phys. Rev. 53, 680 (1938). 

“A. G. Rouse, Phys. Rev. 52, 1238-1244 (1937). 


same gas pressure a greater scattered beam 
intensity may be obtained. Furthermore, by 
proper arrangement and control of slits, it is 


Fic. 1. Positive ion scattering apparatus. F is the 
filament; K the catalyst; S:, S2, and S; the collimating 
slits; A; a nickel cylinder surrounding the catalyst; As a 
nickel collar; 71, JT: tubes to gauge and pumps; 55; slits for 
path of scattered beam; D the collecting disk; O Pyrex 
cylinder; M,, M: electrical leads; W Apiezon “W”; V mer- 
cury well; Ci, C2, C3, Cy, nickel cylinders with slits; C; the 
collecting cylinder; G the lead to grid of amplifier tube. S; 
was made from 1 to 15 volts negative with respect to the 
filament, while S; was made from 3 to 20 volts more 
negative than S2. C2: was kept 2.5 volts negative with 
respect to ground in order to prevent electrons, produced 
by ionization, from reaching the collector. It has been 
found that such electrons have little energy. (Rouse, 
reference 4, p. 1242.) 
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possible to prevent electrons, whether scattered 
or liberated in the scattering region, from reach- 
ing the collector. 


APPARATUS. PROCEDURE 


With these advantages in view, an investi- 
gation of the scattering at 90° of positive alkali 
ions in mercury vapor was carried out for 
potassium, rubidium, and cesium. The source of 
ions was a block of Kunsman catalyst,’ K in Fig. 
1, heated by a tungsten flament F and mounted 
above a set of three circular collimating slits by 
which the ions were accelerated and directed 
down along the common axis of five coaxial nickel 
cylinders where mercury vapor was maintained 
at a known pressure. The slit Sy, 1/16 in. wide, 
extended almost the entire distance around the 
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Fic. 2. Scattering of potassium ions. Ordinate is ratio of 
(current to collector cylinder)/(current to elliptic disk). 
The unit, 5.7 10-4, is the ratio of the sensitivity of the 
amplifier to that of the galvanometer. Abscissae represent 
the total accelerating voltage, i.e., between the filament and 
the third slit. 7; = 76.5° to 77.5°C, T2=65° to 66°C. 


four inner cylinders, Ci, Cs, Cs, Cs, a short 
distance from the upper end. The outermost 
cylinder C;, the collector for the scattered ions, 
was connected by a well-insulated and well- 
shielded lead wire to the grid of a D-96475 elec- 
trometer tube in a Penick amplifier circuit® of 
sensitivity 5.7X10-'5 amp./mm. The mercury 
vapor, the scatterer, was supplied from a well in a 


5C. H. Kunsman, Science 62, 269-270 (1925). 
°D. B. Penick, Rev. Sci. Inst. 6, 115-120 (1935). 
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Fic. 3. Scattering of rubidium ions. See Fig. 2 for ordinate 
and abscissa. 7; =90.5° to 93.5°C, T2=78.5° to 79.5°C. 


side tube connected to the scattering chamber. 
The temperature of the well was used to measure 
the mercury vapor pressure within the cylinders 
where the temperature was maintained at a 
slightly higher value to prevent excessive conden- 
sation of the vapor. D, an elliptic disk placed in 
the path of the incident beam to collect the 
incident ions, was set with its normal at an angle 
of 45° to the beam to decrease the possibility of 
reflected ions reentering the incident beam, and 
was connected to a galvanometer of sensitivity 
10-" amp./mm. The incident current to D and 
the scattered current to the outermost nickel 
cylinder were measured simultaneously, thus 
making unimportant any fluctuations in ionic 
emission. The solid angle subtended by the 
collecting cylinder was 1/28 of the total solid 
angle and made possible the collection of practi- 
cally all the ions scattered at 90°. 


RESULTs. DISCUSSION 


Figures 2-4 are the results of a study made of 
the scattering ratio (the scattered current/inci- 
dent ion current) in relation to increase of the 
total voltage accelerating the ions. The voltage 
plotted on the curves is the net accelerating 
voltage for the incident ions and was measured 
between the tungsten filament which held the 
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Fic. 4. Scattering of cesium ions. See Fig. 2 for ordinate and abscissa. 7; =80.5° to 81°C, T:=60° to 62.3°C. 


catalyst and the earthed third slit. The data were 
repeatable, and averages of several runs for the 
same range of temperature have been used in 
these curves. It is noted that 90° scattering 
always decreased with increase of accelerating 
voltage for all pressures used. To investigate the 
possibility of multiple scattering, measurements 
were made of variation of ratio of scattering with 
increase of pressure of the scattering mercury 
vapor. Results (Fig. 5) show that for the three 
alkali ions studied the scattering ratio increased 
linearly with the pressure of the scattering vapor, 
thus giving evidence that no onset of multiple 
scattering occurred. 

The possibility of scattering other than that 
caused by mercury vapor from the supply well 
was tested. Evidence for its occurrence was the 
existence of a considerable current to the collector 
cylinder when the scattering chamber was heated 
to 70 or 80°C while the mercury well was kept at 
room temperature or lower, whereas this was not 
found when both were at room temperature; but 
if during a temperature run the chamber temper- 
ature was maintained only a few degrees above 
the mercury well temperature, the ratio of 
scattering increased linearly with the vapor pres- 
sure as mentioned above. To investigate the 
cause of this condition, pumping was continued 
for 30 hours with the supply well packed in ice at 
0°C and the scattering chamber kept at 90°C or 


above except for one period of several hours 
during the night when the temperature was 
allowed to fall to 70°C. Scattering ratio readings 
taken at intervals during this test gave results as 
shown in Fig. 6, where it is seen that, for 100-volt 
rubidium ions, in 30 hours the ratio fell from 
10.5X10~ to 1.88X10-. This suggests a supply 
of scattering vapor other than the supply well, 
probably mercury vapor which could not be 
removed easily from side tubes, but as a back- 
ground the residual current was not large enough to 
cause any serious masking of the scattering results. 
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Fic. 5. Representative vapor pressure curves for two 
types of ions studied. See Fig. 2 for ordinate. Vapor pressure 
plotted is obtained from temperature as read on a thermom- 
eter mounted beside the mercury well in an oil bath sur- 
rounded by heating coils. 
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RATIO OF SCATTERING 


26 


HOURS OF HEATING CYLINDERS 


Fic. 6. Scattering of rubidium ions in relation to time of heating. See Fig. 2 for ordinate. During the time recorded 
on this graph the mercury well was packed in ice while the cylinders were heated to 90°C during the day and 70°C during 


the night. 


Other sources of possible errors in incident 
current readings may have been caused by: 
(1) reflection of ions from the walls of the inner- 
most cylinder or from the disk itself; (2) ioniza- 
tion products such as electrons and mercury ions; 
(3) scattering throughout the length of the beam; 
and (4) uncertainties in scale readings. The 
scattered current values would be less likely to be 
affected by the first two of these sources of error 
since the slit system would minimize errors due to 
reflection, the negative potential on the second 
cylinder retards electrons, and the probability of 
mercury ions reaching the collector cylinder is 
small except for those formed beyond the last 
slit. The probability of scattering from the 
scattered beam is about equal to that of scattering 
from the incident beam since the paths are 
comparable. Uncertainties in readings were 
chiefly those of the occasionally unstable zero of 
the galvanometer which was connected to the 
disk and incomplete knowledge of the extent of 
temperature equilibrium. The former was of im- 
portance only when the incident beam dropped 
to small values where it may have caused an 
error as high as 10 percent. The latter did not 
affect the data for the voltage curves as the 
temperature varied little during the runs. 


CONCLUSION 


The present experiments show greater 90° 
scattering for low energy ions. There is not 
complete agreement between this result and that 
of some former investigations of scattering with 
angle in a different type of apparatus, where it 
was found that scattering did not always decrease 
with increasing voltage but in some cases showed 
an increase,‘ whereas here it decreased for the 
three ions studied. These results are not entirely 
compatible with rigid sphere scattering, which 
does not depend on the kinetic energy of the 
impinging particle, but lie somewhat between 
this and Rutherford scattering, and in general 
suggest the action of other forces, such as 
polarization of the particles. There is noted also 
a decrease in scattering with increase of mass of 
the incident ions. 

In conclusion, the author wishes to acknowledge 
gratefully the guidance and valuable suggestions 
of Doctor Arthur G. Rouse, under whose direc- 
tion this work was done; also courtesies extended 
by the Reverend James I. Shannon, S.J., director 
of the Department of Physics, Professor A. H. 
Weber, and other members of the staff and 
department. 
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The 1940 data of Drinkwater, Richardson and Williams 
(DRW) on the He and Da lines have been completely 
recalculated, to get values of the Rydberg constant for 
hydrogen and deuterium, and from these, for helium and 
for infinite mass. Explicit equations, some of which have 
not previously been published or used, are derived for all 
the necessary calculations. The value of R.. now obtained 
is appreciably higher than that published by DRW, but it 
is surprisingly low, compared to the value commonly used 
for the past decade, which is based on the 1927 work of 
Houston on H and He’. His work, as well as later work by 
Chu, is also completely recalculated, with the use of the 
latest values of all auxiliary constants. About one-half of 
the proposed change in the value of R, is due to new values 
of the index of refraction of air. Recent work on this 
subject is summarized and put in a form suitable for use by 


spectroscopists. The data of DRW, as well as that of 
Houston and Chu, lead also to values of E (atomic weight 
of the electron) and e/m. The explicit equations needed for 
such calculations are derived. All recent precision work on 
e/m is also recalculated. Thus one gets finally twelve 
different precision values of e/m, obtained by seven 
distinctly different experimental methods. Although the 
data are not as consistent as could be desired, there are no 
really serious discrepancies. If the adopted value of the 
Faraday were raised by about 30 coulombs, almost complete 
consistency of the e/m values would be obtained. There is, 
however, no direct experimental evidence to justify such a 
change. The recommended values of R, e/m and E are listed 
at the end of the paper. That for e/m is (1.7592 +0.0005) 
X10? abs. e.m.u., as contrasted with the value 1.7591 
+0.0003 recommended three years ago. 


INTRODUCTION 


HE simplest known type of atom is one that 
consists of a nucleus and a single external 
electron. Examples of this type are neutral 
hydrogen, neutral deuterium, and singly ionized 
helium. Since such an atom represents a two- 
body problem, it is possible to obtain a complete 
solution for the frequencies of the resulting 
spectrum. The solution involves the important 
constants R, (the Rydberg constant for infinite 
mass), @ (the fine-structure constant) and E (the 
atomic weight of an electron). Hence, from 
sufficiently precise and detailed frequency meas- 
urements in the spectrum of H, D and Het, one 
can evaluate these three constants. 

Toobtain a value of a, certain measurements in 
any one of the three spectra are sufficient, whereas 
to obtain R, and E, measurements in any two 
spectra are required. In this latter case the only 
auxiliary constants needed are the two relevant 
isotopic masses (in atomic weight units). If, in 
addition, one assumes a value of the Faraday F, 
the value of e/m for an electron may be calculated. 
Furthermore, in the case of the spectra of H and 
D, only certain small frequency differences are 
needed to evaluate E or e/m; but in order to 
obtain R,, the actual frequencies must be 
measured. On the other hand, when He? is one 
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of the two spectra used, the actual frequencies 
must be known to obtain E, (and e/m), as well as 
to obtain R,. Finally, it should be noted that 
it is a wave-length in air that is ordinarily 
measured, and hence the value of the corre- 
sponding index of refraction of air is another 
necessary auxiliary constant. In fact the results 
of recent work on the index of refraction of air are 
chiefly responsible for the appreciable change 
now recommended for the hitherto accepted 
value of the Rydberg constant. 

The primary purpose of the present paper is to 
consider in some detail the numerical values of 
R,, and e/m given by recent experimental work in 
this field, as recalculated with a consistent set of 
auxiliary constants, and as directly as possible 
from the original data. The calculation of a 
involves so many additional considerations, and 
so much additional experimental data, that a 
separate paper is needed for any adequate presen- 
tation. Such discussion of @ as occurs in the 
present paper is therefore purely incidental to the 
calculation of R and of e/m. 

In the course of these calculations I have come 
across several small but not negligible errors, 
some in the theory and others purely arithmetic. 
Furthermore, the method of calculation em- 
ployed in the literature has often not been as 
direct as possible. A second purpose of the paper 
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is therefore to present explicit equations, in what 
appears to be the most convenient form, for the 
several desired quantities in terms of the experi- 
mentally determined quantities. Many of these 
equations, so far as I know, have not previously 
been published or employed. Aside from their 
convenience for calculation, such explicit equa- 
tions are needed in order to determine correctly, 
by propagation of errors, the probable error of 
the final result, as a function of the probable 
errors assigned to the various directly measured 
quantities. 

The original purpose of this work was to 
calculate a new value of R,. Such a calculation 
incidentally results in an evaluation of e/m, and 
in order to make the paper reasonably complete, 
I have accordingly examined several additional 
papers that yield a value of e/m (but not of R.), 
from measurements in the spectra of H, D and 
Het. The paper concludes with a brief discussion 
of all precision values of e/m that have been 
obtained in recent years, and a calculation of the 
present most reliable value of this important 
constant. 


THEORY 


In a two-particle atom such as neutral hydro- 
gen, the electron and nucleus each revolve 
around the common center of mass. If the mass 
of the nucleus were infinite, the common center of 
mass would coincide with the center of the 
nucleus. In this hypothetical case we are con- 
cerned with R,, the so-called Rydberg constant 
for infinite mass, for which Bohr derived the 
theoretical expression 27°e!m/h’c, where R, is 
expressed, as usual, as a wave number (cm™ 
units) and e is in e.s.u. The variation in the value 
of any actual Rydberg constant, such as Ry or 
Rp, from that of R,, is due solely to the failure of 
the center of mass of such an atom to coincide 
with the center of the nucleus. 

R, is now given in terms of the three specific 
Rydberg constants of interest to us by the 
equations 

R,=Ru(1+m/muy), 


R,=Rp(1+m/mp), (2) 
Re = Ru-(1 +m/mue) (3) 


where m= mass of electron, my = mass of hydro- 
gen nucleus, etc. Since only a ratio of masses is 


(1) 
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involved, we may choose any desired unit of 
mass. For convenience we accordingly express 
mass in atomic weight units, on the physical 
scale, and Eqs. (1)—(3) then become 


R.=Ru(1+£/H*), 
R.=Rv(1+E£/D*), (2’) 
(3’) 


where E=atomic weight of electron, H+ =H 
D+=D-—E, and He++=He-—2E, all on the 
physical scale of atomic weights. 

Any two of these last three equations are 
sufficient, in order to evaluate E and R,. 
Furthermore, knowing E and R, from any two 
equations, we can use the third equation to 
evaluate its own specific Rydberg constant. We 
can also express this third Rydberg constant 
explicitly in terms of the other two specific 
Rydberg constants. Such equations are given in 
this paper. 

We first present various explicit equations for 
E. Since we have 


(1’) 


e/m=F/E, (4) 


where F= Faraday, each explicit equation for E 

has a corresponding explicit equation for e/m. 

The unit of e/m depends on the system of units in 

which F is expressed, and we shall use e.m.u., 

since e/m is customarily given in such units. 
From Eggs. (1’) and (2’) 


H*+D*(Rp— Ru) 

D+Ru—H*Rp 

Similarly, from Eqs. (1’) and (3’) 
Het*+xRu—H*Rue 


(5) 


(6) 


Now the actual measured wave numbers (s), after 
proper correction to give the ‘‘Balmer line’’ wave 
numbers, are directly related to R. Thus for the 
H, and the D, lines, 


= (5/36)Ru, Sp= (5/36)Rp. (7) 
Therefore, putting Eq. (7) in Eq. (5), we get 
H+D*(sp— su) 
E= 


H*sp 


(8) 
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This last equation has been used by R. C. 
Williams! in his evaluation of E and e/m. A 
similar equation for He* and H is not of interest, 
since ‘‘corresponding” lines of He+ and H have 
not been precisely measured or, as a more direct 
statement, the lines of He+ and H that have been 
measured precisely are connected with their 
Rydberg constants by different proportionality 
factors. Hence Eq. (6) must be used, in 
evaluating E. 

Equations (5), (6) and (8) are not in the most 
convenient form for calculation, nor for showing 
clearly how the assumed errors in the experi- 
mental quantities affect the resulting error in E. 
To get more convenient equations, we combine 
Eqs. (1’) and (2’) in a way that employs the 
familiar relations, if A/B=C D, then 

A-B C-D A-B C-D 
= , or = : (9) 
B D A Cc 


We thus get, from Eqs. (1’) and (2’), in place of 
Eq. (5), 


H+D(Rp— Ru) 
= (10) 
Ru(D—H) 
or, using Eq. (7), we get 
(11) 
H) 


The proportional errors in D, H, D—H, and 
Su are much smaller than the error in the very 
small difference sp—su. Thus Eq. (11) shows 
that, to a close approximation, the proportional 
error in E equals merely the proportional error in 
Sp—Su. This equation also shows that, for the 
purpose of calculating E, it is not necessary to 
make a precise measurement of sp and of sy. Any 
approximate value from the literature may be 
used, and the entire experiment consists in 
measuring, as precisely as possible, the small 
difference sp —Su. The same statement applies to 
e/m, whose value follows from that of E, by 
merely assuming a value of F in Eq. (4). Thus, 
putting Eq. (4) in Eqs. (10) and (11) we obtain 


RuF(D—H) 
m= 
-H+D(Rp— Ru) 
1R. C. Williams, Phys. Rev. 54, 568 (1938). 


(12) 
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and 
suF(D—H) 


~ H+D(sp—su). 


(13) 


By combining Eqs. (1’) and (3’), using the 
relations of Eq. (9), we get, in place of Eq. (6), 


H*He*(Rue— Ru) 


(14) 
Ru(Het — H) 
The corresponding equation for e/m is 
Ru F(Het* — H) 
em= (15) 


H*+He*(Rue— Ru) 


This last equation is Eq. (8), p. 46 of G.C.1929,? 
while Eq. (14) is Eq. (7) of the same page. In 
that paper the value Of e/m was calculated from 
Houston’s data* on H and Het. In Eq. (15), the 
probable error in e/m is due almost entirely to the 
probable errors assigned to F and to the small 
difference, Rue— Ru. The errors in Ru, He+, H+ 
and He+—H are comparatively negligible. It 
should also be noted that in the case of Ht = H—E, 
He*+ = He-—E, etc., only an approximate value of 
E need be used. The precise value of E is then 
given by Eq. (14) or Eq. (11), ete. 

Just as Eqs. (1’) and (2’), or (1’) and (3’), have 
been solved simultaneously for E, so they can be 
solved for R,. But in the literature it has been 
standard practice to get R,, from Eq. (1’), with a 
value of E obtained by assuming values of F and 
e/m, in Eq. (4), i.e., 


(16) 

(e/m)H* 
I used Eq. (16) on p. 49 of G.C.1929, but one 
does not thus get the most reliable value of R., 
since the measurements on Het (or on D) are 
thereby ignored—unless one adopts just the 
value of E yielded by the experiment. The 
following equations for R,, which give its value 
directly in terms of the experimental quantities, 
have not, so far as I know, previously appeared 
in the literature. 


2 R. T. Birge, Rev. Mod. os 1, 1 (1929). This paper is 


always referred to as G.C.192 
3 W. V. Houston, Phys. = "30, 608 (1927). 
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From Eqs. (1’) and (2’), by eliminating E, we 
get 
RuRp(D—H) 
R 


o= (17) 
D*Ru—H*tRp 


This equation corresponds to Eq. (5) for EZ, and 
again by using Eq. (7) we get 


36 SuSp(D H) 


(18) 
5 D*sy H*sp 


_ which corresponds to Eq. (8) for EZ. Just as in the 
case of Eqs. (5) and (8), we can get a more con- 
venient form of Eqs. (17) and (18) by using the 
relations of Eq. (9). Thus from Eqs. (1’) and (2’) 
we get 

Ro— Ru 


R, Ruy = 
1—(RpH*/RyD*) 


(19) 


In order to evaluate R, we therefore need to 

calculate accurately only the small difference 

R,,— Ru, instead of the large factors of Eq. (17). 
In terms of wave numbers Eq. (19) becomes 


36 Sp— Su 


(20) 


Sat+ 
(spH*+/s_D*) 


Equations (19) and (20) thus represent more 
convenient forms, for accurate calculation, of 
Eqs. (17) and (18), and in fact Eqs. (19) and 
(20), when multiplied out, go directly into Eqs. 
(17) and (18). 

In a similar way, from Eqs. (1’) and (3’), we 


R,—R (21) 


which corresponds to Eq. (19). 
One can, of course, get from Eqs. (1’) and (3’) 
the less convenient form, corresponding to Eq. 

(17). It is 
RuRue(Het H) 


He++Ry—H*+Rue 


(22) 


This completes the equations for evaluating R.. 
But, as already noted, one can also express any 
one of the three Rydberg constants explicitly in 
terms of the other two. We shall give the result 
only for Ru. in terms of Ry and Rp, and then 
only in the form that uses the relations of Eq. 


(9). This result is 
D(Rp— Ru) (He* —H) 
Het(D—H) 


Rue— Ru= (23) 


or, in terms of wave numbers, 


Sp—SuH 
36 
. (24) 


D(Het—H) 


We have already noted how the error in E 
depends upon the errors in sy and sp. Let us now 
determine, from Eq. (20), the corresponding 
error in R,. We note in the first place that the 
denominator, lying under (sp—sy) as numerator, 
has a value close to 0.5, since spH*/sq”D* equals 
approximately 0.5. Hence the assigned probable 
errors in Sp and Sy, as they occur in this denomi- 
nator, produce a negligible effect on R,, as 
compared to these same errors, applied to the 
small difference sp—Su. 

For a similar reason the small probable errors 
in H+ and D* may be neglected. Thus for the 
purpose of calculating the probable error of R, 
one can, to a very close approximation, merely 
assume that the denominator lying under sp— su, 
in Eq. (20), is free from error, and has the value 
0.5. The equation then becomes 


(25) 


From the assigned errors in sp and sq, we now 
calculate by propagation of errors the resulting 
error in 2sp—Sy, and the proportional error in 
this quantity is seen to equal the proportional 
error in R,. 

In particular, if both sp and sy are assigned the 
same absolute error e, then the absolute error in 
2sp — Su is ey/5, and since sp~Su, Or 25p—Su~Su, 
the proportional error in 2sp—sy (and in R,) 
equals e(5)!/sy. In other words, the proportional 
error in R, is, to a close approximation,‘ \/5 

‘This method of calculating the probable error in 2, is 
essentially equivalent to writing out the rather complex 
expression for the probable error, as given rigorously by the 
law of propagation of errors [R. T. Birge, Am. J. Phys. 7, 
351 (1939) and then neglecting various terms the 
expression that are found to be negligibly small compared 
to the remaining terms. It must remembered, in this 


connection, that Eq. (25) is not sufficiently accurate for the 
calculation of a precision value of R... 
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times the assigned proportional error in sy 
(or in Sp). 

In a similar way, the denominator of Eq. (24) 
is very closely 3, and Eq. (24) is therefore 
approximately 


36 
(26) 


from which one can quickly calculate the 
probable error in Rye, resulting from the assigned 
probable errors of sp and sx. In particular, if both 
Sp and sy have the same assigned error, the 
proportional error in Rye is (5/2)! times the 
proportional error of either, and hence only (1/2)! 
as great as the resulting proportional error in R,. 

One can obviously express either R, or Rue in 
terms of its divergence from Rp, instead of from 
Ry. Thus Eq. (19), or its more convenient form, 
Eq. (20), may be replaced by 


Sp— 
36 | 
R,=— Dt (27) 
5 —1 
Sp:H* 


But Eqs. (20) and (27) are merely different forms 
of the same equation and hence necessarily yield 
identical values of R,, and also identical probable 
errors. 

The calculated values and calculated probable 
errors given in the following sections of the paper 
have been obtained by the use of the equations 
given in this section,—in particular, Eqs. (11), 
(13), (14), (15), (20), (21), (24), (25), and (26). 


INDEX OF REFRACTION OF AIR 


In G.C.1929 the adopted value of R, was 
109,737.42+0.06 cm-'!. This value, as already 
stated, was based on Houston’s data.* With the 
same data but with more recent values of the 
auxiliary constants, I obtained,’ in August 1939, 
the value 109,737.45+0.06 cm. Since then there 
has appeared an important paper by Drink- 


5 Asa preliminary to a new detailed article on the general 
constants, I compiled, in August 1939, a mimeographed list 
of constants, together with a brief explanation of the more 
important changes. This list will be called G.C.1939. As a 
result of the appearance of new data, as well as of numerous 
suggestions received since then, many changes have been 
made in the list. The present article gives the detailed 
explanation of the change in R,.. 


water, Richardson and Williams,* who obtained 
R,=109,737.264 cm~', with the auxiliary con- 
stants of G.C.1939. I have now recalculated their 
data and obtain 109,737.303+0.017 The 
causes of this revision appear in the detailed 
discussion to follow. It is, however, immediately 
obvious that the new value of R, deviates by a 
surprisingly large amount from the 1929 value. 
In fact most of the rather large probable error of 
+0.06 cm~ in the 1929 result arises from the 
necessary uncertainty of the value of the inter- 
national meter in terms of the red cadmium wave- 
length. This wave-length is 6438.4696 I.A., 
defining the I.A. system of wave-lengths, but it is 
(6438.4696+0.0020) X10-§ cm in the c.g.s. 
system.’ 

The probable error in the 1929 value of R., in 
the I.A. system, was believed to be only about 
0.02 cm~', as contrasted with the change of 
0.12 cm now suggested. One important cause 
for this large change, as DRW point out, is the 
new value of the index of refraction of air used by 
them. I shall therefore, before proceeding further 
in the calculation of R,, summarize the present 
experimental data on the dispersion of air. Such a 


6 J. W. Drinkwater, Sir Owen Richardson, and W. E. 
Williams, Proc. Roy. Soc. A174, 164 (1940). To be referred 
to as DRW. 

7The +0.0020 estimate of error is the value that is 
needed to give Houston's estimate’ of the absolute error in 
Ry and Rye, as contrasted with the error in their difference, 
used on p. 46 of G.C.1929. J. E. Sears and H. Barrell, 
Trans. Roy. Soc. A233, 143 (1934) give a list of four values 
of the wave-length of the red cadmium line, as measured 
directly in terms of the international standard meter. The 
unweighted average is (6438.4696+0.0004) x 10-8 cm. In a 
later report J. E. Sears, Science Progress 31, 209 (1936), 
adds two more recent measurements by Késters and 
Lampe, at the Reichsanstalt. The new unweighted average 
is 6438.4690 (Sears gives 6438.4693, and the origin of his 
average is not evident), but the probable error is still 
+0.0004. L. J. Briggs, Rev. Mod. Phys. 11, 111 (1939), 
reproduces Sears’ table, but without the original references. 
In spite of the remarkable agreement shown by the various 
measurements of the red cadmium line (the average 
deviation from the mean is only one part in seven million}, 
it seems necessary to adopt the larger uncertainty of 
+0.0020 x 10-* cm, because of the unsatisfactory character 
of the defining scratches on the international meter. 
W.E. Williams, Nature 135, 496 (1935), gives a photograph 
of one of the terminal scratches on substandard meter 
No. 26, magnified 300 diameters. The scratch, which is 
actually a series of rather irregular parallel furrows, is over 
two mm wide (i.e., 0.006 mm on the actual meter bar, or 6 
parts in 108 of its length). Thus the adopted probable error 
in the wave-length of the red cadmium line corresponds to 
an error of a 1/20 of the width of such a terminal 
scratch. The fact that the average deviation of the results 
of different observers corresponds to only 1/40 of the width 
is chiefly evidence of a uniformity of interpretation of the 
terminal scratches on the international meter bar. 
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summary seems especially advisable, because 
thus far no tabulated data have been published 
in the form needed in spectroscopy. 

Since 1918 all measured wave-lengths in 
spectroscopy have been reduced to vacuum by 
the use of Meggers and Peters’ formula’ for the 
dispersion of air. The well-known tables of 
Kayser® are based on their data. But since 1918 
the national physical laboratories of Germany, 
France and England have made new measure- 
ments, with almost identical results, but the new 
results differ considerably from those of Meggers 
and Peters. The German work'® has been pub- 
lished only in abstract form, and the paper giving 
the French work" is not readily available. 
Fortunately, however, all of the results are sum- 
marized in the report on the English work” 
although, as just noted, not explicitly in the form 
needed by spectroscopists. 

All standard spectroscopic wave-lengths are 
based on the wave-length of the red Cd line, in 
“normal” air, at 15°C, 760 mm. ‘“‘Normal” air is 
at present taken to mean air free from water 
vapor, but containing 0.03 percent of COs. Now 
all of the recent published results on the index of 
refraction of air refer to CO>-free air. To reduce 
them to “normal” air I have added 0.044 
to the value of u, in accordance with the data of 
Barrell and Sears."* Their published equation for 
the 20°C value of for CO>-free has been 
reduced to 15°C by the use of their own measured 
value of a=0.003674. Perard’s results" have been 
changed from 0°C to 15°C by means of his 
measured value of a=0.003716. Késters and 
Lampe’s results'® have been changed from 20° to 
15° by means of their quoted a=0.00367. The 
exact constitution of the air used by Meggers and 
Peters* was not measured, and it is now merely 
assumed to have been “normal” air, so that no 
correction for CO, content needs to be made. 

One further recent investigation on the disper- 
sion of air is that of Bender.'® His results differ 
from those of references 10, 11 and 12 even more 

8 W. F. Meggers and C. G. Peters, Bull. Bur. Stand. 14, 
697 (1918). 
om 3. Kayser, Tabelle der Schwingungssahlen (Leipzig, 

10 Késters and Lampe, Physik. Zeits. 35, 223 (1934). 

"A, Perard, Trav. Bur. int. Poids Mes. 19, 1 (1934). 

12H. Barrell and J. E. Sears, Phil. Trans. 238, 1 (1939). 

13 Reference 12. See p. 54. 


™ Reference 12. See p. 30. 
18D, Bender, Phys. Rev. 54, 179 (1938). 
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than do those of Meggers and Peters, and in the 
opposite sense. His equation for » at 0°C has been 
reduced to 15°C by the Barrell and Sears value of 
the observed a for air (0.003674), in place of the 
coefficient of a perfect gas (0.003661) used by 
Bender himself. 

Remembering that all published results, except 
those of Meggers and Peters, refer to CO>-free 
air, we then get, finally, for the dispersion of 


“normal” air at 15°C, 760 mm pressure, 
(28) 


where A is the wave-length in air, in microns, and 
the coefficients have the following values. 


a b c 
Reference 8 (M and P) 272.643 1.2288 0.03555 
Reference 10 (K and L) 272.742 1.5017 0.01834 
Reference 11 (Perard) 272.860 1.4002 0.02994 
Reference 12 (B and S) 272.581 1.5453 0.01268 
Reference 15 (Bender) 273.531 1.4946 0.03990 


One can get a proper comparison of these 
various results only by noting the resulting 
values of yw, at various parts of the visible 
spectrum. Such typical values of (u—1)X 10° are 
as follows. 


Reference 8 10 11 12 15 


6438.47 275.814 276.471 276412 276.383 277.369 
5085.82 277.925 278.822 278.711 278.745 279.906 
44358.32 280.097 281.156 281.061 281.068 282.505 


It will be noted that references 10, 11 and 12 
agree almost perfectly, whereas the results of 
reference 8 are low and those of reference 15 are 
high. In all of the calculations of this paper I have 
used the Barrell and Sears’ 1939 dispersion curve 
(reference 12), since their work is the most 
recent, and appears to be also the most thorough. 

In 1934 Sears and Barrell,’ in connection with 
their measurement of the red Cd line in terms of 
the standard meter, made a determination of u 
for air, for the red Cd line only. The value that 
they obtained, for “‘normal”’ air, is 276.494.'® 
This gives 1.7802A for the wave-length correction 
of the red Cd line to vacuum. Hence if, in 
agreement with Sears and Barrell,’ we adopt 
6438.4696 X10-* cm as the best value of the 
wave-length in “normal” air, the best value in 


16 Sears and Barrell give u only for CO,-free air, but their 
published values of \6440.25099 in vacuum and 6438.4708 
in “‘normal” air correspond to this value of (u—1) 10°, as 
contrasted with their 1939 value of 276.383, as already 
given. 
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vacuum is 6440.2498 X 10-8 cm on the basis of 
their 1934 work. On the other hand, the 1939 
work of Barrell and Sears® gives 1.7795A as the 
correction, and hence 6440.2491 X 10-8 cm as the 
wave-length in vacuum. This 1939 work evidently 
appeared just as the DRW paper‘® was ready for 
the press. Hence all of the published DRW 
results are based on the 1934 value of \6440.2498 
as the vacuum wave-length, and to these they 
have added adjusted values corresponding to the 
1939 value of \6440.2491. In discussing their 
data, it is more convenient to leave it in its 
original form, and finally to make this small 
adjustment, just as DRW have done. It amounts 
to an increase of each wave number by 1.687 parts 
in 10’. This increase is 0.001655 cm~! for the H, 
and D, lines. 


EXPERIMENTAL Data or DRW, 
AND CALCULATIONS 


In order to obtain a value of R,, DRW® 
measured in vacuum with a reflection echelon the 
wave-lengths of the H, and D, lines. This 
instrument consisted of 40 plates, each 6.88 mm 
thick, and had a calculated resolving power of 
0.018 cm~', or of 846,000(=A/AA) at the wave- 
length used. The chief objection to this type of 
instrument is doubtless the limited field (not over 
two spectral orders) and the consequent radical 
warping of the true intensity distribution by the 
“envelope” of the diffraction pattern. However, 
when the theoretical correction had been made 
for this warping, the resulting corrected micro- 
photometer curves were not only extremely self- 
consistent, but they also satisfied the various 
known or expected intensity relations. Hence the 
final experimental results appear to be quite 
reliable. The wave-lengths in vacuum of the H, 
and D, lines were measured directly in terms of 
the red Cd line, so that, as already noted, only 
the vacuum wave-length of that line need be 
known. DRW publish only the resulting wave 
numbers in vacuum of the various lines measured. 

The theoretical pattern of the fine structure 
components of the H, (and D,) line and the 
energy level diagram from which the pattern is 
obtained, are given by DRW on page 165 of their 
paper. Each of the levels lies below the ‘‘Balmer”’ 
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level by an amount,’ 

j+} 4n 
where Ry should be used for Hz and Rp for Daz, 
with Z=1. In the case of He 4686 one uses Rue 


and Z=2. The “Balmer” lines are given by the 
simple formula 


AT (29) 


s=RZ?( —-— }. (30) 


The problem is thus to obtain in some way the 
wave number of some one fine structure com- 
ponent, and then to apply the necessary cor- 
rection to give the wave number of the Balmer 
line, which for H, and D, corresponds to n’=3, 
n'’=2. The component chosen by DRW® and by 
Houston* is number one, for which j’=5/2, 
j'' =3,/2. This, the strongest of the predicted five 
components, cannot be resolved from compo- 
nents 4 and 5, which should lie to the long wave- 
length side at distances of approximately 0.036 
and 0.144 cm~', with theoretical intensities 1/9 
and 1/45 of that of component 1. The “‘center of 
mass”’ of these three components accordingly lies 
0.00637 cm~' to the red of component 1, if 1/a@ is 
taken as 137. Houston*® used a correction of 
0.0056 cm~', and Williams and Gibbs'® pointed 
out that, on the foregoing simple assumptions, 
0.0063 cm~! should have been used. 

The problem is, however, much more involved, 
since each fine structure component has its own 
intensity distribution, due mainly to the Doppler 
effect. The half-width depends on experimental 
conditions, and must be determined in each case 
more or less empirically.'® But both Williams and 
Gibbs'* and DRW‘* agree that if the upper part of 
the intensity curve is used to locate its ‘‘center 
line,’”” component 5 has a negligible influence and 
this center line lies almost exactly 0.0036 cm! to 
the red of component 1. Incidentally this is just 


7See H. E. White, Introduction to Atomic Spectra 
(McGraw-Hill, 1934), p. 137. 
(934) C. Williams and R. C. Gibbs, Phys. Rev. 45, 491 

19 Thus R. C. Williams and R. C. Gibbs, Phys. Rev. 45, 
475 (1934) find half-widths of 0.180 cm™ and 0.171 cm™ 
for the components of Hg and Da, respectively. Their paper 
gives a detailed account of the method used by them in 
analyzing the microphotometer curves. Different investi- 
gators have used somewhat different methods. 
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the correction gotten on the simple “center of 
mass” theory, when the very weak component 5 
is ignored. 

In subsequent work R. C. Williams! found a 
half-width of 0.165 cm and a correction of 
0.0041 cm for the effect of component 4 of Hg. 
He then makes an additional correction of 
—0.0012 cm~! for the effect of component 3, 
which consists really of two coincident com- 
ponents, of total intensity about } that of 
component 1, and which lies 0.221 cm~ to the 
violet of component 1. DRW give only the 
position of their 1,4 complex and state that the 
effect of component 3 has already been taken into 
account in determining the position of this 
complex. Williams! uses slightly different cor- 
rections for H, and D, whereas DRW use the 
same correction for both. A general study of the 
evidence indicates that 0.0036 cm is probably 
the best correction to be applied to the data of 
both Houston* and DRW,,° and for both the H, 
and the D, lines. As a matter of fact, whether 
0.0036 or 0.0056 cm~ is chosen for the correction 
makes very little difference in the resulting value 
of the Rydberg constant but, as will be shown 
presently, Houston’s data are internally more 
consistent when 0.0036 cm is used. 

Houston also used component 1 of the Hg line 
to make an additional evaluation of Ry. The 
simple “center of mass’’ procedure here gives 
0.00269 cm~'! as the necessary correction to the 
1,4,5 complex. Williams and Gibbs'* erroneously 
publish this last result as 0.0029 cm~', but they 
find by the type of empirical analysis used for Ha 
that 0.0024 cm~ is the required correction. This 
latter figure happens to be just the correction 
used by Houston,* and I have accordingly 
adopted it. 

We now start the explicit calculation of the 
DRW data on H, and D,. The wave number of 
the 1,4 complex of H, was obtained from 10 
different plates, with an unweighted average of 
15,233.0654+0.0003. cm~', where the stated 
uncertainty is merely the standard probable 
error of a mean, and hence is based wholly on 
internal consistency of the data. To allow for 
possible systematic errors I have arbitrarily 
included an additional error of 0.0010 cm, 
which may add to or subtract from the preceding 
error. Hence the resulting uncertainty (square 


root of sum of squares of the errors) is +0.0010; 
cm. DRW base all uncertainties, throughout 
their paper, on the root-mean-square (r.m.s.) 
error of a single result (or single plate in this case), 
with no further allowance for errors of any kind. 
There seems to be no particular logic in such a 
procedure since according to it a series of obser- 
vations yields a result no more reliable than a 
single observation. In this case their published 
uncertainty is +0.0014 cm~. 

We now add 0.0036 cm to the wave number 
just calculated, and thus get 15,233.0690 cm= 
for the wave number of component 1. This value 
is now to be corrected to get the Balmer line. The 
necessary (negative) correction” is, by Eq. (29), 
Rya?(1/64—1/324) =0.07327 if 1/a=137. 
It seems best, wherever possible, to use the value 
of a determined by the data themselves but, as is 
well known, it has not been possible to get a very 
reliable value of a from the analysis of H, and 
D,. On the other hand, the magnitude of the 
correction is so small in this case that it is almost 
immaterial what value of a is used. As a matter 
of fact, these last two facts are connected logi- 
cally. Thus in the case of \4686 the correction is 
large and it is important what value of @ is used. 
But there the fine structure components are more 
widely spaced and it is possible to get a fairly 
reliable value of a from their analysis. 

With the Balmer line correction just given, we 
get 15,232.9957; cm~ for the Balmer line itself. 
This result, however, is in terms of the 1934 
index of refraction of air.’ It becomes 15,232.9974 
cm in terms of the 1939 index," and this last 
quantity I denote sq. Thus* 


$y = 15,232.9974+0.0003, 
or +0.0010; cm~. 


Multiplying sq by 7.2-(the reciprocal of 1/4—1/9 
in Eq. (30)) we get 


Ry = 109,677.581,+0.002, 
or +0.007; 


DRW give Ru = 109,677.583+0.010 


20 DRW write 0.072-0.0036 as the total correction to be 


subtracted, but give this difference as 0.0694, which happens 
to be nearly correct. (0.0733-0.0036 = 0.0697.) 

11 give the calculated results of DRW in general with 
two probable errors. The first is based merely on the internal 
consistency of the data, whereas the second (which is 
adopted) includes some allowance for systematic error. 
— only one probable error is stated, it is the latter of 
the two. 
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The D, line is now to be similarly treated, in 
order to evaluate Rp. The 1,4 complex of this line 
is measured by DRW on 7 plates, with an 
unweighted average of 15,237.2096+0.0003, 
cm~', The probable error happens to be identical 
with that for H, and I have again allowed 0.0010 
cm~' for systematic error. As already discussed, a 
correction of +0.0036 cm is again to be applied 
to get component 1, and the correction for the 
Balmer line is 0.07329 cm, the small change 
from that for Ha being due to the use of Rp in 
place of Ry. With the additional correction for 
the 1939 index of refraction, one thus gets finally 


Sp = 15,237.1416+0.00032 
or +0.0010; cm 
and 
Rp = 109,707.419;+0.0023; 
or +0.007; 
DRW give Rp = 109,707.421+0.009. 


We now calculate R, directly from the ob- 
served values of sy and sp, by the use of Eq. (20). 
For this purpose we need values of the atomic 
weights of ‘H and D, and for convenience I give 
at this point all the auxiliary constants that will 
be needed in the remainder of the work. They are 
(all on the physical scale) 


H = 1.00813+0.000017 
D =2.01473+0.00001., 
He = 4.00389 +0.00007, 
C =12.01465 +0.00023, 
F=9651.49+1.0 abs. e.m.u., 
E=0.0005486 (used in getting H*+, D* etc.), 
p=1.00048 +0.00002 (one int. ohm 
=p abs. ohm), 
qg=0.99986+0.00002 (one int. amp. 
=q abs. amp.), 
pq =1.00034+0.00003 (one int. volt 
= pq abs. volt), 
c= (2.99776+0.00004) cm/sec. 


The values of the atomic weights and of p are 
those used in a previous letter” on e/m, but the 
value of g has been decreased by a significant 
amount.” 


# R. T. Birge, Phys. Rev. 54, 972 (1938). 

23 The values of p and g are now in a very satisfactory 
state. An entire paper would be required to cover the 
details, but the values here adopted are obtained as follows. 
(1) H. L. Curtis, R. W. Curtis, and C. L. Critchfield, 
J. Research Nat. Bur. Stand. 22, 485 (1939) obtain 
q= 0.99986 in terms of the N.B.S. int. amp. This becomes 


Equation (20) then gives 


R, 109,737.302;+0.005. cm 
or +0.0163 


The method of calculating the probable error in 
R,, from the assigned probable errors in sy and 
Sp, has been discussed earlier in the paper. 
Equation (20) gives just the result that would be 
obtained by calculating e/m from Eq. (12), and 
then using this value of e/m in Eq. (16). This 
latter procedure was used by DRW, but they 
obtained 109,737.26,;+0.020 although they 
used the latest values of the auxiliary constants, 
as just listed. The discrepancy is due to two 
errors in their use of Eq. (16). The atomic weight 
of neutral H was used in place of the correct Ht, 
and the atomic weight of the mixture of hydrogen 
isotopes (H and D) was used in place of the 
correct isotope of mass 1. I noticed the latter 
error, when the DRW paper first appeared, and 
calculated corrected values, which DRW later 
published.** Unfortunately I did not notice the 


0.999868 in terms of the mean int. amp. (2) P. Vigoureux, 
N.P.L. Collected Researches 24, 173 (1938) gets ¢ =0.99986 
in terms of the N.P.L. int. amp. which becomes 0.99985 in 
terms of the mean int. amp. From these two results I get 
as a final average, g=0.99986+0.00002, in terms of the 
mean int. amp. (3) H. L. Curtis, C. Moon, and C. M. 
Sparks, J. Research Nat. Bur. Stand. 21, 375 (1938) obtain 
as the weighted average of this and a previous investigation, 

= 1.000468 in terms of the N.B.S. int. ohm, or 1.000472 
in terms of the mean int. ohm. (4) Two investigations at 
the N.P.L., one by L. Hartshorn and N. F. Astbury, Phil. 
Trans. A236, 423 (1937), and the other by P. Vigoureux, 
N.P.L. Collected Researches 24, 277 (1938) each give 
p=1.00050 in terms of the N.P.L. int. ohm, or 1.000508 in 
terms of the mean int. ohm. A slightly higher accuracy is 
claimed for the N.B.S. work, and I accordingly adopt as a 
final weighted average, p= 1.00048+0.00002, in terms of 
the mean int. ohm. 

The value of F given in reference 22 is 9651.31 +0.80 abs. 
e.m.u., on the physical scale. A slightly more precise 
calculation shows that this figure should have been given as 
9651.33 e.m.u. The present adopted value is only slightly 
larger, but actually two large and almost compensating 
changes have been made. In the first place the value of q 
has been changed from 0.99993 to 0.99986. This change 
decreases the value of F by 0.67 e.m.u. On the other hand, 
I have now adopted a certain weighted average of the 
values of F given by the silver and the iodine voltameters, 
in place of merely the silver result. This very involved and 
controversial subject will be discussed in detail elsewhere. 
The change increases the value of F by 0.73 e.m.u. A very 
small change (from 1.00027 to 1.000272) has also been 
made in the ratio of the physical and chemical scales of 
atomic weights. The net result of these various changes is 
an increase of 0.07 e.m.u. in the value of F. 

The atomic weights are the standard values given by 
M.S. Livingston and H. A. Bethe, Rev. Mod. Phys. 9, 245 
(1937), except that for carbon, which I have derived from 
the latest experimental results. 

24 DRW, Proc. Roy. Soc. A175, 345 (1940). 
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other and more serious error until I had started 
the systematic calculations of this paper.*® 

In addition to R,, the data of DRW allow the 
calculation of E and of e/m. The value of E may 
be obtained directly from Eq. (11), with the 
result 


E= (5.48647 +0.00059) 10-4 
or +0.00099. 


The error in E depends directly on that in 
(sp—Su) and any systematic error in this differ- 
ence is likely to be less than that in either sp or 
Sy. Hence the second quoted probable error in E 
is based on an assumed systematic error of only 
0.0006 cm! in sp — Su, in place of the 0.0010 cm~! 
systematic error already adopted for sp and su, 
separately. 

Then from Eq. (4), or directly from Eq. (12), 


e/m = (1.75913+0.00027) X 107 e.m.u. 
or +0.00037. 


DRW, using Eqs. (8) and (4), get E = (5.48646 
+0.00050) X10-*, and e/m=(1.7591+0.0004) 
X<107 e.m.u. For this last result they used 
F= 9651.31 abs. e.m.u. 

A discussion of the causes for the deviation of 
DRW’s value of R, from my 1929 value is 
postponed until Houston’s data* have been 
recalculated. But it should be noted at this time 
that there is a second method for calculating R, 
from data such as that of DRW. Instead of using 
Eq. (20), we may use Eq. (1’), or its equivalent, 
Eq. (16), with some ‘‘outside’’ value of e/m and 
of F (to give E of Eq. (4)). R,z may also be 
calculated from Eq. (2’), and Eqs. (1’) and (2’) 
will then yield different values of R, unless the 
adopted value of e/m happens to be identical 
with the value yielded by Eq. (12). In the case of 


*6 This confusion regarding the isotopes of hydrogen is a 
direct consequence of the unfortunate situation that exists 
in respect to the symbols in use. In writing equations such 
as those used in this paper, it is obvious that H and D are 
more convenient than H' and H?, since double superscripts 
are thereby avoided. On the other hand, the symbol H 
represents, to the chemist, the atomic weight of the mixture 
of hydrogen isotopes H' and H?, just as C represents the 
mixture of C” and C®. In fact, in setting up a partial table 
of atomic weights (on the physical and on the chemical 
scale), as I have done in G.C.1939, it seemed really 
necessary to use H' and H?®. For if H and D are used, how 
is the atomic weight of the mixture to be symbolized? Thus 
there are definite objections to each of the two sets of 


symbols. 


DRW, the value of e/m given by Eq. (12) does 
happen to satisfy this condition very closely.*® 

Even in this situation, the resulting probable 
error in R, will be different, according as Eq. (1’) 
or (2’) is used, and each error will be different 
from that resulting from the use of Eq. (20). 
Furthermore, with equally accurate values of Rp 
and Ruy, the former—used in Eq. (2’)—will lead 
to the more accurate value of R,. This necessary 
fact is shown clearly by the actual results. 
With the “outside” values of e/m=(1.7591 
+0.0005) X10? and F=9651.49+1.0, and with 
Ru = 109,677.5812+0.0075, we get from Eq. (16) 
—or from Eqs. (1’) and (4)— 


Rz = 109,737.303¢+0.019; cm—, 


whereas with Rp = 109,707.4193+0.0075 we get, 
from Eqs. (2’) and (4), 


109,737.3032+0.0113 em. 


The value of R,, obtained from Rp alone is thus 
apparently more accurate than that obtained by 
Eq. (20) from sx and sp simultaneously, even 
though the ‘outside’? value of e/m has been 
assigned a slightly larger error (+0.0005) than 
the calculated error (+0.00037) of the “‘inside”’ 
value, and the use of Eq. (2’) involves also the 
error in F, whereas Eq, (20) does not. It must 
be realized, however, that the use of Eq. (20), 
rather than of Eqs. (2’) and (4) with ‘“‘outside”’ 
values of e/m and F, eliminates the possibility of 
certain types of systematic error (in particular, 
errors in F), and hence I believe that the result 
yielded by Eq. (20), i.e., 


=109,737.303+40.017 


is probably the best to adopt. It is, in fact, the 
final result, as adopted in this paper, except that 
the probable error should be increased to +0.05 
cm~! to allow for the uncertainty in the relation 
of the I.A. to the c.g.s. system of length. 

Finally, from DRW’s values of Ry and Rp it is 
possible to calculate any other R, such as Rye. 
This last constant is given by Eq. (24), directly in 
terms of the original data sq and sp, with the 
result 


Rue = 109,722.263,;+0.012 


26 The present most probable value of e/m is considered 
later in this paper. 
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In addition to the work on R, and e/m, DRW® 
also make a detailed study of the fine structure of 
H, and D,. Such work should lead to a fairly 
precise value of the fine structure constant. 
However, as is well known, repeated investi- 
gations of this subject have shown rather defi- 
nitely that the position and intensity of certain 
components, notably component 3, are not in 
agreement with theory. Pasternack?’ has tried to 
explain the observed discrepancies by assuming a 
displacement of the 22S level by about 0.03 cm~. 
DRW, however, believe that a slight contami- 
nation of the fine structure components of D. 
by H, components, and of both D, and H. by 
molecular lines, is sufficient to account for all 
observed irregularities. It is just because of such 
irregularities, which have been known for many 
years, that so much work has been done on the 
fine structure of H, and D,. There is as yet no 
reason to suppose that the method we have just 
used in calculating Ru, Rp, and R, is affected by 
the observed irregularities, but because of them 
it has not been possible to get a reliable value of 
a, as already stated. Furthermore, because of the 
great amount of experimental work on the 
subject, a separate article is needed to cover it 
adequately, as also stated earlier in this paper. 
Fortunately no similar irregularities have been 
found in the case of the \4686 line of Het, and 
from this line one can obtain a reasonable value 
of a, as will be shown presently. 


DaTA OF HousTON, AND CALCULATIONS 


As already noted, the value of R, given in 
G.C.1929 was based entirely on Houston’s ex- 
perimental value*® of Ru, with the use of “‘out- 
side" values of e/m and F in Eq. (16). Houston’s 
experimental value*® of Rye was not used at all, 
but it should have been used. It is therefore of 
interest to calculate the most reliable value of 
R., by the use of all of his data, and with the 
latest auxiliary constants, and to compare such a 
revised value with that obtained by DRW. In 
1929 I merely accepted Houston’s calculated 
value of Ru, but I will now recalculate it, just as I 
have done in the case of the DRW work. 

Houston calculated Ry from the measured 
wave-length, in air, of component 4 (also called 


27S, Pasternack, Phys. Rev. 54, 1113 (1938). 


I,) of Ha and of Hg, and he calculated Ry. from 
the measured wave-lengths of components I, and 
II, of the 44686 line of Het. The work was done 
with a Fabry-Perot interferometer, and all wave- 
lengths were measured against the 45015 helium 
line, whose wave-length in air was taken as 
5015.6750 I.A. This is also the present accepted 
wave-length, but the discrepancy that still exists 
—as will presently be shown—between his value 
of R, and that of DRW, must be ascribed in the 
main, according to DRW, to some error in this 
accepted wave-length. 

Houston, in his reduction to vacuum, naturally 
used the dispersion curve of Meggers and Peters. 
We shail now use Barrell and Sears" 1939 
dispersion curve and this change in the dispersion 
of air alone makes each calculated wave-length 
in vacuum about 0.004A larger, a by no means 
negligible change. As already noted, I will use 
+0.0036 cm~', in place of Houston’s +0.0056 
cm~', for the correction to the 1,4,5 complex of 
H, to get the wave number of component 1, but 
for Hg I retain his correction of +0.0024 cm™. 
The only additional quantity needed is the 
“Balmer” line correction, which for component 1 
of Hg is 3072) =0.0780 cm. 

Two further points should be noted. (1) Hous- 
ton measured the H, line on 5 plates, and the 
average wave-length in air of the 1,4,5 complex, 
as given correctly in his Table I, is 6562.8473 I.A. 
This wave-length, however, is written and used as 
6562.8475 in the calculations of his Table III, 
thus introducing a very small error. (2) The 
corresponding complex of Hs was also measured 
on 5 plates, and the average, as correctly stated 
in his Table I, and correctly used in his Table III, 
is 4861.3578 I.A. Because of a purely typo- 
graphical error this wave-length appears in his 
Table III as 4861.3678 and DRW, overlooking 
the true situation, point out (reference 6, p. 166) 
that his stated wave number of Hg in vacuum, in 
Table III, does not correspond to the stated 
wave-length in air, on the basis of the Meggers 
and Peters’ dispersion curve. 

With these changes I then obtain, by the same 
process as was used for DRW, 


Sy = 15,233.0113 from Hag, 
Su’ = 20,564.5643 cm~' from Hg, 


where Sq and Sy’ denote the ‘‘Balmer’’ lines. 
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Multiplication of sy by 7.2 and of sy’ by 16/3 
gives Ry, and the results are 


Ruy = 109,677.6813 cm from H, 
= 109,677.6762 from Hsg. 


The close agreement of these two values is a 
really remarkable verification of the entire theory 
of the hydrogen atom. The discrepancy is only 
+0.0051 cm (5 parts in 10%), as contrasted with 
a negative discrepancy of 0.021 cm resulting 
from Houston’s own calculated values of 
109,677.754 and 109,677.775. With his own cor- 
rection of +0.0056 cm in the case of Ha, and 
with Meggers and Peters’ dispersion curve, I 
obtain 109,677.755 from H, and 109,677.771 from 
Hg (discrepancy —0.016 cm~') as the values of 
Ru that Houston should have gotten from his 
data. If now the correction is changed to +0.0036 
cm, the discrepancy is increased in magnitude 
from —0.016 to —0.030 cm~. With the 1939 
dispersion curve, however, the corresponding 
discrepancies are +0.0195 cm~! and +0.0051 
cm~!. Thus the use of the +0.0036 cm~ correc- 
tion decreases the discrepancy between the values 
of Ry from Ha and Hg only when the 1939 
dispersion curve is used. With Meggers and 
Peters’ curve the discrepancy (now of opposite 
sign) is increased. Our final conclusion is, then, 
that in order to get a consistent value of Ru from 
Houston’s data on H, and Hsg, it is necessary to 
use the new 1939 dispersion curve and also the 
new +0.0036 cm~ correction. This fact may be 
considered as internal evidence in favor of the use 
of both. 

Houston gives relative weights of 3 to 1 to the 
two values of Ry from H, and Hg, and with this 
weighting I get the weighted average 


Ry = 109,677.6800 


which is 0.079 cm~ smaller than his published 
average, but 0.099 cm larger than the DRW 
result. If now we use the Meggers and Peters’ 
dispersion curve, but leave everything else 
unchanged, the weighted average rises to 


Ru = 109,677.7493 


which is 0.168 cm larger than the DRW result. 
The use of Houston’s +0.0056 cm~ correction 
raises this figure by only 0.0098 cm, and in- 
creases the 0.168 cm~' discrepancy by a like 


amount to 0.178 cm. Thus the use of the newer 
values for the dispersion curve of air decreases 
the discrepancy between Houston’s and DRW’s 
values of Ru only from 0.168 to 0.099 cm-', and 
the other changes that have been made in 
Houston’s calculations are considerably smaller. 

The final unexplained discrepancy between the 
results of DRW and Houston may be stated 
more clearly as follows. Houston measured 
6562.8473 I.A. as the wave-length in air of the 
1,4 complex of Hz. DRW’s published wave 
number in vacuum, when reduced to wave-length 
in air by the 1939 dispersion curve, is 6562.8533 
I.A.—just 0.006A larger. This difference is many 
times the apparent probable error of the measure- 
ment but, as already noted, Houston’s result is 
based on the 5015 helium line, whereas the 
DRW result is based on the \6438 cadmium line. 

This completes the discussion of Houston's 
value of Ry, and we now consider his value of 
Rue, which is based on the measured wave- 
lengths of components I, and II, of the \4686 
line. For this line R= Rue, Z=2, n’=4 and n” =3, 
in Eqs. (29) and (30). For the I, component, 
j' =7/2, =5/2. For the II, component j’=5/2, 
j'' =3/2. These are the two strongest compo- 
nents, out of a total of 8 predicted components, 
and are of nearly equal intensity. Chu®* gives a 
diagram of the fine-structure pattern. Because of 
the great intensity of the two components used, 
as compared to the remaining components, 
Houston* used merely the observed wave-lengths, 
without any additional correction (but see foot- 
note 30). Chu*®* later made a more detailed 
analysis of the microphotometric curves. His 
results will be discussed presently. 

Houston obtained component I, on 7 different 
plates, with an average wave-length in air of 
4685.8030 I.A., and a mean deviation of 0.0026 
I.A. The II, component was measured also on 7 
plates, with an average of 4685.7030 I.A., and 
mean deviation 0.0012 I.A. In this case we can 
use the difference of these two wave-lengths to 
obtain an experimental value of the fine structure 
constant a. With the usual dispersion curve," one 
obtains for the wave number difference, in 
vacuum, As=0.45533 cm. This should equal 
(101/1296)Ruea, thus giving 1/a=137.038, 


Y, Chu, Phys. Rev. 55, 175 (1939). 
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which all recent experimental evidence indicates 
is close to the true value. With this value of a, the 
negative correction of the I, component to give 
the “Balmer” line is 
=0.19723 cm~!. The negative correction for 
the Il, component is 16Ry-a?(1/108 — 7/3072) 
= 0.65256 cm~. When one thus uses the value of 
a derived from the data, the two resulting values 
of the ‘Balmer’’ line are necessarily identical, 
and the result is 


Ste = 21,324.89082 


Multiplying by 36/7 (the reciprocal of 4(1/9 
—1/16) in Eq. (30)), we get Ru. =109,722.295, 
cm, which is 0.0322 cm~ larger than the value 
already calculated indirectly from DRW’s data. 
The agreement is, accidentally, much better than 
that between the two directly calculated values of 
Ru, by Houston and by DRW. Houston himself, 
using the Meggers and Peters’ dispersion curve 
and an “outside” value of a, got from the two 
components slightly different values of Ru. and a 
weighted average value’ of 109,722.403+0.004 

From Houston's values of Ry = 109,677.6800 
and Rye=109,722.2957, as just recalculated, we 
can now obtain E and e/m by means of Eqs. (14) 
and (15). The results are E=5.47829X10-‘ and 
e/m=1.7617,X10" e.m.u. No probable error is 
calculated for this value of e/m, since it is not 
used in the final average, as discussed later. If 
Houston’s H, correction of 0.0056 cm were 
used, we would get 1.76213. In G.C.1929 (p. 47) I 
got, from the same data, e/m=1.7608+0.0008. 
In 1938 I obtained” 1.7601;+0.0008. These 
changes indicate very well the effect of varying 
values of auxiliary constants, and varying 
methods of reduction of the same data. 

We can also obtain R,, by means of Eq. (21). 
The result is R,=109,737.3125 cm~', which is 
only 0.010 cm~ above the adopted DRW result. 
Houston* himself gives R,=109,737.424 cm™ 
and this value was adopted in G.C.1929 (p. 49). 

29 Each of Houston’s stated errors appears to be based on 
external consistency, but since only two values are available 
in each case, internal consistency would yield a far more 
reliable estimate. I have not carried through a calculation 
of probable error such as I did in the case of the DRW 
results, because I am adopting the DRW results. The sole 
purpose of the present sation is to get a correct com- 
parison with the DRW results, by the use of the same 


auxiliary constants and the same method of calculation in 
each case. 


In the case of DRW we have calculated R, also 
from Ry and from Rp separately, using Eqs. (1’) 
and (2’) and outside values of e/m and F to get E 
from Eq. (4). The result from Rp is the more 
accurate. In a similar way, we can get R, from 
Houston’s value of Ry and from his value of Rue, 
by the use of Eqs. (1’) and (3’). These two new 
values are 109,737.4024 cm and 109,737.3351 
cm~', respectively. Both of these values are 
higher than that obtained from Eq. (21), and 
differ considerably from the DRW result. This 
shows that the previous close coincidence in the 
values of R, from the DRW and the Houston 
work is purely accidental, and is due to the fact 
that real discrepancies in the values of Ry and 
e/m (or E) happen nearly to cancel each other, in 
the calculation of R,. 


Data OF CHU, AND CALCULATIONS 


Chu,”* as already noted, has analyzed the fine 
structure of He \4686 and has calculated a value 
of Rue. His work was done under the direction of 
Houston, with similar apparatus to that used by 
Houston. The wave-length measurements of both 
men are based on the 45015 helium line. Hence it 
seems legitimate to combine Chu’s value of Rue 
with Houston’s value of Ru, to get a value of R,, 
and also of e/m. 

In Houston’s work the two strong maxima of 
the intensity curve of \4686 were identified with 
components I, and II. In Chu’s work a careful 
analysis of the intensity curve was carried out, in 
much the same way as was done by Williams and 
Gibbs'® for H, and D,. In the case of 44686, the 
fine structure extends over so large an interval, 
2.43 cm~', that with the usual orders of interfer- 
ence of the Fabry-Perot interferometer, there is a 
great deal of overlapping of orders. By varying 
the order of interference, the resulting intensity 
pattern can be made to vary markedly. The 
corrections necessary to apply to the two strong 
maxima were calculated by Chu, from his analysis 
of the intensity curves, for each order of inter- 
ference used, and these corrections, as might be 
expected, also varied greatly. In fact, for the IT, 
component the corrections varied from +0.0012A 
to —0.0029A. This variation is really an advan- 
tage, for if the final corrected wave-length is 
constant from plate to plate, it is very good 
evidence that the procedure employed is funda- 
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mentally sound. This result was actually attained. 
Thus the resulting wave-lengths*®® of the II, 
component, on the four plates used, varied only 
by 0.0009A. It was more difficult to measure the 
I, component and on one plate (No. 126) its 
wave-length differs by 0.0037A from the average 
of the remaining three plates. Chu states that the 
position of the maximum on the microphotometric 
tracing of this plate was particularly uncertain. 

We shall first follow Chu in accepting the 
results of all four plates, and we shall also use his 
assigned weights, which run from 4 to 1. Then 
for the I, component the average wave-length in 
air is 4685.8012A, giving s=21,335.09626 cm 
for the wave number in vacuum (with the 1939 
dispersion curve of air). For the Il, component 
the corresponding results are 4685.7017A and 
21,335.54931 cm~'. These wave-lengths are re- 
spectively 0.0018A and 0.0013A larger than 
Houston’s values. The wave number difference 
is, from Chu’s data, 0.45305 cm", from which we 
get 1/a=137.383. As we have seen, Houston’s 
corresponding result is 1/a= 137.038. Chu does 
not give explicitly his value of 1/a@ (137.383), but 
I had calculated it and until the recent work of 
Christy and Keller* I had considered it the best 
experimental determination, with an assigned 
probable error of +0.4. It is to this value by Chu 
that I had reference in my own recent letter® on 
the subject. 

With Chu’s value of a, the “Balmer’’ correc- 
tions for the I, and II, components are now 
0.19625 and 0.64930 cm", respectively. The 
necessarily identical value of the Balmer line, 
from either component, is SHe=21,334.90001 
cm~', and from this we get 


Rue = 109,722.3429 cm™, 


which is 0.047 cm~ larger than Houston’s value, 
and 0.079 cm larger than DRW’s indirect value. 
Chu himself gives 109,722.43+0.04 cm, using 
Meggers and Peters’ dispersion curve. 

We can now combine Houston’s value of 


*® It is only fair to remark that Houston’ also recognized 
clearly the need for making corrections on each individual 
plate, to get the true positions of the I, and II, com- 

nents. He tried, however, to vary the order of inter- 

erence in such a way that the various resulting corrections 
would, on the average, be expected to cancel out. 
a ono, F. Christy and J. M. Keller, Phys. Rev. 58, 658 

® R. T. Birge, Phys. Rev. 58, 658 (1940). 


Ry =109,677.6800 and Chu’s value of 
Rue = 109,722.3429 to get, from Eqs. (14), 
(15), and (21) 


E=5.484086 x 10-4, 
e/m=1.7598,X 107 e.m.u., 
109,737.3756 


It has just been noted that component I, was 
observed to be very ill-defined on one of Chu’'s 
plates (No. 126), and the results given by this 
plate are in definite disagreement with those of 
the remaining three plates. The disagreement is 
shown clearly by calculating 1/a from each plate 
separately. The results are 


Plate 1/a 
118 137.314 
124 136.971 
126 139.790 
134 136.971 


The weighted average value, omitting No. 126, is 
137.079 and if I had noticed these facts earlier, I 
would certainly not have taken Chu’s weighted 
average of all four plates (137.38) as the best 
experimental value of 1/a. This seems to be a 
situation where it is clearly justifiable to discard 
one plate, and hence it would appear that Chu’s 
best value of 1/a is more nearly 137.1+0.2, as 
compared to the value 136.95+0.05 obtained 
from x-ray data.*" Although Houston’s value 
(137.04) is also very satisfactory, it is presumably 
not as reliable as Chu’s, because of Houston's 
less detailed analysis of his intensity curves. 

With plate No. 126 omitted, and with 
1/a=137.079, one gets, from Chu’s remaining 
three plates 


Rue = 109,722.3280 


which is only 0.0156 cm larger than Houston's 
value. Combining this last result with Houston’s 
Ru, we obtain 


E=5.48225X10-, 
e/m=1.7604s X10? e.m.u., 
R, = 109,737.3556 cm™. 


These results are, I believe, the most reliable that 
can be obtained from the combined work of 
Houston and Chu. Their work and the work of 
DRW furnish all of the experimental data from 
which one can obtain a precision value of R,. In 
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adopting a final value of Rz one should really 
consider some weighted average of the two 
results. But because of the character of the 
discrepancies, as already discussed, such an 
average appears to have little meaning. For this 
reason, and for this reason alone, I have adopted 
as the final value of R, merely the DRW 
result.** 

On the other hand, the value of e/m just given 
is to be combined with many other values, later 
in this paper, to get a final best value. Hence it 
is necessary to adopt a probable error for it. This 
value of e/m has been calculated from Eq. (15), 
and the only probable errors that need be con- 
sidered are those in F and in the difference 
Rue— Ru. The adopted value of F is 9651.49+1.0 
abs. e.m.u., as already stated. 

Houston's two values of Ry (from H, and Hg) 
differ by only 0.0051 cm~', and hence 0.010 cm~ 
seems a liberal estimate for the probable error of 
the weighted average (109,677.6800 cm~). Chu’s 
value of Ry., from 3 plates, differs from that 
calculated from all 4 plates by 0.015 cm~. Since 
the omitted plate is obviously poor, 0.010 cm~ 
again seems a liberal estimate for the probable 
error of the 3-plate result (109,722.3280 cm~'). 
Now both Houston and Chu base all their wave- 
lengths on He 5015. Hence the most serious 
sources of possible systematic error do not appear 
in the difference Rue—Ru. We therefore get 
Rue— Ru = 44.648)+0.014 and finally 


e/m = (1.7604s+0.0005s) X 107 abs. e.m.u. 


as the adopted result from the Houston-Chu 
work, 


33 Professor Houston has called my attention to the fact 
that the relations between Ry, Rp, Rue and R.., as obtained 
by various observers, can conveniently be shown on a 
graph of the same general character as the Birge-Bond 
diagram. Thus in Eqs. (1’)—(3’), if we use an assumed value 
of E(=Epo), each experimental value of R, will yield an 
apparent value of R..(= R..’). If now these values of R..’ are 
plotted against the reciprocal of the nuclear atomic weight 
(i.e., against 1/H*, 1/D* etc.) as abscissa, all discrepancies 
are immediately obvious. If there are no discrepancies, all 
of the points will lie on one straight line. The extrapolation 
to zero abscissa of the straight line through any two points 
gives the value of R.. resulting from the use of these two 
data. On such a diagram, the DRW data give one straight 
line, and Houston's data give another line with a quite 
different slope, corresponding to the quite different value of 
e/m yielded by his data. The Houston-Chu data give a line 
of intermediate slope. The diagram shows immediately that 
no compromise line will give a satisfactory representation 
of the data of any one of the different observers. 


EXPERIMENTAL VALUES OF ¢/m FROM FINE 
STRUCTURE ANALYSIS 


a. Work of R. C. Williams 


Although we have now covered the experi- 
mental data leading to values of R., and inci- 
dentally of e/m, there are further similar 
investigations from which one can obtain experi- 
mental values of e/m but not of R,. I have not 
examined this additional work with the thor- 
oughness of the preceding, and the following 
discussion is rather brief. 

One important additional investigation on the 
fine structure of H, and D, is that of R. C. 
Williams, which was carried out under the 
direction of R. C. Gibbs. In addition to the 
papers giving the final results on e/m (reference 
1) and on a@,* and the papers already quoted,'* '* 
there are other earlier reports.**-** 

In order to calculate e/m, the difference of 
wave-length of component 1 of H, and D, was 
measured with an interferometer of over 250,000 
resolving power (3-mm etalon silvered to reflect 
over 90 percent). The method of analyzing the 
fine structure is discussed in references 34 and 19. 
Williams’ published average ‘“‘peak-to-peak”’ 
interval in air for the H, and D, lines is 
4.14716+0.0004 cm~'!. Division by the index of 
refraction of air (1.00027624) gives 4.146014 cm~! 
for the interval in vacuum. The use of Meggers 
and Peters’ 1 = 1.0002757 changes this result only 
to 4.146017, but Williams gives 4.14610 cm~ as 
his value in vacuum. 

Full details of the corrections to be applied to 
each complex “‘peak”’ to get the position of 
component 1 are given by Williams in his Table 
III of reference 1. He finds a differential correction 
of +0.0009 cm~', but most of this correction is 
due to a differential correction of +0.0012 cm 
for the effect of component 3. As noted earlier, 
DRW make no differential correction to their 
H,—D, interval, but their printed results include 
the correction for component 3. Hence there is no 
real disagreement between Williams and DRW. 

Williams’ corrected interval is thus 4.14691, 

* R. C. Williams, Phys. Rev. 54, 558 (1938). 

36 R. C,. Gibbs and R. C. Williams, Phys. Rev. 44, 1029 


(1933). 
3° R. C. Williams and R. C. Gibbs, Phys. Rev. 48, 971 


(1935). 
37 R. C. Williams and R. C. Gibbs, Phys. Rev. 49, 416 
(1936), (two abstracts). 
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cm and Eq. (11) then gives E. For convenience 
in comparing the results of different observers, 
Eq. (11) may be written as 


E= (1.3238985 X (sp — sn) (11’) 


with the DRW value of sq (15,232.9974 cm-) 
and the atomic weights already listed. Therefore 
E=5.49009;X10~ from Williams’ data. Then 
from Eq. (4), e/m=1.7579;X 107 e.m.u. Williams 
gives (5.4902+0.0005)X10-* and (1.75794 
0.0004) x 107. 

In the case of DRW we have already found the 
result e/m= (1.75913+0.00037) X10’, where the 
probable error was calculated with some care, and 
includes a certain allowance for systematic error. 
In order to get a final weighted average value of 
e/m, it is necessary to assign relative weights to 
the various results. Hence it is only the relative 
probable errors that are important, and in the 
case of all the remaining results, we ask ourselves, 
“If 0.00037 X 10’ is taken as a fair estimate of the 
probable error of the DRW value of e/m, what 
is a proper estimate for the new investigation, 
when one considers the relative accuracy of the 
instruments used, the detailed measurements 
made, etc.?”” The probable errors adopted here 
are therefore frankly purely personal estimates. 

For the R. C. Williams result I adopt 


e/m = (1.7579;+0.0005) X 107. 
b. Work of C. F. Robinson 


Robinson has as yet published only an ab- 
stract*® of his investigation, which is a continu- 
ation of work started by Houston*® on an analysis 
of the fine structure of H, and D, by a new 
method. I am, however, indebted to Professor 
Houston for further details. Robinson’s measured 
separation of component 1 of H, and Da, in air, is 
As = (4.1453+0.0010) cm~'. This becomes 4.14415 
+0.0010 cm in vacuum, as contrasted with 
Williams’ 4.14691 and with Houston's 
preliminary values*® of 4.1440 and 4.1399 cm}. 
Accepting Robinson's estimate of his own 
probable error in As, and using the probable 
error in F already given, we obtain 


E=(5.4864,+0.00132) X 
e/m= (1.75914+0.00046) X 10’. 


and 


38 C, F. Robinson, Phys. Rev. 55, 423A (1939). 


8° W. V. Houston, Phys. Rev. 51, 446 (1937). 


The resulting probable error in e/m is essentially 
that adopted for Williams’ result, and the data 
just given thus show directly the error one needs 
to assume in the H,—D, interval to get my 
adopted error in e/m. I assume the same accuracy 
for the work of Williams and of Robinson and 
adopt for Robinson’s result 


e/m=(1.75914+0.0005) X 10’. 


c. Work of Shane and Spedding 


The last investigation on the fine structure of 
H, and D, to be discussed, but the earliest in 
point of time, is that of Shane and Spedding.” 
They used an interferometer with a plate sepa- 
ration of about 7.8 mm. They make an elaborate 
analysis of their intensity curves, and give their 
final separations of H, and D, components in 
terms of orders of interference, but not in terms 
of wave-length or wave number. Each result is 
accompanied by its corresponding value of e/m. 
Under the circumstances the simplest procedure 
appears to be merely to accept their final average 
value of e/m, which is (1.7579+0.0003) x10’, 
and then correct it for the change in values of 
auxiliary constants. They used Eq. (12) and 
in this equation their adopted values give 
F(D—H)/H+*D =4786.383, whereas my present 
adopted values give 4785.7526. The value of e/m 
is then to be reduced in this ratio. But according 
to Williams and Gibbs,** Shane and Spedding 
did not reduce their wave number interval from 
air to vacuum, although in Eq. (12) they used the 
wave number in vacuum of Ry. To correct for 
this error, one should multiply their result by the 
index of refraction of air (1.000276), thus giving 
1.0001444 as the resulting factor for the two 
almost compensating corrections. The revised 
result, with my own adopted probable error, is 


e/m= (1.7581;+0.0006) X 107 e.m.u. 


EXPERIMENTAL VALUES OF e/m BY 
OTHER METHODS 


In 1938 I obtained,” as the best average value 
of e/m from all methods, (1.7591+0.0003) x 10’ 
abs. e.m.u. The first four determinations listed 
then comprise the work on H, D, and Het, and 
revised values for all four have now been ob- 


“C. D. Shane and F. H. Spedding, Phys. Rev. 47, 33 


(1935). 
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tained, together with the new value by DRW.® 
The remaining six determinations given in the 
1938 letter, together with one new determination, 
will now be considered. No attempt is made at 
critical discussion, and in fact my only object 
now is to explain, as briefly as possible, why my 
present adopted value, or its assigned probable 
error, deviates from that given by the author. In 
general it is difficult, in the case of the methods 
now to be considered, for an outsider to make any 
intelligent evaluation of the probable error, and 
hence I have, in general, merely accepted that 
published by the author. It is quite possible that, 
relative to the adopted probable errors in the 
preceding sections of the paper, all of the 
probable errors to follow should be slightly 
increased. 

When e/m is obtained from fine structure 
analysis, the value of Fig is involved, where F; is 
the value of the Faraday in international electric 
units, and F(= Fg) its value in absolute units. 
Now as already discussed in footnote 23, I have, 
subsequent to the 1938 on raised the 
value of F; and lowered that of g by almost equal 
amounts. Hence the resulting values of e/m are 
virtually unchanged, although additional small 
changes have now been made in certain instances. 
On the other hand, in the case of five of the six 
remaining methods that have been used in a 
precise determination of e/m, F; does not appear 
at all in the equation, but g appears in the 
denominator. Hence in all these cases e/m is 
raised by about 0.00012 X10’ e.m.u., as a result 
of the lowering of gq from 0.99993 to 0.99986. In 
the sixth method (refraction of x-rays), Fig 
occurs in the denominator and again the value of 
e/m is virtually unchanged. 


a. Zeeman effect 


A number of investigations on the value of e/m 
from Zeeman effect have been carried out by 
Houston and co-workers. The results of all of their 
work are tabulated by Kinsler and Houston." 
These comprise the final results of work on Zn, 
Cd, He, and Ne. They give as a final weighted 
average, e/m=(1.7570+0.0007) X10’. It is not 
apparent how the probable error was obtained. 
Accepting their own probable errors for the four 


“tL. E. Kinsler and’ W. V. Houston, Phys. Rev. 46, 533 
(1934). 
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individual results, I obtain as the weighted 
average 1.75694.+0.00017 by external con- 
sistency, and +0.00047 by internal consistency. 
Since, however, it is well to make some allowance 
for systematic errors, I will increase this latter 
error to +0.0007, to agree with their own 
assigned error. The values of the auxiliary con- 
stants are not stated, but in an earlier paper® the 
values g=0.99995 and c=299,796 km/sec. were 
used, and I believe were used also in all suc- 
ceeding papers of the series. In 1934 I recom- 
mended** c= 299,776+4 km/sec., and the inclu- 
sion of more recent results leaves unchanged both 
this average and its probable error. In the 
Zeeman effect work the change in ¢ decreases e/m 
by almost the same amount that the change in q 
increases it. The exact corrected result is 1.7569s;. 
Hence I adopt 


e m= (1.7570+0.0007) X 107 e.m.u. 
b. Refraction of x-rays 


Bearden," in an elaborate investigation, ob- 
tained e/m=(1.7601+0.0003) 107. All of his 
auxiliary constants agree with those I was using 
at the time. Hence the only essential change now 
to be made is that resulting from the very small 
changes in F and C. His stated probable error is 
very small, but is the result of a detailed con- 
sideration of the errors of all factors in the final 
expressions for e/m. His exact probable error in 
e/m, based on F=9651.3+0.80, is +0.0003;:. 
With the new F=9651.49+1.0 and C=12.01465 
+0.00023, his result becomes 1.7600¢+0.00033. 
In order to make at least a slight allowance for 
systematic error I adopt 


e/m= (1.7600¢+0.0004) X 107 e.m.u. 


c. Direct velocity measurements 


Two careful investigations by this method 
have been made, the first by Perry and Chaffee,* 
and the second by Kirchner.** Perry and Chaffee 
obtained six series of results, with an actual 
average of 1.7609+0.0002, where the error is 


“ 3; S. Campbell and W. V. Houston, Phys. Rev. 39, 601 
(1932). 

*R. T. Birge, Nature 134, 771 (1934). 

“J. A. Bearden, Phys. Rev. 54, 698 (1938). 
( 930) T. Perry and E. L. Chaffee, Phys. Rev. 36, 904 
1 

46 * Kirchner, Ann. d. Physik 8, 975 (1931) and 12, 503 
(1932). 
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I. 
REFERENCE 

e/m P. E. IN TEXT METHOD S(Fipq) 
1.75913 | 3.7 6 Haq and Da Fiq 
1.75797 | 5.01074 1 Hg and Da iq 
1.75914 | 5.01074 38 Haq and Da Fiq 
1.75815 | 6.01074 40 Hg and Dg Fiq 
1.76048 | 5.8 xX1074 3,28 Ry and Rye Fiq 
1.75700 | 7.01074 41 Zeeman effect 1 
1.76006 | 4.0 X10~4 44 X-ray refraction 1/Fiq 
1.76110 |10.0 X10~4 45 Direct velocity 1/pq 
1.75900 | 9.0 X10~¢4 46 Direct velocity 1/pq 
1.75982 | 4.01074 48 Magnetic deflection 1/q 
1.75820 -|13.0 X1074 49 Crossed fields p/aq 
1.75870 | 8.0 $1 Busch method p/a@ 


based on external consistency. They however 
adopt a final probable error of +0.0010, and I 
shall do the same. The value of c enters into their 
formula in a rather complex way and I find, by 
actual trial, that the change of c from 299,796 to 
299,776 produces a change in e/m of only about 
one part in 10°, which is negligible. They used 
pq= 1.00046 to reduce from int. to abs. volts, 
and with the new value 1.00034, one gets 
e m=1.7611,. I therefore adopt 


e m= (1.7611+0.0010) X10? e.m.u. 


Kirchner** obtained, in two separate investi- 
gations, 1.7585+0.0012 (revised from the original 
1.7598+0.0025) and 1.7590+0.0015. From these 
results I obtained‘? a weighted average of 
1.7587+0.0009, where the error is based on 
internal consistency. This average has since been 
quoted in the literature as Kirchner’s final result. 
Kirchner gives corresponding figures in int. and 
abs. volts that imply an assumed pg=1.0005. 
With our new value, the weighted average 
becomes 1.7589, and I adopt 


e/m = (1.7590+0.0009) X 107 e.m.u. 
d. Magnetic deflection 


Using a new magnetic deflection method 
devised by E. O. Lawrence, Dunnington** has 
made a most extensive and refined investigation 
of the value of e/m. His final published result 
is 1.7597+0.0004. This result is based on 
q=0.99993. With the new value of g one obtains 


e/m=(1.7598.+0.0004) X 107 e.m.u., 
which I adopt. 


‘7 R. T. Birge, Phys. Rev. 42, 736 (1932). 
‘8 F. G. Dunnington, Phys. Rev. 43, 404 (1933) and 52, 
475 (1937). 


e. Crossed fields 


Shaw** has used a method of crossed electric 
and magnetic fields to obtain a published value of 
e/m=1.7571+40.0013. He, however, failed . to 
make the reduction from international to abso- 
lute electric units, and since the factor V/H? 
appears in his equation ( V = potential difference, 
H=magnetic field), the proper reduction factor 
is p/q= 1.00062 for the values of p and 
now chosen. Hence his corrected result is 


e/m=(1.7582+0.0013) X10? e.m.u., 
which I adopt. 


f. Longitudinal magnetic field (Busch method) 


In 1929 the best work on e/m with the Busch 
method was that by Wolf.®° His investigation is 
discussed on pp. 43-44 of G.C.1929. At that time 
I assumed that he had not made the correction 
from international to absolute units, an assump- 
tion that was later verified. As in the case of 
crossed fields, the factor V/H* again appears in 
the equation. In 1929 I corrected Wolf's result 
from 1.7679 to 1.7689. With the latest p/g value, 
it becomes 1.7690. 

A recent and much more accurate determi- 
nation by this method has been made by 
Goedicke.*! He does not give the values of p and g 
used, but from certain published figures, his 
value of g is evidently about 0.99991. I assume 
that his value of p is 1.00050, since that value is 
commonly used. On these assumptions, his value 
of e/m=1.7586 should be raised only to 1.7587. 
He gives 0.0023 as a maximum possible error, and 
following my usual custom I will take 1/3 of this 
as the probable error. Hence I adopt 


e/m= (1.7587 +0.0008) X 107 e.m.u. 


for his final result. 


SUMMARY OF VALUES OF e/m 


We have now discussed and adopted values of 
e/m from twelve different investigations, covering 
seven essentially different methods. All of the 
adopted results and probable errors are listed in 
Table I. So far as values of the auxiliary constants 


49 A. E. Shaw, Phys. Rev. 54, 193 (1938). 
%® F. Wolf, Ann. d. Physik [4 83, 849 (1927). 
5t Egon Goedicke, Ann. d. hysik [5] 36, 47 (1939). 


1 
4 
1 
q 
s 
g 
V 
ll 
is 
l- 
il 
n 
5 
3. - 
or 
d 
al 
is 
04 


784 


are concerned, the chief sources of recent change 
and controversy are the values of F; (the Faraday 
in international units), p and qg. The last column 
shows the explicit dependence of e/m on these 
three constants. 

The weighted average of the 12 results is 
e/m=1.75920.. The probable error by external 
consistency (r,) is 0.000192, and by internal con- 
sistency (r;) is +0.000153. The ratio r./r; is 1.213, 
which is very satisfactory. In 1938 I got” for the 
corresponding results 1.75909+0.00024 (r.) and 
r./r;=1.51, and recommended 1.7591+0.0003 as 
the best value to adopt. Although the new 
average (1.759206) is only slightly larger than the 
value 1.7591, which has been rather widely used 
for the past three years, it seems best to adopt it. 

It is, however, desirable to consider in some 
detail the internal consistency of the data since, 
as noted in my 1938 paper,” this is a subject that 
has been under discussion for many years. In 
that paper I concluded that the annoying previ- 
ous discrepancy between the ‘‘spectroscopic”’ and 
the ‘‘deflection’”’ values of e/m had virtually 
disappeared. Bearden®™ definitely disagreed with 
this conclusion, and even disagreed with my 
classification of his x-ray refraction value of e/m 


as a “spectroscopic” result. So far as this last. 


point is concerned, I think it is mainly a question 
of how the two groups of results are defined. 
Although in 1929? I gave two values of e/m, 
commonly referred to as the “‘spectroscopic’’ and 
the ‘‘deflection”’ values, I never for one moment 
believed that there are actually two different 
values. I assumed then, and I still assume, that 
any apparent discrepancy is due to errors of 
theory or of experiment, or both. Certain values 
of e/m are obtained with spectroscopic equipment 
and involve measurements of spectral wave- 
lengths. These I have grouped together to obtain 
a weighted average ‘‘spectroscopic’”’ result, and 
Bearden’s result certainly belongs in this class. 
The other experiments are electric and magnetic 
in character, and involve linear acceleration and 
deflection of free electrons over comparatively 
long paths. The technique and quantities meas- 
ured are totally different from those in the 
spectroscopic work. Hence the systematic errors 
are also likely to be totally different, and to such 


8 J. A. Bearden, Phys. Rev. 55, 584 (1939). 
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errors I attribute any outstanding difference in 
the values of e/m that may exist. 

Bearden® gives a rather close argument to 
show that his value of e/m is concerned with 
‘free’ electrons, and hence should be classed 
with the deflection experiments. With this as- 
sumption he gets two values of e/m (1.7583, 
+0.0002s and 1.7598;+0.0001;), and he states 
that there is less than one chance in 10° of such a 
discrepancy being due to chance. His data, 
however, give for the discrepancy 0.00157 
+0.0003276, or 4.79 times its probable error, and 
there is one chance in 826 of this occurring. If 
these last odds are taken at face value, a really 
serious discrepancy exists, but there seems to be 
no point in exaggerating the situation. 

In considering the question now, I will omit 
Bearden’s x-ray result from both groups. Then the 
first six results in the table comprise the spectro- 
scopic group, and these six give 


e/m=1.75880+0.0002829 (r-) 
ratio 1.355, 


or +0.00020s; 


whereas the last five results of the table comprise 
the deflection group, and give 


e/m=1.75959+0.0002407 (re) 
ratio 0.785. 


or +0.00030¢; 


With the probable errors based on external 
consistency (r,), there is one chance in 6} that 
this discrepancy is accidental (the discrepancy 
being 2.13 times its probable error), and an 
almost identical result is obtained from the r; 
values. Although one would naturally wish that 
the odds were more nearly even, a discrepancy of 
this magnitude is not considered by statisticians 
as of any real significance. 

If now we compare Bearden’s result (point 7 of 
the table), with the two averages just found, and 
use r, errors, we find that the odds are just even 
that his result is consistent with the deflection 
average, but there is only one chance in 12 that 
his result is consistent with the spectroscopic 
average. One possible conclusion from these facts 
is that point 7 is definitely too high. It is, 
however, worth while to explore the matter 
further. 
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It has already been mentioned in footnote 23 
that the value of F now adopted is almost 
identical with that advocated in 1938, but the 
coincidence is only accidental, and is actually the 
result of two almost compensating changes that 
have been made. 

Since F appears in different equations for e/m 
in different positions, it is possible to evaluate it 
merely from the e/m data. The resulting value 
has, however, a larger probable error than the 
value obtained from electrolysis. Thus the 
weighted average of points 1-5 inclusive is 
e/m=1.75897 +0.0002645 (r-) or +0.000217 
The ratio 7./r; is 1.214, which is satisfactory. If 
now we calculate the value of F that will bring 
this result into agreement with Bearden’s point 7, 
the result is 


F=9651.40(1.76006/1.75897)! 


= 9654.39 abs. e.m.u., . 


and the new value of e/m is [(1.76006) (1.75897) |}! 

= 1.75955. This is very close to the ‘‘deflection”’ 
average of 1.75959, which does not involve F at 
all. In fact, the value of F that will force the fine 
structure average (1.75897) to agree with the 
deflection average is 


F= (9651.40) (1.75959) /(1.75897) 
= 9654.80 abs. e.m.u. 


The two new values of F just computed corre- 
spond to 96,531 and 96,535 int. coulombs 
(chemical scale), values that are far above the 
direct results given by either the silver voltameter 
(96,494) or the iodine voltameter (96,511). Such 
high values of F are clearly ruled out, on the 
basis of all experimental work with voltameters, 
but it is certainly worth noting that the value 
F=96,533 int. coulombs (chemical scale) will 
bring into almost perfect agreement all values 
of e/m, except that from Zeeman effect (point 6), 
which is certainly too low. 

My final conclusion is then that, on the basis 
of all experimental work to date, 1.7592 is the 
best value of e/m. But if it becomes possible, as a 


result of the future discovery of new sources of 
systematic error in F, to raise its value some 30 
coulombs above the present adopted result of 
96,501.3 int. coulombs (chemical scale), then the 
consistency of the present e/m results will be 
greatly improved, and the best value of e/m will 
become 1.7595 or 1.7596. The weighted average 
of all 12 results, as already given, is 1.7592; 
+0.0002, but just because of this possible 
uncertainty in regard to the value of F, I shall 
enlarge the probable error to 0.0005. Hence the 
final recommended value is 


e/m=(1.7592+0.0005) 107 abs. e.m.u. 
or 
(e/m)c = (5.2736.+0.0015) X10" abs. e.s.u. 


Then, from the adopted values of e/m and F 


E=(5.4862,+0.0017) 
(physical scale) (derived). 


SUMMARY 


As a matter of convenience, the final values 
adopted in this paper are now listed. A list of 
adopted values of auxiliary constants has already 
been given, as well as the individual values of 
e/m. All of the following results are based 
entirely on the work of DRW,* except the value 
of e/m and the derived value of E. 


Ru = 109,677.5812+0.007; (I.A. scale), 
Rp = 109,707.419;+0.007; (I.A. scale), 
Rue = 109,722.263 +0.012 cm (I.A. scale), 
R,, = 109,737.303 +0.017 (I.A. scale), 
or +0.05 (c.g.s. system), 
e/m= (1.7592 +0.0005) X10’ abs. e.m.u., 
(e/m)c = (5.2736¢+0.0015) X10" abs. e.s.u., 
E=(5.4864;+0.0009,) 
X10- (physical scale) (direct), 
E=(5.4862,+0.0017) 
X10- (physical scale) (derived). 


In general calculations, the derived value of E 
should be used, in order that it may be consistent 
with the adopted values of F and e/m. 
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Further Resolution of Two Parallel Bands of Ammonia and the Interaction between 
Vibration and Rotation 


Hsi-yin SHENG, E. F. BARKER AND D. M. DENNISON 
Harrison M. Randall Laboratory of Physics, University of Michigan, Ann Arbor, Michigan 
(Received October 9, 1941) 


The fundamental low frequency parallel band at 10u and 
its upper stage band at 16u have been re-examined under 
high dispersion. It is found that most of the fine structure 
lines hitherto reported as single are resolved into multiplets. 
The positions of the lines may be accounted for accurately 
by assuming that the energy levels have the form, 


Wy" */he+A + Cy"*K?+ F(R). 


The superscripts s or a@ indicate whether the level in 
question is the lower or the upper of the two levels produced 
by the double minimum potential. F(JK) is the small 
centrifugal distortion term calculated by Slawsky and 
Dennison. The constants Ao’, Ao*, Ai*, Cot 
and C,*—Co* are determined from the experimental 
data. The following vibrational energies are also found, 
=932.24hc, Wi* =968.08hc and = 1597.4hc. 

In the discussion of the rotational constants, it is pointed 
out that two quantities are involved, (a) the average 


constants, averaged over the two levels of the double 
minimum, as }(A,*+A,*), and (b) the difference of the 
constants, as A,*—A,*. Recently Shaffer has calculated the 
theoretical expressions for the averaged constants and has 
shown how these depend upon the anharmonic terms in the 
potential energy. A treatment of the difference of the 
constants is equivalent to a treatment of the change in the 
splitting of the two levels of the double minimum due to the 
vibration-rotation interaction. This may be computed 
from the expression for the splitting of the double minimum 
levels given by the W-K-B method of approximation. The 
calculated result for the change in splitting of the first 
vibration levels is 


(64/hc) theo = —0.162(J?+ J) +0.222K?, 
which is in excellent agreement with the observed change 
(64/he)exp = 


The significance of the agreement is discussed. 


INTRODUCTION 


HE general problem of the evaluation of the 

vibration-rotation energy levels of poly- 
atomic molecules has received considerable atten- 
tion in recent years and much progress has been 
made. The present paper contributes to this 
problem and is concerned with the high resolu- 
tion of the low frequency parallel bands of 
ammonia. It will be shown that certain novel 
features appear here which are produced by the 
double minimum potential of ammonia and 
which are not likely to occur in measurable 
degree in other spectra. 

The Hamiltonian function representing a poly- 
atomic molecule is separable only in zeroth 
approximation and the usual procedure is to 
develop #7 as a series involving a number of 
orders of approximation and to treat these by 
the methods of perturbation theory. In zeroth 
approximation the vibration-rotation energy 
levels of an axially symmetrical molecule, such 
as ammonia, are given by the expression, 


Wo/he= >i wil 
+A(P+J)+ (C.—A.)K?. 


Here the w; are the normal vibration frequencies 
(expressed in cm~') and d; are the degrees of 
degeneracy of the various normal vibrations. 


A,=h/8r'cI 4°, 


where J4° and J¢° are the equilibrium values of 
the moments of inertia for axes perpendicular to 
and parallel to the symmetry axis, respectively. 

When the change of the electric moment due 
to the vibration lies along the symmetry axis, 
the selection rules for the rotational numbers 
are AJ=0, +1 and AK =0. As is well known this 
produces a simple band of the so-called parallel 
type which consists of a positive, negative and 
zero branch. The lines of the positive and nega- 
tive branches have the spacing 2A, and are 
multiple, that is, each line is composed of a 
number of superimposed components. It is to be 
expected that in the next order of approximation 
these component lines may separate slightly 
from each other. 

In the case of the ammonia molecule there is 
one further circumstance to be recalled. The 


‘For a review of the ammonia spectrum see D. M.~ 
Dennison, Rev. Mod. Phys. 12, 175 (1940). 
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nitrogen atom possesses two equivalent positions 
of equilibrium, one on either side of the plane of 
the hydrogens. This fact produces a splitting of 
each energy level into two levels. The doublet 
separation is small (about 0.66 cm~') for the 
ground state but becomes rapidly larger for the 
excited states of the vibration »;. The wave 
functions characterizing these two states are 
either symmetrical or antisymmetrical with re- 
spect to the plane of the hydrogen atoms, the 
lower of the two levels of the doublet being sym- 
metrical. Parallel type bands, such as those to 
be discussed here, correspond to transitions con- 
necting a symmetrical with an antisymmetrical 
level. 

The methods for treating the interaction be- 
tween the vibration and rotation of polyatomic 
molecules have been discussed by Wilson and 
Howard? and have recently been applied in 
detail to the pyramidal XY; molecule by Shaffer.’ 
This treatment leads to an expression of the 
following type for the energy of a level not in- 
volving a perpendicular vibration.‘ 


/he= 
(1) 


The first term Wy* represents the vibrational 
energy. The quantities Ay* and Cy* are functions 
of the vibrational numbers, in first approximation 
linear functions. The superscript s indicates that 
the level in question is the lower member of the 
doublet pair caused by the double minimum 
character of the potential; the upper level will 
be designated with the superscript a. Shaffer's 
formula does not contain the superscript s or a, 
of course, since he did not consider effects 
arising from the double minimum potential. 
The fact that Ay* and Ay* may differ from 
each other is indicated by the following qualita- 
tive argument. Consider the case where K =0. 
This corresponds to the classical motion in 
which the molecule is rotating about an axis 
perpendicular to the axis of the pyramid. The 
centrifugal force tends to force the hydrogen 
atoms away from the nitrogen and to make the 
apex angle of the pyramid more acute. In other 


2 E. B. Wilson and J. B. Howard, J. Chem. Phys. 4, 260 
(1936). 

3 W. H. Shaffer, J. Chem. Phys. 9, 607 (1941). 

‘In the case of the perpendicular frequencies there occur 
additional terms due to the Coriolis forces. 


words the centrifugal distortion raises the barrier 
separating the two equivalent equilibrium posi- 
tions of the nitrogen atom. Accordingly, the 
doublet separation, which is a very sensitive 
function of the barrier height, will be decreased 
and therefore Ay*— A y* must be negative. On the 
other hand when K=J the energy levels corre- 
spond to the classical motion of rotation about 
the axis of the pyramid. The centrifugal forces 
will clearly tend to make the apex angle of the 
pyramid more obtuse and hence in effect will 
lower the barrier height. Thus we may expect 
Cy*—Cy* to be positive. It will appear that the 
experimental data are in accord with these 
predictions. 

The last term in Eq. (1), namely F(JK), 
represents what may be called the centrifugal 
distortion of the molecule. It is small, depends 
principally upon quartic combinations of the 
quantum numbers J and K and is in first approxi- 
mation independent of the vibrational state; 
hence no superscript s or @ is necessary. F(JK) 
has been calculated by Slawsky and Dennison.® 
Recently Foley and Randall® have measured 
under high dispersion lines of the pure rotation 
spectrum of ammonia. They were thus able to 
determine F(J, K)—F(J—1, K) experimentally 
and their results are in almost perfect agreement 
with the calculated formula, particularly when 
the new values of the molecular force constants‘ 
are used. 


F(J, K)— F(J—1, K) = —0.00294J% 
+0.00279J7K*. (2) 


The frequencies of the lines constituting a 
parallel type band may be written out easily 
from the energy expression Eq. (1). For example, 
the positive branch lines of the band correspond- 
ing to the change 0,—1, in the vibrational 
number are given by 


0,3, J-1K 


+(C\°— Co" K* 

+F(J,K)—F(J-1,K). (3) 

son ZidggawokY and D. M. Dennison, J. Chem. Phys. 7, 


*H. M. Foley and H. M. Randall, Phys. Rev. 59, 171 
(1941). 
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Fic. 1. Curves A and B: The double 10y band with partially resolved Q branches. The 0,1. component with center 
at 968.08 cm is marked A. Center of the 0.—1, band lies at 931.58 cm~. Curve C: A part of the P branch of the 16u 


band 1,—2.. 


Thus the fact that the quantities A and C are 
functions of the vibrational state introduces a 
convergence in the series through the term 
proportional to J?. Moreover since K may take 
on the values, 0, 1, ---, J—1, each line which 
was formerly single will be split into a multiplet 
of J components. The spacing of the lines in a 
multiplet will be quadratic. Since, as has been 
pointed out, the A and C depend upon whether 
the superscript is s or a, the convergence and 
the multiplet spacing of the complementary 
vibrational band 0,—+1, may be very different 
from that of the band 0,—1,. The intensities 
of lines in a multiplet are interesting. Because of 
the proton spin, when K is a multiple of three 
the line in question is enhanced by the factor 
two. The state where K=0, J is even and the 
vibration superscript is s does not exist and 
consequently lines involving a transition to it 
are absent. Likewise, states where K=0, J is 
odd and the superscript is a do not occur. 

The bands which have been chosen for the 
present investigation are the long wave parallel 
bands corresponding to the following changes in 
the vibrational quantum number V3, 0,—1,, 
0.1, and 1,—2,. The first two bands lie at 
about 10u while the third has its center at about 
16u. The next section describes the experimental 
procedure and results, and shows that the ob- 
served lines may be adequately represented by 
means of a formula of the type of Eq. (3). In the 
last section an attempt is made to calculate 
(Ay*—Ary’*) and (Cy*—Cy*) theoretically. 


EXPERIMENTAL’ 


The spectrometer employed for these observa- 
tions has already been described.* For the 10u 


7 During the spring of 1936, Dr. M. V. Migeotte pomp | 
in this laboratory mapped a part of the positive branch 


band an absorption cell 12 cm long with KBr 
windows was found convenient. The gas pressures 
ranged from 2 to 5 cm. The resolving power 
was about 1000 with the slit widths used (about 
0.45 mm, corresponding to 0.35 cm~'). For the 
weaker 16u band, a cell having an absorption 
path of two meters was available. Because of its 
small aperture much wider slits (about 0.9 mm) 
were necessary, and the resolving power was 
somewhat less. 

The results are shown in Fig. 1. The upper 
curve is the double 10u band, its high frequency 
component, namely 0,—1,, being indicated with 
A. The lower curve is the negative or P branch of 
the 16u band (1,—2, transition). It is consider- 
ably obstructed by many water lines. The posi- 
tive or R branch of this band is not accessible on 
account of the strong absorption of atmospheric 
CO:. The numbers below the lines in the figure 
are the initial J values, and those above the lines 
are the K values. The line frequencies, measured 
in wave numbers reduced to vacuum, are given 
in Table I. 

A study of Fig. 1 reveals that the 104 band 
possesses all the qualitative features predicted 
by Eq. (3). The lines of this parallel type band 
are clearly multiple and the multiplet spacing is 
very different for the band 0,—1,, where it is 
large, and for the band 0,—1, where it is so 
small as to be observable only in the lines of 
higher ordinal number. The enhancement of the 
lines where K is a multiple of three is also clear, 
but the absence of half of the lines where K =0 is 
less apparent since the resolution of the spec- 


the 10u band, and was able to resolve the fine structure of a 
number of the lines. His results were in substantial agree- 
ment with those we have obtained. We wish to acknowledge 
his priority in resolving this portion of the band, and to 
thank him for allowing us to examine his measurements. 

8 J. D. Hardy, Phys. Rev. 38, 2162 (1931). 
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trometer is not sufficient to distinguish between 
the lines K=0 and K=1. However, a study of 
the positions of the lines shows without question 
that the theoretically predicted missing lines are 
indeed absent. Specifically, the separation be- 
tween the line K=2 and the first line of the 
multiplet, which is either the single line K=1 or 
the merged line of K=0 and K =1, is systemati- 
cally larger in the latter case. 

‘The quantitative treatment of the data is not 
difficult. It is desired to show (a) whether the 
energy levels may be expressed in the form given 
by Eq. (1), and (b) to evaluate the constants 
Ay*, Ay*, Cy* and Cy*. The first step is to elimi- 
nate the centrifugal distortion term F(JK). 
This may be done by subtracting the numerical 
values given by formula (2) from the observed 
line positions of the positive branch and adding 
them to the line positions of the negative branch. 
One thus obtains lines corresponding to a ficti- 
tious molecule which possesses a vibration rota- 
tion energy expression in which the centrifugal 
distortion term F(JK) has been omitted. The 
second step is to use the combination relation to 
obtain the difference in energy between two 
states of the molecule having the same vibra- 
tional excitation, 1, or 1,, the same K but with 
values of J differing by two units. From these 


TABLE I. Observed positions of lines in the bands at 10u 


and 
R BRANCH P BrRancu 
TRANSITION J K ’vac P K "vac 
O.—1s 5-6 3 1053.11 2-1 0.1 892.05 
4 1051.42 3-2 1 872.52 
4-5 3 1032.05 2 871.70 
3-4 1 1013.06 4—3 0,1 853.66 
2 1012.47 2 852.70 
2-3 0,1 992.55 3 851.30 
5—4 1 834.73 
2 833.96 
3 832.61 
4 830.60 
6-5 0,1 816.48 
2 815.55 
3 814.13 
4 812.53 
5 809.68 
5-6 1084.59 1-0 948.21 
4-5 1065.49 3-2 908.15 
3—4 1046.34 4-3 887.93 
1-2 1007.51 5-4 867.80 
9-8 6 783.39 
7 782.75 
8 782.14 
10-9 6 763.55 
7 762.87 
9 761.35 
la—2. 2-1 0,1 590.32 5—4 1 536.14 
3-2 1 571.41 2 534.74 
4-3 0,1 553.96 3 531.99 
2 551.93 4 528.49 
6-5 3 $15.38 


TABLE II. Values of constants A. 


J K 2A J 
5—>6, 5-4 3 20.111 | 5—>4 19.780 
4 20.121 | 4-+5, 4-3 19.772 
4-+5, 4-3 3 20.120 | 3-+2 19.773 
1 20.113 | 1+2, 1-0 19.773 
2 20.137 


2-1 0,1 20.116 
Average 20.120 


Average 19.774 


2A;° or 2A;* are found by dividing by (2/+1). 
The results of this procedure are shown in 
Table II. Note that in the case of the band 
0,—1,, which determines 2A,*, the multiplet 
spacing is so small that up to J=6 the lines are 
observed as single. It is seen that the various 
determinations of 2A ,* and 2A;* are quite con- 
sistent although the last figure in the average 
values is no doubt meaningless. 

In a similar way the combination relations 
may be used to determine the 2A and 2A 
these are given in Table III. 

The spacings of the lines composing a given 
multiplet furnish the means of determining such 
quantities as The pro- 
cedure of course is to divide by 2K+1 the 
interval between two lines having the same J 
but with K values differing by unity. The 
results appear in Table IV. 

For the 16x band which corresponds to the 
vibrational transition 1,—2,, it is not possible to 
use combination relations since only one branch 
could be measured experimentally. However, 
since the energies of the 1, levels are known, 
those of the 2, levels may be found. From them 
the following numerical values are obtained: 


2A 2* = 20.36, 


Finally the centers of each of the three bands 
may be located with considerable precision. The 
interval between the ground levels 0, and 0, is 
taken from observations of Wright and Randall.® 


Wo'/he 0, 
Wo*/hc=0.66, 
= 932.24, 
W,2/he = 968.08, 
W2!/he= 1597.4. 


*N. Wright and H. M. Randall, Phys. Rev. 44, 391 
(1933). 
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TABLE III. Values of constants Ao. 


J K 2Ao* J 2Ao° 
4-5 3 19.876 | 5—>4, 19.887 
3-4 1 19.873 | 3-—»2, 19.890 

2 19.884 
2-3 0,1 19.877 Average 19.888 
Average 19.877 


It will appear in the next section that certain 
combinations of A and C are susceptible of 
theoretical estimation. For convenience we list 
these below. The data do not seem to warrant 
the retention of more than two significant 
figures. 


Ao*—Ao' = —0.00;, 

A,;*—A;*= —0.17, 

=0.232(C,*— Ci") — A). 


In the next section it will be shown that there 
are theoretical arguments for supposing that 
both Ao*—Ao* and Co*—Co* are small. Experi- 
mentally Ao?— A," is indeed small and it is safe 
to conclude that Co*—C,* is also negligible. 
Hence the approximation in the last equation is 
probably quite good and C,*— C,* = +0.06. 

The experimental data and their analysis 
appear to be very satisfactory. The consistency 
shown in Tables II, III, and IV demonstrates 
that these vibration rotation energy levels of 
the ammonia molecule may be expressed by a 
formula of the type of Eq. (1). Moreover the 
qualitative predictions regarding the effect of 
centrifugal forces upon the rotational constants 
A and C are verified since A,;*— A,’ is observed 
to be negative, whereas C,*— C;' is positive. 


DIscussION 


In considering the interaction between vibra- 
tion and rotation it is best to divide the effects 
into two classes: (a) the effects which are inde- 
pendent of the double minimum nature of the 
potential energy, and (b) effects which may be 
attributed solely to the double minimum poten- 
tial. Of these, the former will determine the 
average values of A and C (averaged over the 
two levels caused by the double minimum), that 
is, }(Ay*+Ay*) and 3(Cy'+Cy*). The latter will 
fix the interval between the two levels of any 


pair resulting from the double minimum. From 
Eq. (1), this separation, Av, is given by 


Av/he= (Wy*— Wr") /he 
+ (Cy*— Cy) K?. 


(a) The methods of treating the vibration- 
rotation interactions which are independent of 
the double minimum potential are well known 
and, as already stated, have been applied in 
detail to the pyramidal molecule XY; by Shaffer.’ 
They consist in developing the Wilson and 
Howard form of the Hamiltonian in a series of 
approximations and of treating each term by 
perturbation theory. This procedure does not 
distinguish between the levels s and a and 
consequently will describe the average values of 
the A and C for these levels. Shaffer finds, 
essentially, 


3(Ay*+A v’) =A,.-> a;( Vi+3d,), 
3(Cy°+Cv2) =C.— Vit3di). 


The subscript 7 ranges over all the normal 
frequencies and the V; are the corresponding 
quantum numbers. The quantities a; and y; are 
functions of the masses, the normal frequencies 
and of those anharmonic terms in the potential 
which are cubic in the normal coordinates. There 
are eight of these unknown cubic constants and 
it is thus evident that the fine structures of at 
least eight of the ammonia bands must be 
measured before this part of the problem can be 
solved in detail. Before turning our attention to 
(b) we may remark that the average values of 
Ay and Cy depend essentially upon the potential 
function in the neighborhood of the equilibrium 
configuration, that is upon the region where the 
wave function is large. On the other hand, the 
doublet separation is determined by the behavior 
of the potential near the barrier. A study of 
the doublet separation is thus capable of furnish- 
ing information about the molecule at points 
which are distant from the equilibrium con- 
figuration. 

(b) Since the Hamiltonian function proposed 
by Wilson and Howard is exact, it might furnish 
a starting point for the calculation of Ay*—Ay* 
and Cy*—Cy*. It is very complicated however 
and the particular series development used in (a) 
would probably not be applicable here. Certain 
approximations may be made which greatly 
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simplify the work and which appear to be well 
justified physically. The principal assumption is 
that the normal coordinate x3 associated with 
the band v; at 10u, leads directly over the lowest 
part of the pass connecting the two equivalent 
equilibrium configurations of the molecule. This 
assumption must be very nearly true since it is 
known experimentally that the doublet separa- 
tions for the excited states of the vibrations 7, v2 
and v4 differ only slightly, if at all, from the 
doublet separation 0.66 cm— of the ground state. 
On the other hand, the doublet separation of the 
first excited state of v3, is 35.9 cm~'!, about 50 
times that of the ground state. We may conclude 
that the motions associated with x1, x2 and x, do 
not lead from one equilibrium configuration to 
the other except by way of a very high potential 
barrier. This is also a plausible conclusion when 
the actual motions associated with these coordi- 
nates are considered. 

Since the subject of our study is the doublet 
separation and since this depends almost wholly 
upon x; we shall ignore the remaining coordi- 
nates. The actual Hamiltonian will be approxi- 
mated by one which is a function of x3 alone, 
the other coordinates being given their equi- 
librium values, namely zero. The normal coordi- 
nate x3 describes a deformation in which the 
hydrogen atoms move symmetrically with re- 
spect to the axis of the pyramid, and in first 
approximation the N—H bond distances are 
unchanged. Thus, whatever the value of x3, the 
molecule retains its axial symmetry. The mo- 
ments of inertia are of course functions of x3 but 
the products of inertia are always zero. Under 
these assumptions, the Hamiltonian takes on a 
very simple form. We express it in terms of y 
rather than of x; where y is the actual displace- 
ment of each hydrogen relative to the nitrogen 
atom, taken, however, along the path of the 
normal coordinate x3. 

V. 
P,, Ps and P, are the components of the total 
angular momentum of the system. The combina- 
tions occurring here may be diagonalized at once 
leading to a Hamiltonian which contains only 
the vibrational coordinate y and its conjugate 
momentum p. 


H= p?/2p+ 


It will prove convenient to rearrange the 
terms in the Hamiltonian; let, 


H= p?/2p+V+6V+ We’, 


where W,’° is a constant equal to A,(J?+J—K°*) 
+C,.K? and 


I—K?) (1/224 —1/2T4% 


The problem has now been brought to such a 
form that it differs from the usual double mini- 
mum problem in only two respects, (1) the 
potential is increased by 6V and (2) the energy 
constant of the whole system is increased by W,’°. 
This last point will, of course, have no effect 
upon the quantity to be calculated, namely the 
splitting of the levels caused by the double 
minimum potential. 

This splitting has been calculated by Dennison 
and Uhlenbeck'® using the W-K-B approxi- 
mation. 

They find 

A=2hv/A?, 
where 


(2/0) f My] 
0 


The origin for the coordinate y is chosen to 
be the equatorial plane through the nitrogen 
atom, that is, at the center of the barrier sepa- 
rating the two equilibrium positions for the 
hydrogens. y, is the first root of V-E=0. 

In the present calculation V must be increased 
by 6V and E by 6E, the latter to take account of 


TABLE IV. Values of constants. 


# K Cy —Co* J K Cit —Ai* +A’ 
4 — 0.300 10-9 7 —0.020 
3 —0.213 8 —0.024 
2 — 0.266 9-8 7 —0.016 
1 —0.240 6 —0.024 
5—4 3 —0.273 
2 — 0.256 Average —0.021 
1 —0.243 
43 2 —0.270 
1 —0.251 
3—2 1 — 0.266 
3 —0.258 
34 1 —0.210 
Average —0.254 


10). M. Dennison and G. E. Uhlenbeck, Phys. Rev. 41, 
313 (1932). 
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the fact that the average of the levels s and a 
may be slightly shifted by 6V. A* can be de- 
veloped since 5V is small. Thus A+6A=2hv/A? 
and 


A*=exp] (2/h) f 
0 


=exp(2/t)| f [2u( Jay 


+ f 8E)(2u(V—E) 


The first term in this exponent gives A while 
the next two determine 6A. The third term may 
be expected to be negligible in comparison with 
the second since the integrand is so small in 
the region y; to yi+éy:. A numerical calculation 
indicates that it contributes less than one per- 
cent to 6A and consequently we shall ignore it. 
If we assume that 6A is small, the exponential 
may be developed" and we obtain, 


6A = —(24/h) fiw V—8E)[2u(V—E) 


It does not appear to be possible to evaluate 
this integral by any means other than numerical 
quadrature. We have chosen to use for V the 
double minimum potential proposed by Man- 
ning.” 


V/he=a sech*(y/2p) sech*(y/2p). 


This function has the correct general form but 
there is no reason to suppose that it is the best 
function for describing the ammonia molecule. 
However, it is well known that the splitting of 
the levels is rather independent of the particular 
form of the potential. The constants obtained 
by Manning from the data for ammonia are 


1 The approximation of developing the exponential is 
only justified as long as the integral occurring in the 
exponent is small compared with unity. A substitution of 
the numerical values to be obtained later shows that this 
is true for the lines actually used in the analysis. For lines 
of large ordinal numbers where the multiplet spacings are 
correspondingly large, this will no longer be the case and 
it is probable that the energy levels cannot be expressed in 
the form of Eq. (1). < 

12M. F. Manning, J. Chem. Phys. 3, 136 (1935). 


a=66,551, b=109,619 and p=0.04793(h/cu)!. 
The mean energies E are not given directly by 
Manning but only the energy differences. A cal- 
culation yields as the value for the first excited 
state (E,°+E,*)/2hc = —43,692.0. 

The quantity 6V may be found by writing 
out the moments of inertia J, and J¢. These 
will be expressed first as functions of two vari- 
ables, 7, the distance between a hydrogen and 
the nitrogen atom, and a, the angle between a 
N—H bond and the equatorial plane through 
the nitrogen. 


I4=[3mM/(3m+ M) sin? a+ 3mr? cos? a, 
Ic¢=3mr* cos? a. 


The .notion corresponding to the frequency v3 
is, to . high degree of approximation, one in 
which r remains constant and a alone varies. 
The relationship between a@ and y is then, 
a=aoy/yo where yo is the distance along the arc 
from the equatorial plane to the equilibrium 
position and ap is the equilibrium value of the 
angle. From previous work! these are known to 
be, 0.38610-8 cm and 21° 49.7’, respectively. 
Actually the assumption that 7 is a constant is 
not quite justified and will be modified later. 
It is very instructive, however, to carry out the 
computation under this assumption and that 
will be done first. 5V is then determined as a 
function of y and the only quantity not yet 
discussed is 5E. 5E is the change in the average 
position of the levels s and a caused by the 
rotation. Since this change, for example the term 
3(Ay*+Ay*)—A,, is proportional to V+3, it is 
evident that for the first excited state, 


] 
X 
—}3(Co*+Co*) ]K*. 


Introducing the observed values we obtain, 
bE /hc =0.05(J?+ J — K*) —0.16K?. 


The contribution due to 6£ turns out to be 
6A/hc = 0.007 (J?+J) —0.029K?, which is con- 
siderably smaller than that due to 6V. 

The numerical evaluation of the integral deter- 
mining 6A for the first excited state is not 
difficult. Near the point y=y, the integrand 
becomes infinite but this portion could be ob- 
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tained through an exact integration where, how- 
ever, the potential was approximated by a 
straight line. Multiplying the integral by the 
observed splitting in the absence of rotational 
interaction, namely A/hc=35.8 we find, 


6A/he= —0.105(J?+J)+0.194K°. 


It will be recalled that the observed change of 
splitting is 


(6A/hc) obs = —0.17(J?+J) +0.23K?. 


Clearly this first theoretical estimate of 5A is 
only moderately satisfactory. It was based upon 
the assumption that the distance r between each 
hydrogen and the nitrogen atom remains con- 
stant throughout the motion. This point may 
now be reexamined. The normal coordinate 
treatment of the problem yields not only the 
frequencies but also the motion associated with 
each frequency. Thus it is easy to show from 
the normal frequency determinant that, in the 
motion v3, the hydrogen atoms oscillate along 
lines which make an angle Bo with the symmetry 
axis of the pyramid where 


cot Bo= V3(md3—b) /c. 


The quantities } and ¢ are potential constants 
the latest values of which are given in reference 1. 
A substitution yields 8)=16° 7.2’. The lines 
along which the hydrogens would travel under 
the assumption that r remains constant are 
easily found from the known equilibrium con- 
figuration of the ammonia molecule. These lines 
make an angle ap=21° 49.7’ with the axis of 
the pyramid. 

Thus as the hydrogen atoms move toward the 
equatorial plane through the nitrogen, the dis- 
tance r decreases and at such a rate that the 
angle between the line along which they move 
and the line perpendicular to the N—H bond is 
ay — Bo= 5° 42.5'=0.10 radian. If this motion 
were to continue from the equilibrium position 
all the way to the equatorial plane r should be 
set equal to 


rol 1 —0.10(a0— a) 1 —0.10a0(1 — y/o) ]. 


This last approximation depends upon the fact 
that the distance along the actual path is very 
nearly equal to the distance along the arc for 
such a small angle ao— fo. 


The spiral path postulated for the hydrogen 
particles has one unsatisfactory feature, namely 
that it makes an inward cusp at the equatorial 
plane. This may be eliminated by adding a 
further term to the expression for r; let, 


r=roL1—0.10(ao— a) +(0.05/a0)(ao—a)*] 
=r[1 


This function satisfies the condition imposed 
by the normal coordinate treatment, namely 
(dr/da)a—«,=0.10, and has no cusp at the equa- 
torial plane since (dr/da),.o=0. We are well 
aware that very little reliance can be placed 
upon this estimate for r. It must be emphasized 
that the normal coordinate treatment can only 
give information concerning the displacements 
near the equilibrium configuration. Our function 
contains this information and is otherwise merely 
plausible. 

The integral determining 6A may now be 
recomputed. It yields the result, 


(6A/hC) = —0.162(J?+ J) +0.222K?. 


The agreement between (6A) theo and (6A)obs is 
almost perfect and some remarks should be 
made concerning its significance. The change in 
the splitting of the levels is an effect which 
depends upon the potential energy function of 
the molecule at some distance from the equi- 
librium configuration, namely in the region near 
the potential barrier. From the motion near the 
equilibrium configuration we have estimated the 
manner in which the molecule most easily passes 
from one equilibrium position to the other, that 
is, we have estimated the distance r as a function 
of a. This estimate was admittedly very rough 
but the fact that it leads to a correct value of 6A, 
indicates that it must be of the correct order of 
magnitude. 

The agreement may be interpreted as furnish- 
ing information on the trustworthiness of the 
potential constants which have been used. It will 
be recalled that these were obtained! from esti- 
mates of the normal frequencies of the ammonia 
molecule and not, as in the usual procedure, from 
the positions of the fundamental bands. This 
latter method was employed by Slawsky and 
Dennison and led to rather different values for 
the constants. A substitution of the older con- 
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stants yields the angle B)=28° 40’ and thus 
would predict that as the hydrogens approach 
the equatorial plane, the distance r increases. 
This result would have destroyed the agreement 
between the calculated and observed separation 
of the double minimum levels. One further 
point may be noted. If the potential of the 
ammonia molecule was of the valence type and 


would not lead to the very satisfactory agreement 
which we have obtained. 

The change in the splitting of the ground level 
may be calculated in a similar fashion. It will 
be considerably smaller, partly because the 
integral is smaller but principally because the 
splitting itself Ao/hc =0.66 is so small. We find 


—0.0011, 


contained no cross product term between the 


change in the valence angles and the change in 
valence distance, the angle By may be computed The first of these values was determined experi- 


readily and is 22° 1’. This also predicts an in- mentally and found to be Ao*—Ao*= —0.005. 
crease of r, although only a slight one, as the The agreement is satisfactory although not very 
hydrogens approach the plane and consequently _ significant. 


Co* — = 0.0005. 
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The quantum-mechanical Hamiltonian function, H, for the general vibrating-rotating poly- 
atomic molecule of N atoms has been expanded by the method of Wilson and Howard to second 
order of approximation. It has been practicable to obtain solutions of the Schrédinger equation 
(H-—E)’=0 for asymmetric molecules, axially symmetric molecules and linear molecules. 
Expressions for the anharmonic terms occurring in the vibration energy, the effective moments 
of inertia, the antplitudes of the internal angular momentum of oscillation and the centrifugal 
distortion coefficients are derived so that when the normal coordinates are known the vibration- 
rotation energies of a polyatomic molecule may be calculated from the results given. Tetra- 
hedrally symmetric molecules and models in which internal rotation occurs are regarded as 
anomalous cases and are not treated. 


I. INTRODUCTION 


ECENTLY a number of papers have appeared in which the form of the vibration-rotation 
energies of certain polyatomic molecular models! have been derived to a second order of 
approximation. The method followed has in each instance been equivalent to that delineated by 
Wilson and Howard,? although the details have varied slightly from case to case. The method consists 
essentially in approximating by the method of the perturbation theory to three items, namely: the 
energy of the atomic nuclei oscillating anharmonically about their positions of equilibrium in a set 
of body-fixed coordinates; the effective moments of inertia of the molecule regarded as a semi-rigid 
rotator and the distortion of the molecular energies due to the centrifugal forces. 
An inspection of the final results in these papers reveals that they are nearly always of the same 
form. Thus, for example, the oscillational energy always consists of a set of terms linear in the 
1A. Adel and D. M. Dennison, Phys. Rev. 43, 716 (1933); W. H. Shaffer and H. H. Nielsen, Phys. Rev. 56, 188 (1939); 
W. H. Shaffer, H. H. Nielsen, and L. H. Thomas, Phys. Rev. 56, 895 (1939); W. H. Shaffer, H. H. Nielsen, and L. H. 
Thomas, Phys. Rev. 56, 1051 (1939); B. T. Darling and D. M. Dennison, Phys. Rev. 57, 128 (1940); Ta-You Wu, J. 
Chem. Phys. 8, 489 (1940); S. Silver and W. H. Shaffer, J. Chem. Phys. 9, 599 (1941); W. H. Shaffer, J. Chem. Phys. 


9, 607 (1941); W. H. Shaffer and A. H. Nielsen, J. Chem. Phys. 9, 847 (1941). 
2 E. B. Wilson, Jr., and J. B. Howard, J. Chem. Phys. 4, 262 (1936). 
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vibration quantum numbers, the coefficients of which are the normal frequencies; a set of terms 
quadratic in the vibration quantum numbers, the coefficients of which depend in a complicated 
fashion on the normal frequencies, the non-vanishing constants occurring in the cubic and quartic 
parts of the potential energy function and upon the dimensions of the molecule. Also the effective 
moments of inertia are in each case linear functions of the vibration quantum numbers and depend 
in an involved manner upon the size and the shape of the molecule, the normal frequencies and the 
coefficients of the anharmonic part of the potential energy which is cubic in the coordinates. The 
coefficients of the centrifugal distortion are to second order of approximation always independent of 
the quantum numbers and depend only upon the normal frequencies and the dimensions of the 
model. One may summarize this by stating that to a second order of approximation the vibration- 
rotation energies of a molecule depend always only upon the following items: the normal frequencies, 
the size and shape of the molecule and the non-vanishing coefficients of the anharmonic terms in the 
potential energy cubic and quartic in the normal coordinates. 

This suggests that it should be practicable to set up and solve the problem of the vibration- 
rotation energies of a perfectly general molecule of N atoms oscillating and rotating in space. 
Precisely as before the final formulation should consist of three parts; the vibration energy, the 
energy of a semi-rigid rotator and the energy due to the distortion of the molecule by the centrifugal 
forces. To obtain the energies in a particular case it would then be necessary only to solve the 
problem of the normal frequencies and determine from the symmetry of the model which of the 
coefficients of the anharmonic terms in the potential energy function are non-vanishing and then 
insert these results into the general formulation. 

We have undertaken to do this and have found that the problem of determining the vibration- 
rotation energies of the general molecule is in reality much simpler to carry out in detail than that 
of the special cases. In the following sections we shall set down the results of this investigation. It 
will be found that it is far less time consuming to arrive at the vibration-rotation energies of a par- 
ticular molecule by the general method than it is to derive them independently for a particular model. 


II. THe QUANTUM-MECHANICAL HAMILTONIAN 


The problem of determining the vibration-rotation energies of a polyatomic molecule reduces to 
that of finding the solutions and the eigenvalues of the quantum-mechanical equation 


(11— E)y =0, (1) 


characteristic of the motion in space of the molecule. 

In discussing the motion of a molecule it is convenient to make use of a set of coordinates x, y 
and z fixed in the molecule and a set of space fixed coordinates x’, y’ and z’. The origin of the two 
sets of coordinates shall be the same and they shall be related to each other by the Eulerian angles 
6, ¢g and y, defined as follows: @ is the angle between the fixed 2’ axis and the axis z; ¢, the angle 
between the x axis and the nodal line and y the angle between the nodal line and the x’ axis. The 
general polyatomic molecule will contain N atomic nuclei all of which may have different masses. 


The kinetic energy of these nuclei when expressed in terms of the coordinates x’, y’ and 3’ fixed in 

space will be - 

D Mie ty +21). (2) 
t=1 


The molecule is so oriented in the coordinates x, y, and z that the equilibrium position of the ith 


particle will be x, y and z, subject to the restrictions: 
i=N i=N i=N 
i=l t=1 i=1 

and t=N 00 i=N °° i=N 00 

M x,y; =0, M x,2,;=0, My 2;=0. (3b) 


i 
> 
i=l i=l i=1 
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The first condition states that the center of gravity of the molecule coincides with the origin of the 
coordinates x, y, when the atomic nuclei are in their equilibrium positions. The second condition 
states that the axes x, y, and z coincide with the principal axes of inertia of the molecule. The actual 
positions in the body fixed axes of the ith particle will be designated by x;, y; and 2; where 


xj;=x,+4x,, w= yt dy;, etc., 6x;, dy; and 62; being the displacements of the ith nucleus from its 
position of equilibrium. Expressed in these coordinates the kinetic energy may be written: 


i=N 
i=l a=z,y,z B=z,y,z 


where Jaa are the instantaneous moments of inertia* }> MA(g:+7.) and Jas are the instantaneous 
i=N 
products of inertia which by virtue of Eq. (3b) become equal to >} M. (Byb0r,-+.02;68;+ 5a,68;). In Eq. 
t=1 
(4) wa are the components of the angular velocity along the x, y and z axes‘ and Q, take the form: 
i=N 
Qa= 
i=1 
The first set of terms in (4) represents the kinetic energy of the atomic nuclei oscillating about their 
positions of equilibrium in the force field of the electrons. The form of this force field will depend 
upon the nature of the assumptions made concerning it. For our purpose we shall take for the 
quadratic part of the potential energy a function which may be regarded as perfectly general in the 


coordinates: 
Vo=3 K ij, (5) 


where Ajj, mn may be regarded as force constants and the 6s;; is the variation of the distance between 

the atomic nuclei 7 and 7. The displacements 6s;; can be expressed in terms of the coordinates x;, 
yi; and z; so that (5) becomes: 

Vo(6s;;) = Vo(xi, Vis 3;). (6) 

The coordinates x;, y; and 2; which occur in (4) and (6) are not all independent of each other, but 
are subject to the restrictions: 


i-N 
> AMWia;=0 (@ taking the values x, y and 2z), (7) 
t=1 

i=N 

M — =0 (a and taking the values x, y and z but (8) 

t=1 


the first of which states that the center of gravity of the molecule remains at rest in these coor- 
dinates and the second of which stipulates that to zero order of approximation the internal angular 
momentum is equal to zero. With the aid of (7) and (8) it is possible to determine the normal coor- 


dinates of the system. Of these there will be 3N—6, one normal coordinate, Qi, for each normal 
frequency w,. Not infrequently, however, a molecule will have so high a degree of symmetry that 
in certain of its modes it will be oscillating in force fields which are isotropic in two or three dimen- 
sions. In such cases two- and threefold vibrational degeneracies will arise, i.e., there will be 


3 Certain molecular spectroscopists [see for example R. S. Mulliken, Phys. Rev. 59, 873 (1941)] have found it con- 
venient to designate the moments of inertia of the molecule in ascending order of magnitude by Ja, J» and J-. In a dis- 
cussion of the vibration-rotation energies of a general polyatomic molecule it is inconvenient to assign definitely any one 
of the moments of inertia as the largest or the smallest. In this work the z axis will be taken as the approximately unique 
axis of the rotator. In some cases (for example, the methyl halide molecules or the formaldehyde molecule) J,, will 
therefore be the smallest moment of inertia, i.e., Ja; in other cases (for example, the BF; molecule or the H2Se molecule) 
Tes — be the largest moment of inertia, i.e., J,. It is probably always convenient to identify the moment of inertia J,, 
with 

‘ Explicit expressions for wa may be found in books dealing with classical mechanics. 

5 See, for example, E. T. Whittaker, Analytical Dynamics, fourth edition, p. 177. 
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respectively, two and three coordinates Q’ associated with the same frequency w. It is convenient, 


for this reason, to denote a normal coordinate by Q’,, rather than simply by Q/, where s will denote 
the particular frequency to be associated with the coordinate and o@ will take the values 1; 1, 2 and 
1, 2, 3 depending upon whether the frequency is non-degenerate, two- or threefold degenerate, 
respectively. Thus, for example, the non-degenerate frequency of a harmonic oscillator oscillating 


with a frequency w, will have associated with it the coordinate Qj, 1, while the twofold degenerate 
frequency of a two-dimensionally isotropic oscillator oscillating with a frequency w» will have 


associated with it the two coordinates Q;,1 and Qj, », ete. 
It is always possible to express the 6x,, 6y; and 6z; in terms of the normal coordinates Q,, as follows: 


so that when these are substituted into the first term of (4) and into (6) we obtain for the kinetic 
and potential energies of the oscillational motion the following relations: 

where the y,, are the reduced masses and k, are the generalized forces constants, equal, respectively, 
to Aa C bdon which will be complicated relations of the K;;, mn. It is convenient to replace the coor- 


dinates Q’. by (use)~!Qse So that instead of (10) we shall have: 


2T=CQie and 2Vo=¥E (Av = (10’) 


where bice™ (tse) Vises etc. 
In terms of the coordinates Q,, the expressions for the cubic and the quartic terms in the poten- 


tial energy expansion may conveniently be taken to be: 


8,0 


8, oO o’ o’’ er 
The appropriate Hamiltonian form for the kinetic energy of the general molecule in the quantum 
mechanics will, according to Wilson and Howard,* be the following: 


T=} haPat} Patas pat} UD s, ot (12) 
where ha=} Ds | 2uasPst+(Psuas) +Masu'(Psu')} and a, 8 are summed over x, y and z. In Eq. (12), 
aa = and (4 (13) 
where Jig are respectively defined in the following manner: 
Korot =X Mile Mine) 
Quer 


* Darling and Dennison have given a form which is slightly different from that derived by Wilson and Howard. To 
second order of approximation, however, they are entirely equivalent. 
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and where 


Explicitly the quantities etc., will become: 
(e) 


oe 


so 


where 


eso= My M (silico +X Mise), 


s’o’ 


= M + > M; i(Nisolis’ — Nis’ 


(15) 


| 
Ty =u. 
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The quantities P, are the operators for the three components of the angular momentum referred 
to the moving coordinate system and they are defined as follows: 


P.=(—iy(sino)| (= cos) +sine= |sino)-, 

sind \dy 00 

P,=(-—th) (= cose) [sind (16) 
sind \dy 00 


P.= (—th)d/d¢, 
6, g and wy being the Eulerian angles referred to above. The operators p., py and p, which occur in 
(12) are the operators corresponding to the internal angular momentum due to oscillation and their 


definition will be given at a later point. 
It is convenient here to replace the Q,. by (h?/d.)'g.-. After considerable algebraic manipulation 
it is possible to set down in orders of magnitude the quantum-mechanical Hamiltonian for the general 


polyatomic molecule. It will be: 


(e) (e) (e) 


Ho=(h/2) + + + + (17a) 


(e) (e) (e) (e) (e) 


(e), 2 (e) (e) 


+ ) ) — (Ise )) 


(e) (e) (e) (e) 


—(fse(P2P:+P:P2)/ (Lez )) — )) (17b) 


— (3b — ) — sre + 2d Lee 
(e) (e) (e) (e) 


s’’a’’ 


(e) 


so s’’a’’ 


(e) (e), —2 


(e) (e) 


(Ased ) = (Deed /Tyy ) 


(e) (e) (e) 


(e) (e) (e), .—1 (e) 


(e) (e) (e), 


— (declarer — (PP: + (Ly ) (Lee 


(e) (e) (e) 


+ ) — ) — ) \(PyP2+P2Py)}. (17) 


| 

so 
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In terms of the coordinates g,, and their conjugate momenta p,,= —ih(0/0q,.<) the operators pa 
will be the following: 


owe 


p, and p, being obtained from (18) by permuting /;,., mise and ms. in cyclic rotation. 


III. THe ENERGY VALUES OF THE MOLECULE 


In this section we shall endeavor to obtain the eigenfunctions of the Hamiltonian, 7, which was 
arrived at in Section II. We shall approximate these by resorting to the methods of the perturbation 
theory of quantum mechanics proceeding after the manner of Shaffer, Nielsen and Thomas! in 
dealing with the vibration-rotation energy levels of the tetrahedrally symmetric XY, type of 


molecule. 
This method, recognizing that all terms in JZ; save the Coriolis terms arising from the degenerate 


oscillations can contribute first in second order of approximation, depends upon the effectiveness 

of carrying out a contact transformation on H, THT =H)+Hi+H:+--- such that H{ will to 

second order contain only the above Coriolis terms. It is possible always to find a linear combination 


of the zero order wave functions such that the matrix H; will contain elements only along the 
principal diagonal. In effect this reduces the problem of obtaining the second order energies to a 
first order perturbation calculation where the wave functions required are the stabilized wave 


functions of Hp, i.e., the linear combinations of the zero-order wave functions which will diagonalize 


We shall denote the transformation function e** by 7(A) which to second order of approximation 
may be written: 


(19) 


Transformation of the Hamiltonian J7=H)+H,+H2+--- in the manner indicated above leads to 


THT =H’ =Hi+Hi+Hit- -+, After equating like powers of \ we obtain: 


The portions of JZ, which it is desired to remove consist of terms each of which is a function of the 
normal coordinates q,. (or the conjugate momenta p,.) multiplied by a coefficient which either is a 


constant or a function of the operators P., for example aq, or QseP.. Such coefficients as occur in 
the latter example may also be treated as constants for to second order of approximation the error 
incurred because of the non-commutability of functions of P. with Hp will not contribute to the 
energy. The complete S function will consist of a sum of terms, each element being so chosen that 
it will remove a single term in //;. In orders of magnitude the transformed Hamiltonian H’ may be 


written: 
Hy=Ho, (21a) 


where r assumes those values of s for which the frequencies w, are twofold degenerate; 
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(e), 4 (e), 2, (e), 2 


(e), 2, 2 


(e).2 


+3 > ) )L(Sa sok sss/Xs) + +L MM — Mis 
2 


+ (bere M (Nisel iste: —Nis'o'lise) vol it — Nir se) 


X (dst Der) EL! (h/M (12) + 


so 
2 


so 


so 


{ (2A /h?*) —QseP | 


where Af contains terms which can contribute in higher approximations only. A prime after a sum- 
mation sign indicates that the indices can here not assume the same values as those which precede 
them ; a double prime following a summation sign preceded by a summation sign with a single prime 
indicates that the third set of indices do not take the same values as the first two sets and that the 
second set do not take the values of the first set. 

Since the Hamiltonian is entirely general it is possible in principle to proceed to find the eigen- 
values for any model in which there may be non-degenerate, twofold degenerate and even threefold 
degenerate oscillations. In practice, however, it seems advisable to exclude the latter class of molecules 
from this discussion because the stabilized wave functions in such cases are quite complicated to 
work with. Such molecules together with molecules in which internal rotation takes place or where 
other anomalies occur had better be treated as special cases. 

In the case of molecules where there are oscillations isotropic only in two dimensions there will 
be two coordinates gai and ga2 associated with a normal frequency wg. It has proved convenient when 
dealing with such oscillations to introduce the cylindrical polar coordinates pz and x, so that 
Jai = Pa ANd Ga2= Pa COSXa. When this is done the zero- and first-order Hamiltonian operators 
become: 


Ho=(h/2) +42) + — pr — ps 


(e) (e) (e) 


(22a) 


i i’ 
t 
| 
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(e) 


where di M ire —lirgmiri). 


The solution to the zero-order Schrédinger equation is known’ to be: 


where /7,,(g.) are the Hermite polynomials and R,,1,(p,) = pF (pr), F(p-) being the associated 
Laguerre polynomial. In the latter case /, may be regarded as a quantum number of internal angular 


momentum. It takes the values v,, v-—2, ---0, or 1. Ve in (23) is the characteristic function of the 


last part of (22a). The function Ve is most conveniently expressed as an expansion 


K, M) 


where ¥°(J, K, M) are the known solutions to the problem of the symmetric rotator which forms a 
complete orthogonal set of functions. In the above functions J, K and M are, respectively, the 
quantum numbers of total angular momentum, the quantum number of angular momentum about a 
body fixed axis and the so-called magnetic quantum number. 7 is an index number running from —J 


(1) 
to +J and Ase m are constants. 


To find the eigenvalues of /7) and the values of ee u it is easiest to follow Wang’ and set up the 
matrix Ho in terms of the basic wave functions ¥(v)¥°(J, K, M) where, of course, ¥°(J, K, M) will 
not in general be a solution to the rotation problem. The resulting matrix will in general not be a 
diagonal matrix, but will have the elements: 

e e 

(v, 1, J, M, K|Ho|0, 1, J, M, K)=h + 1) ] 

2.2 (e) (e) (e) 
+(K h /2){ (1/22 ) J}, (24) 

(v, 1, J, M, K| Hol, 1, J, M, — —K(K 1) (K £1)(K£2)}, 

where v and / are taken to embrace all the vibration quantum numbers and g, is the order of the 
degeneracy of the vibration frequency w,. An inspection of the elements (24) will reveal that for 
each value of the quantum number v and / the matrix will consist of steps, one step for each value of 
J. Since the elements (24) are independent of the quantum number M which may take all values 
from —J to +J it is evident that each step will be repeated 2/+1 times. The energies Eo(v, J, J, 7) 
for the molecule may now be had simply by diagonalizing Ho, or which is the equivalent, by setting 
equal to zero the secular determinant of this matrix | (K|/7Zo|K’)—Eéxx’|, 6xx- being the Kronecker 


symbol, and determining the roots. The function Ve is now uniquely determined also, for the coef- 


ficients Apne are simply the first minors of the above determinant. 
Using the same basic wave functions we may arrive at the first-order corrections to the energy 
matrix. The only non-vanishing elements are the following: 
2,,(e) 

(v, 1, J, M, K|Hi\», J, M, K)=X, , (25) 
which lie along the principal diagonal. Since 17; has elements only along the principal diagonal the 
second-order corrections can be determined by a first-order perturbation calculation using the wave 
functions (23). The function ¥(J, K, M), as we have seen is really not a solution to the rotation 


problem; hence H: will contribute also to the elements of H which do not lie along the principal 


7D. M. Dennison, Phys. Rev. 41, 304 (1932). 
8S. C. Wang, Phys. Rev. 34, 243 (1929). 
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diagonal. H} is, however, diagonal in all the quantum numbers except the quantum number K and 
the only non-vanishing elements are the elements (K|K), (K|K+2) and (K|K+4). When the zero-, 
first- and second-order contributions are combined one obtains for the elements of the entire energy 
matrix accurate to second order of approximation the following: 

(K|H|K) = 
(26) 
(K | |[f—K(K+1) J[f—(K+1)(K+2) J[f—(K+2)(K+3) J[f—(K+3)(K+4) ]}'Re, 
where 
Ro=(E,/he) +3(X Y.)I(J +1) Ri=2Z. 
R2=Z,—3(X Y,) —DyxJ(J+1), R3= — Dx, R= (1/4)(X,— Y,) —6,J(J+1), 


(e) (e) (e).2, (e).2 


(1/8w.) { ) — ) (Ise) (27) 


(e), 2, .(e), 2 


+2(bseCse— 2€se)/ ) ) a 


2 2.2 (e), 2, _(e). 2 


The quantities occurring in (27) are defined in the following manner: 


where 


X.=h/8e cle, Y.=h/8ecl, Z.=h/8x 


and where 
(ce), 2 (e), 2.2 (e), 2, 2 


Xo = (1/8we)[ ) ) ) ) ) 


(In) (Ise) J, 
48 (Iss) Use) Js 
2 (e).2, 2 2 (e), 2, _(e), 2 
16(feo/ ) (Ii: ) ) (Jez ) (28) 
(e) (e) (e) $208 


Xee=Xe (h/Aw ws) {LA sone — (Gee/ — — (fre/ Tea) ]— 


s’o’ 


X 
Yeo= VoL! cgelyy {[ Boose — (dee/ Ise) — ) — sz) J— 3b 


C Ws — > M —Nis'a'lise) My (Nir se” se) 


X 


| 
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s’a’ 


X 


(e), 2, 2 


Dx — (Ise) (Ise ) 


(e), 2, 2 


— (Ize) 


(e).4 (e).4 


The quantity EZ, in (27) is the vibration energy of the molecule and may be written: 


# r 


in which 


22 2 222 


Xs3° = ) {Ress's’ c Ws) eg c 


co’ (30) 


22 2,.(e) 


x> M ire —Lirgm iri) M — }, 


= (h/8x Is.) M (Liem ir2 —L irom ir) M Mm 
i 
(Eo/hc) is a constant term depending in a very complicated manner on the potential energy constants 
and since in practice one is only interested in differences between energies and not in their absolute 
values this term is not set down. 
As in the case of Ho, for a given value of the quantum numbers v and / the above matrix will 
consist of steps, one step for each value of J. As before, the energy is independent of the quantum 
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number M so that each step will be repeated 2/+1 times. For a given value of J the actual energy 
values may now be obtained as in the zero-order problem, i.e., by diagonalizing the matrix H. This 
again is accomplished by setting the secular determinant |(K|H|K’)—Eéxx:| equal to zero and 
solving for the roots. This determinant will resemble the Wang-Klein determinant for the asymmetric 
rotator, but will be considerably more complicated in general. The determinant can at once be 
factored into sub-steps corresponding to given values of the quantity }°,/, and each step will be 
repeated twice, once for the positive value and once for the negative value of >°,/,. Each of these 
sub-steps can further be factored simply by arranging the rows and columns so that the even values 
of K stand together and the odd values of K stand together. The elements of the Wang-Klein 
determinant are independent of the algebraic sign of K so that a further factorization is possible. 


( ) . . 
It may be seen that because of the term >-,(¢J,K/I my in the diagonal elements this is no longer 
true in the general case and no further factorization is at once possible. 
Darling and Dennison have pointed out that for a molecule like the H,O moleculate a relation 


A=I,, (Tes <I exists such that A is independent to this approximation of the 
anharmonic constants. That such a relation exists in general for all planar molecules may be verified 


from the relations (28) by setting z;=0 in which case @s.+)5. — Cs =0. Such relations will not exist in 
general for non-planar molecules. 


IV. THE VIBRATION-ROTATION ENERGIES IN SPECIAL TYPES OF MOLECULES 


A. Asymmetric molecules 
In the case of asymmetric molecules the atoms assume such geometric configurations that the three 


moments of inertia J,,, J,, and J,, are all different from each other. In such cases, however, it is 


true almost without exception that there are no degenerate frequencies. There will therefore be only 
a single coordinate g.; associated with a frequency w,, i.e., the index o takes only the value o=1. 


When this is so the terms which would occur in Hj will all vanish and hence the elements in the 
matrix H will become entirely independent of the sign of K. Except in cases where two vibration 
frequencies become accidentally degenerate the sub-steps of the secular determinant of the matrix 
H can be further factored as pointed out by Shaffer and Nielsen with reference to the triatomic 
non-linear molecule. When results of this kind are to be applied to experimental data it is frequently 
desirable to expand the secular determinant of the energies into a set of algebraic equations. This 
has been done by Nielsen for the triatomic non-linear molecule for the values of J up to and including 
J=6. It is found that these same equations are valid in the case of the general asymmetric molecule 
except that the more general definitions of the R; given in (26) must be used and that ¢ is here equal 


to (2E/h?— Ro). 


B. Axially symmetric molecules 
In molecules belonging to this important classification the geometric arrangement of the atoms is 


such that the two moments of inertia ae and oi are equal to each other while J ny will in general 
be different. In the case of such molecules the (K|K-+2) elements will always vanish except for 
second-order terms.’ The (K|K-+4) elements are also of the same order of magnitude. It is found 
that when the secular determinant for the energies is expanded that the elements off of the principal 
diagonal can contribute only in an approximation higher than the second. For our purpose they may 
therefore legitimately be neglected. This is equivalent to stating that to second order of approxima- 

®°H. H. Nielsen, Phys. Rev. 59, 565 (1941). 

1 Shaffer and Silver have found that in the case of the planar XY;, the pyramidal XY; type and the axially symmetric 


XY;3Z type molecules these elements vanish identically even to second order of approximation. Reference to the definitions 
of Xo, Yo, Xs¢ and yz¢ [Eq. (28) ] makes it seem plausible that to second order at least this is always true. 
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tion a molecule which is a symmetric rotator when the atoms are in their equilibrium positions will 
effectively remain a symmetric rotator also in higher vibration-rotation states. The vibration- 
rotation energies of such a molecule are then just the diagonal elements of (26) of the energy matrix. 
For such molecules it is convenient to state them in the following manner: 


(E/hce) =(E,/hce) + (Erot/he) (h cl (31) 
where (£,/hc) is the same as that given by (29) and (30) and where 
(Eyot/hc) = J (J +1)B,—K*(B,— C,) (32) 


in which 
B, =B,—ay- 


C,=C.-yo- (vstgs/2)ys, 


B.=(h/8r lec), (1 


2_(e) 


C.=(h/8x I: ¢), 

2(Xotyo), Yo=rZo, 

and where D,;, Dyx and Dx are defined as in (28). 


(33) 


C. The linear polyatomic molecule. 
Linear polyatomic molecules constitute a group which may be regarded as special cases of the 
axially symmetric molecules. In molecules of this type p:—P.=0, ie., p:=P.=(¥ih) O,1, and 


from this restriction it is easily shown that all terms involving J,,, és, and f,, will vanish from the 
Hamiltonian function, #7, and hence also from the elements of the energy matrix. The moments of 


inertia 84 and will be alike so that there will also here be no elements off of the principal 
diagonal. The vibration-rotation energies of linear polyatomic molecules may therefore conveniently 
be written: 

(E/he) = (Eyin/hc) + (Erot/he), (34) 


where 
(Evin/he) (Eo/hc) +>, +¢,/2) +>. (ve +>, (35) 
in which 
222 22 2 
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and where 
(Erot/he) =(J(J +1) —P]B.—[J(J +1) PD, (36) 
in which 
B,=B.->., (vetgs/2)as, (37) 
2 _(e) 


B.=h/8x , 


(e) 222 


As =B,(h/4e { [A soso 3d C Ws 


Mi (Mj — es) / } ’ 


(e).4 


V. CONCLUSION 


In the preceding sections relations have been derived for the vibration-rotation energies of a 
general polyatomic molecule to second order of approximation. These relations have been used to 
verify the results arrived at independently by various authors for certain molecular models and in 
each case have been found to corroborate their results. Enough cases have been tested by the author 
to verify the contributions from each kind of term in the Hamiltonian, H, and this is deemed a con- 
vincing check on the accuracy of the calculations here reported upon. In each of the cases tested by 
the author it was found that the vibration-rotation energies could be deduced in a small fraction 
of the time required to obtain them by the Howard and Wilson method starting at the beginning. 

It will be seen by inspection of the relations (26-30) that two and only two types of resonances 
occur in the energies to this approximation. These are of the Fermi-Dennison"™ type and of the 
Coriolis resonance type. The first is a resonance between a fundamental frequency w, and a com- 
bination frequency w,+,, i.e., the resonance denominator will contain (w-—w,—w,). Such terms 
will shift the centers of the vibrations bands, but will leave the spacings between rotation lines essen- 
tially unchanged. The second type is a resonance between two fundamental frequencies, e.g., 
(w,—w,;). This kind of term will alter the effective moments of inertia and may therefore give rise 
to anomolous spacings between the rotational lines. Essentially, however, it will not alter the 
positions of the band centers. Whenever the resonance denominators w,—w,—w; and w,—w,, re- 
spectively, become too small the usual methods of perturbation theory fail and other methods of 
procedure must be resorted to. While to this order of approximation these are the only resonances 
which can occur, resonances may, of course, also occur in higher approximations and may probably 
also make themselves felt in the spectrum. 

It is pertinent here to inquire into the practicability of carrying such calculations to still higher 
orders of approximation. It is doubtful whether such calculations would be particularly useful since 
in a higher approximation it appears that the secular determinant for the energies would contain 
elements of the kind (K|K), (K|K+1), (K|K+2), etc. This would mean that the secular deter- 
minant would not factor into sub-steps and that the work involved in determining the energies 
would become prohibitive. 

Moreover, in another approximation interactions between the electronic motions and the motion 
of vibration and rotation might be expected to set in and these certainly would be difficult to take 
satisfactorily into account. Also the approximation used throughout these and other similar calcu- 


1 E. Fermi, Zeits. f. Physik 71, 250 (1931). 
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lations that the oscillators are approximated by simple harmonic oscillators where the variables g,. 
extend from + * to — ~ rather than from go to + © would begin to become invalid in higher approxi- 
mations.” Finally the approximation made earlier that coefficients multiplying functions of qg,. in 
IT,, which themselves are functions of P,, may be treated as constant coefficients becomes inadequate. 
This convenient approximation is satisfactory so long as the energies are not estimated beyond the 
second order, but in still higher orders of approximation the non-commutative character of Hor and 
P., may become significant. It appears, therefore, that in general the work which would be required 
to extend calculations of this kind to approximations higher than the second would become pro- 
hibitive. 
APPENDIX 

It has seemed worth while to insert an example to illustrate how the results derived in the foregoing 
sections may be used to obtain the vibration-rotation energies in a particular case. The example 
chosen pertains to the planar XY; type molecule which already has been studied by Silver and 
Shaffer! (S. and S.). In their work the plane of the molecule is made to coincide with the xy plane 
of the body fixed axes so that the z axis becomes the axis of symmetry. The configuration is so 
oriented that the equilibrium positions of the Y particles and the X particle will be, respectively 
(—a, 0, 0); (a/2, —34a/2,0); (a/2, 34a/2, 0); (0, 0, 0). 

From the foregoing work it is possible to write the displacements of the particles X and Y along 
the x, y and z coordinates in terms of the normal coordinates Q,, as follows: 


where the miss and are the following: 
1% =(3M)-3, (m/(3M+m))'y—6], 


In=—3}(3M)-4, lon =(3M)-3 
6 


la=0, = Ly22=0, 
= L442 =0, 
Lis =123 =133 =143 = 0, (Al) 
m,=0, min =0, = min =0, 


ma=}(M)-, man = M222 >= M241 >= 


M242> (m/(3M-+m)) 


22 T am indebted to Professor L. H. Thomas for this information. . 
13 In the case of the frequencies a and w; which are non-degenerate the index o takes only the value o=1, It is con- 


venient in such cases to omit the index since no ambiguity arises out of so doing. 
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m4,=0, =0, = —m*(3M/(3M+m))y, Mai =0, 
M2 = 


M13 = M23 = M33 = M43 = 0,7 


M and m being the masses, respectively, of the Y and X particles and 6 and y being transformation 


constants. 
The first step in evaluating the energies is to determine the values of d4., dsc, etc. and A sess, Basses 


etc. We shall evaluate only one of these, namely, a2; and Ag222. The first of these is: 
=2(Myyminn+ + yma). 


Since m121=m42,=0 and we obtain: 
=2(M) 'y26= — (Ise) 6. (A2) 


In a similar manner we find that 


Cse=0(s #1). The constant A 2222 is defined by (15) and in our particular case it becomes: 


A 2222= 
— {(m/(3M+m))y+(3M/(3M+m))y}?= 8/2. (A3) 


Likewise one finds that A,,=B,,= }, A 2121 = A 2222 = = Bo222= 6*/2, A3333 = B3333 = 0, =A 


= Bagi Cu=1, Crose=0 (S¥1). 
We shall next compute the values of £2 and £4 which occur in the first order energies. By definition 


4 

f2= D0 Mi(lievmi22—li2emin). With the values of J;,¢ and mise given in (Al) one obtains at once 
i=1 

f2=(7?— 6). In the same way one verifies that {4= —fo. 

The vibration term values are given by the relations (29) where the constants x,, and x1, depend 
upon the coefficients occurring in the anharmonic part of the potential energy. We shall evaluate 
only two of these, namely, x3,4 and x43. Inspection of the anharmonic part of the potential energy 
function adopted by S. and S. shows that the only k’s which can occur in (30) are R4a33, Rags and iss. 


Evidently 


34> (h/ 32m Ras 1k133/' ws( [M (Mo41N23-+ 134133) 
13+ 23+ M3 42M 23) }, 


of which the last two terms reduce to 4x ¢ w36 . The term x4; is entirely like x3, except that the last 
term is equal to 4m c ws . The x34 and x43 are added together and if it is taken into account that 


= and "kiss =hcBiss in the notation of S. and S. we have at once that 


= (8 Ts.) In a similar fashion all the other Ga 


obtained by S. and S. may be derived. 


| 


810 HARALD H. NIELSEN 


We shall then proceed to consider the rotational constants. This consists in evaluating the B, 
and C, as defined by S. and S. We shall here evaluate only B, which is B,.—ao— > (vs+g;/2)as, Be 


being (h/8 Isc) since re =I. The constant ap will in our notation be xo=yo. From the relations 
(28) we have that: 


xo0= (—h'/256r'c ) (1/8) { [(dse } ’ (A4) 


which upon insertion of the values d.., bs, and d,, leads immediately to the result: 


Xo = Yo= Ao = (AS) 


We proceed to evaluate the a;. In this appendix we shall evaluate only a: which will serve to demon- 


strate how the others may be obtained. We have: 


X01 = (8/2 — 20. + (2) Use) 


and 
(e) 


2 2 2 
— (Mm Mmo22n23+ Mmg22M33 + mm 422M 43) *(3w2+01)/ w1) 
After inserting the values for m;,, and simplifying one obtains: 


= = [(—38'/2)+ (212) — 7 ], (A6) 


which is equivalent to the expression of S. and S. for a2. 

Of the centrifugal stretching coefficients D we shall evaluate only Ds which is defined in Eqs. (28). 
It here becomes: 
Dy= { (1/801) (16152 + (1/802) — + 8128) J+ 16 Lee) 


which becomes equal to 
/w2)+(y (A7) 
It is easy to verify the relationship for Dx obtained by S. and S. and also their relation that 
These constants being evaluated the vibration-rotation energies of the XY; planar molecule are 
determined and are given explicitly by the relations (31) and (32) in Section IV. 
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Increase of Residual Magnetism Caused by a Current Flowing Through an Iron Bar 
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If a bar of high grade wrought iron or of permalloy is placed in a sufficiently strong magnetic 
field, and if this field is abruptly reduced to zero, a current set up in the iron causes an increase 
of the residual magnetism. This effect is not obtained when the field is gradually eliminated, nor 
is it obtained when the current flows in a wire along the axis of a bored out iron cylinder unless 
the cylinder is slotted longitudinally. Then the increase of magnetism is again observed. There 
seem to be two effects involved in this phenomenon. One of these results in decreasing the 
residual magnetism in accordance with Wiedeman’s and Villari’s observations, and is apparently 
associated with the circular flux set up in the bar by a current sent directly through it or in a 
separate wire along its axis. The other effect, if it could be completely isolated, would probably 
produce only increases in residual magnetism, and seems to depend on an unstable condition 


of the magnet. 


INTRODUCTION 


T has long been known that a current of 
electricity sent through an iron bar or wire 
previously magnetized tends to reduce the re- 
sidual magnetism much like the effect of a me- 
chanical jar. In 1862 Wiedeman! described this 
decrease of residual magnetism, and also a less 
vigorous increase when the current through the 
iron was broken. In 1865 Villari? published the 
results of extensive investigations in the same 
field, and he too observed the decrease in residual 
magnetism caused by a longitudinal current, but 
does not seem to have observed the increase when 
the circuit was broken. A study of his tabulations 
shows that after the first ‘‘make’’ the galvanome- 
ter deflections (by which alterations of flux were 
observed) caused by subsequent and repeated 
“makes” and “breaks,” rapidly approach the 
same constant value though in opposite senses. 
This is obviously a purely inductive effect and 
has nothing to do with changes in the residual 
magnetism. So, although Villari reports galva- 
nometer deflections suggesting increased mag- 
netization, they are always subsequent to the 
initial ‘‘make,’’ while the significant deflection 
produced when the current was first established 
always indicates the expected decrease. 
The phenomena just cited seem to be closely 
associated with the effect of mechanical stresses 
of various kinds upon the magnetic properties of 


1G. Wiedeman, Poggendorff’s Annalen 117, 215 (1862). 
2 E. Villari, Poggendorff’s Annalen 126, No. 9, (1865). 


iron, but in general a jar or twist or pull tends to 
diminish residual magnetism. However in his 
collected Mathematical and Physical Papers, Vol. 
II, page 351, Sir William Thomson (Lord Kelvin) 
reports an increase in residual magnetism when 
tension was suddenly applied to a previously 
magnetized wire of soft iron, whereas a steel wire 
showed the usual decrease. In a footnote he says 
he later concluded that this could be explained by 
the influence of terrestrial magnetism, but as we 
found the same increase when a bar of wrought 
iron received an end on tap, it would seem as if 
the effect really existed in certain kinds of iron. 

During the summer of 1939 we undertook a 
further study of the effect of a current on residual 
magnetism, using at first a bar of permalloy 
given us by the Bell Telephone Laboratories, and 
later bars of very pure Norway iron, cast iron and 
nickel, the latter being the gift of the Inter- 
national Nickel Company. In the course of these 
observations we were very much surprised to find 
that when a cylindrical bar either of permalloy or 
of Norway iron had been magnetized beyond a 
certain minimum value by a solenoid whose cur- 
rent was then suddenly broken, an increase 
instead of a decrease of the residual magnetism 
was produced by a current sent through the iron 
in either direction. Two effects seem to be 
present. One of them which we shall call “‘effect 
A” operates with striking rapidity, and tends to 
increase the residual magnetism. The other 
“effect B’’ results in decreasing the residual 
magnetism and is definitely slower in its action. 
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Ke 
Fic. 1. Diagram of apparatus. S...... S—magnetizing 
solenoid; B...... B—bored out cylinder of Norway iron; 


C—galvanometer coil. 


The difference in the speed with which these two 
effects take place appears in borderline cases 
where both are of nearly equal magnitude. The 
galvanometer starts vigorously in the direction 
indicating increasing magnetization, but is then 
deflected with less violence though often farther 
in the opposite direction. With still larger residual 
magnetization, the positive throw, indicating an 
increase, completely dominates the negative 
throw and increases with increased magnetization 
and also with an increase in the longitudinal cur- 
rent through the bar. 


EXPERIMENTAL ARRANGEMENT 


The apparatus used in these experiments has 
gone through a variety of modifications, but the 
following arrangement was finally adopted as the 
most satisfactory one. The bar under investigation 
is a hollow cylinder about 20 cm long and 2.20 cm 


diameter with an axial hole 1.2 cm in diameter. In 
a second exactly similar cylinder cut from the 
same bar of Norway iron a longitudinal slot 
about 2 mm wide was cut through the entire 
length of the cylinder, thus making its section 
like a nearly closed letter C. This slot was 
designed to minimize circular magnetization and 
it resulted in a very marked reduction in the 
“B effect.’”” The ends of these hollow cylinders 
were covered with sheet copper caps tightly 
screwed down and the wires supplying the longi- 
tudinal current were soldered to the centers of 
these caps. In this way it was hoped that the 
current density throughout the bar’s section 
would be made fairly uniform. | 

The bar is magnetized by a double layer 
solenoid of 15 turns per cm having a mean diame- 
ter of 3.6 cm and a length of 76 cm. This is set up 
in line with the earth’s field, and Helmholz coils 
having a diameter of 60 cm are used to neutralize 
the rather slight influence of terrestrial mag- 
netism on the phenomena under examination. 
The center of the bar coincides with the center of 
the solenoid, and a coil 5.5 cm long of 4455 turns 
of No. 24 wire surrounds the central portion of 
the solenoid. This coil is connected with a low 
resistance galvanometer having a sensitivity of 
0.35 microvolt per millimeter, but it was used 
only as a ballistic in all our observations. 

Figure 1 shows the essential parts of the appa- 
ratus without the Helmholz coils. The com- 
mutator switches for reversing the longitudinal 
current J,, and the magnetizing current J,, are 
also omitted. 


PROCEDURE 


In order to obtain consistent results it was 
soon found necessary to establish a definite cycle 
of operations which, if rigidly followed gave 
galvanometer throws that could be repeated at 
any time with satisfactory regularity. The pro- 
cedure adopted was as follows: The iron is first 
carefully demagnetized. Then J,, is brought up to 
the value desired by cutting out resistance in its 
circuit until some point as a in the B—// curve, 
shown in Fig. 2, has been reached. If the circuit 
of the magnetizing current is then abruptly 
broken by opening Ky, the flux decreases from a 
to some point } on the B axis. Next the galva- 
nometer circuit is closed by means of K3, and any 
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desired current J; is sent through the bar by 
closing K». In those cases when the flux increased, 
the result is a throw of the galvanometer pro- 
portional to dc, or when the flux decreased, 
proportional to be. The next step is to open K; 
and close K, thus bringing the magnetism back 
to d or f. The cycle is then completed by closing 
K; and observing the galvanometer throw when 
I, is broken by opening K»2. This throw is pro- 
portional to da showing a decrease of B, or to af 
showing an increase. These values were generally 
nearly equal to those corresponding to cb and be 
but of course in the opposite sense. This cycle 
may be repeated as many times as desired with 
only minor changes in the readings. 

In our curves we have given changes in the 
residual flux, Ag, in terms of the galvanometer 
throws to which they are proportional. But in 
order to be able to interpret these throws in 
terms of maxwells we calculated the flux change 
for several deflections when there was no iron 
within the solenoid and found that one millimeter 
on the galvanometer scale corresponds to a change 
of 0.044 maxwell passing through the coil CC. 

Before describing a typical series of observa- 
tions we should point out two important con- 
siderations. In all cases in which the residual flux 
is increased by J;, the magnetizing current was 
broken abruptly.’ If instead J,, is reduced to a 
very small value by gradually introducing resist- 
ance in the circuit, the effect of J; is invariably to 
produce a decrease in the residual magnetism. In 
cases where J; after an abrupt break caused an 
increase, the gradual elimination of the mag- 
netizing current was followed by a moderate 
decrease of residual flux. But in cases where J, 
after an abrupt break, caused a moderate de- 
crease, a slow elimination of I,, was followed by a 
very large decrease of magnetization produced by 
I;. This must mean that an abrupt break brings 
the point 5 (Fig. 2) to a lower point on the B axis 
than when the current is decreased gradually. To 
test this peculiarity still further we tried an 
intermediate procedure of pushing the slider 
along the drum rheostat as fast as possible from 


§ We have found that in spite of the abrupt break of Im, 
the condition of the iron is surprisingly stable. Even when 
a week elapsed between subjecting the bar to a field of 
about forty oersteds and the closing of the J; circuit, the 
galvanometer indicated the same increase of flux as when 
there was no appreciable interval. 


positions of low to high resistance. This gave 
intermediate results. That is in cases where the 
magnetism was strongly increased by J), after an 
abrupt break of J,,, the less abrupt decrease of 
the magnetizing field still resulted in an increase 
of magnetism, though less than before. In cases 
where J;, after an abrupt break caused a decrease, 
the less abrupt elimination of J,, was followed by 
a decidedly larger decrease of flux, but not nearly 
as much as when the field was gradually 
eliminated. 

The other matter to be noted is that there was 
always a small inductive effect between the J; 
circuit and the coil CC. This resulted in a re- 
versible throw of the galvanometer on make and 
break of J; as has already been mentioned in con- 
nection with Villari’s experiments. This dis- 
turbing throw varied with the orientation of the 
wires in the galvanometer and longitudinal cur- 
rent circuits, and by careful adjustment could be 
made as small as two or three millimeters. As this 
appeared as a constant increase or decrease in the 
throws caused by changes of flux, we have cor- 
rected those readings accordingly. 


Fic. 2. Magnetic cycles. 


RESULTS 


Two typical curves obtained with increasing 
values of the magnetizing field (ignoring the 
demagnetization by the poles) from 0.0188 
oersted created by a current of one milliampere, 
up to 56.4 oersteds, (3 amperes), and with J; 
constant at some value as 0.5 ampere up to 4 
amperes (though usually 2) are shown in Fig. 3. 
Here the axis of ordinates represents A¢, while 
the axis of abscissas measures the magnetizing 
field in oersteds. One curve gives the behavior 
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when K¢ was closed and corresponds to da or fa 
in Fig. 2. It is labeled ‘‘make.”’ The other labeled 
“break”’ represents the behavior when Ky was 
opened with J,, flowing, and corresponds to da or 
fa. These curves are almost identical in shape 
though of course reversed. However near the 
point of reversal, a, there is considerable uncer- 
tainty due to double throws already mentioned, 
and its exact position could not be determined 
with precision. But in general, with J; at 2 


+|4¢ 


DEFLECTION IN MILLIMETERS 


Fic. 3. Typical make and break curves. 


amperes, it came when the field was in the 
neighborhood of 10 oersteds. 

The shape of these curves strongly suggests the 
existence of the two effects we have postulated. 
For instance the ‘‘make’’ curve, M, as shown in 
Fig. 4, may be decomposed into two other curves 
A and B indicating the increasing and decreasing 
flux caused by the A and B effects, respectively. 
This hypothesis was strongly substantiated by 
observing the effect of a current J, flowing in a 
wire placed in the axis of the hollow cylinder. 
Such an axial current sets up an even stronger 
circular flux than one flowing in the iron itself, as 
is found from the following calculation: If we 
assume the permeability of the iron to be the 
same when the same current is sent along the axis 
as when it flows through the iron, and if we 
suppose the current to be of uniform sectional 
density, then in the case of the axial current, a 
familiar calculation gives 


R 
2ip loge, (1) 


where R and ¢ are the outer and inner radii of the 
bored-out cylinder and ¢, is the total flux per 
unit length. In the case of the longitudinal cur- 
rent J;, the total flux per unit length is given by 


2r? R 
Xlog-— (2) 
r 


In our cylinders r=6 mm, and R=11 mm. 
Substituting these values and dividing (1) by 
(2) we obtain ¢./¢:=2.5 very nearly, Now when 
H is small (before the A effect has really de- 
veloped) the galvanometer throws indicate 
roughly twice as great a decrease of magnetiza- 
tion when the current is axial as that produced by 
a current through the iron. This is near enough 
the calculated ratio of 2.5 to suggest that the B 
effect tends only to demagnetize the iron. 

This demagnetizing effect can be discounted in 
two ways thus exposing, so to speak, the A effect 
more fully. One way was to send suitably adjusted 
currents of different values and in opposite 
directions through the iron and along its axis. 
Under these circumstances the B effect was very 
much reduced, though we were unable to neu- 
tralize it altogether; while the A effect gave us 
larger deflections than when a longitudinal cur- 
rent of the same value acted alone. Apparently 
the opposite direction of the currents did not 
alter the A effect but by reducing the B effect 
made it appear larger. 

The other and even more effective way of 
reducing the circular flux was to cut a slot in the 
cylinder as already explained. This greatly in- 
creased the reluctance of the bar to circular flux 
and greatly reduced the B effect. 

Our latest series of observations was made with 
the two similar hollow cylinders described above, 
one slotted, the other not. We made a series of 
readings for each of four cases beginning with a 
magnetizing field of less than 0.02 oersted, and 
ending with about 60 oersteds. The longitudinal 
and axial currents were both two amperes. These 
cases are: 

No. 1, current through iron, no slot, 
No. 2, current through axis, no slot, 


No. 3, current through iron, slotted bar, 
No. 4, current through axis, slotted bar. 


The curves obtained are very illuminating, and 
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are shown in Fig. 5 with case No. 2 plotted on a 
scale one-tenth that of the others with respect to 
the axis of abscissas which is laid off in galva- 
nometer deflections for convenience, although it 
measures Ad¢, the change of flux. 

If we contrast curves 2 and 4, it becomes 
evident that the slot has so greatly reduced the 
very large negative deflections of 2 that the A 
effect is again exhibited. A comparison of curves 
1 and 3 also illustrates the effect of the slot 
though in a lesser degree than when 2 and 4 are 
compared. The positive galvanometer throws 
begin earlier and reach higher values in curve 3 
than in the other three curves. It should further 
be noted that in this case the current flows 
through the iron of a slotted bar, so the circular 
flux is weaker than in the other cases. At any 
rate there is no evidence that the A effect is 
stronger in one case than in any other. 

We have not yet studied permalloy beyond 
obtaining a qualitative agreement with the 
Norway iron. That is, when sufficiently mag- 
netized, there is an increase in the residual flux 
caused by a current sent through a solid cylinder 
of the alloy. In the case of a cast iron cylinder 
bored out and slotted like our other cylinders, we 
find only a decrease in residual flux as far as we 
have been able to carry the magnetizing field. 
However it is significant that whereas an axial 
current causes a progressively larger decrease in 
the residual magnetism, as that is raised to 
higher and higher values, the effect of a current 
sent through the iron begins to diminish with 
large values of residual magnetism, and the 
curve, though still negative, slants upward 
toward the H axis which it seems to be ap- 
proaching. This suggests the possibility that 
with sufficiently large magnetizing currents 
(larger than our coil would carry) the curve might 
cross the H/ axis, indicating that the A effect 
exists even in cast iron. The nickel bar gave 
results similar to cast iron, but as yet we have not 
examined its behavior very thoroughly. 

A final experiment was a repetition of what we 
had roughly measured early in our work. This 
was a study of the effect of a mechanical jar on 
residual magnetism. The jar was produced by 
giving a light tap with a brass rod to the upper 
end of the bar which had previously been 
demagnetized both circularly and longitudinally, 


and then brought up to the desired magnetic 
intensity by first closing and then opening K,. 
The results were surprisingly similar to curve 1, 
showing first a decrease of residual flux, then, for 
higher intensities, showing a reversal associated 
with double throws, and finally quite large 
deflections indicating the usual increase of re- 
sidual magnetism. If however the magnetizing 
current was not broken, when H was small, the 
tap produced a rather large increase in mag- 
netization, while with large values of 17, beyond 
about twenty oersteds, the tap caused a decrease. 
Thus we found more or less complimentary curves 
like those of Fig. 3. 


DISCUSSION OF RESULTS 


In explaining our observations we have to 
account for a number of effects which are not 
easily harmonized. It is necessary to find one or 
more causes which will simultaneously account 
for the facts: (a) that a longitudinal current 
through an iron bar tends to decrease small 


a 


Fic. 4. Proposed decomposition of experimental curve. 


residual fluxes and increase large ones; (b) that a 
mechanical jar does the same thing; (c) that 
a current through the axis causes only a decrease 
in the residual magnetism, unless the bar is 
slotted, when results similar to (a) are obtained; 
(d) that these phenomena depend upon a sudden 
break in the magnetizing current, and that a 
more or less gradual decrease in that current 
tends more or less to result in decreases of the 
residual flux; (e) that in borderline cases there are 
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Fic. 5. Four experimental curves, with 2 plotted on a reduced scale of ordinates. 


two galvanometer throws, the first being invari- 
ably the one indicating an increasing flux. 

It might be argued that a circular field applied 
to a bar which has a residual longitudinal flux 
would produce a twist in the bar, and that such 
a twist combined with the circular flux would 
result in an increase in the residual flux provided 
the circular flux was small. This effect, based on 
Wiedeman’s experiments, was observed by 
Nagaoka and Honda‘ in 1902. It amounts to a 
‘Lenz law”’ reaction to the current which origi- 
nally created the longitudinal magnetization. 
But we do not feel that this explanation is 
satisfactory, because it does not explain the 
increased flux caused by a mechanical jar after 
the bar had been carefully deprived of any 
residual circular flux. Nor does it take into 
account the fact that slow demagnetization 
destroys the phenomenon. And further it does 


4H. Nagaoka and K.‘Honda, Phil. Mag. 4, 45-72 (1902). 


not account for the remarkable result of slotting 
the bar which, if anything, should result in a 
reduced twist and therefore less of the proposed 
“Lenz law”’ reaction, as a result of weakening the 
circular flux, whereas we found just the opposite 
to be true. 

It seems to us that a possible explanation is 
that on an abrupt break, the elementary dipoles 
in swinging away from the position of maximum 
magnetic intensity overshoot the mark, so to 
speak, by a kind of momentum, and are then 
brought back to a more stable position either by 
a trigger-like action of the mechanical jar, or by 
the abrupt setting up of a circular flux. This 
might account for our A effect. Our B effect is 
also brought about either by a jar or by abruptly 
setting up a circular field, but it invariably 
involves a decrease of magnetic intensity. This 
would be the expected result of a jar, and there- 
fore it seems necessary to suppose that the abrupt 
creation of a circular field has much the same 


a 
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effect, in that it tends to throw the elementary 
dipoles out of a longitudinal alignment. Thus the 
same two causes (jar and circular flux) seem to 
produce two effects, one tending to bring about a 
more stable alignment of elementary magnets 
that have overshot the mark, and the other 
tending to turn them out of that alignment. The 
fact that in the case of double throws (observed 
both with mechanical jar and abruptly created 


circular field) the trigger effect comes first shows 
that the decrease in longitudinal magnetization 
takes more time to be produced than the increase 
caused by a jar or by circular flux after the 
magnetizing field has been abruptly broken. This 
seems reasonable, because upsetting an unstable 
condition should take less time than a more or 
less permanent rearrangement of the dipoles 
associated with a decrease of magnetization. 
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In the precise molecular beam experiments of Rabi and his collaborators, it is necessary to 
know the value at the nucleus of the magnetic field produced by the diamagnetism of the atomic 
electrons. This has been calculated on the basis of the Fermi- Thomas model, and checked for a 
number of atoms by use of the available Hartree calculations. The statistical treatment gives 
for ratio of induced to external field 0.319 X 10-*Z*/8, while the numerical coefficient on the basis 
of the Hartree model is lower by 19 percent at Z = 19 and by 12 percent at Z = 80. The effect is 
equivalent to a reduction of the nuclear g value by a factor of (1—0.319X 10-#Z*/8), and in this 
form, the correction may be applied in the calculation of hyperfine structure of heavy atoms. 
The influence of the diamagnetic fields on an orbital electron has also been considered, and it is 
shown that it is equivalent to a reduction in the g value of an outer s electron by an amount of 
just the same order of magnitude as the relativistic correction calculated by Margenau. 


N discussing the magnetic properties of solids, 

it is important to know the value of the 
magnetic field produced at one atom due to the 
action of all the other atoms. One uses for this 
some modification of the Lorentz formula.' We 
shall be concerned in this note with the field 
produced within an atom by its owndiamagnetic 
moment when an external field is present. This 
problem has only arisen because of the very 
precise molecular beam measurements of Rabi 
and his collaborators.** 


* Publication assisted by the Ernest Kempton Adams 
Fund for Physical Research of Columbia University. 

1H. Casimir, Magnetism and Very Low Temperatures 
(Cambridge, 1940), p. 56. 

2 P. Kusch, S. Millman, I. I. Rabi, Phys. Rev. 55, 1176 
(1939); P. Kusch, S. Millman, Phys. Rev. 56, 527 (1939) - 
R. H. Hay, Phys. Rev. 60, 75 (1941). 

3 P. Kusch, S. Millman, I. I. Rabi, Phys. Rev. 57, 765 
(1940); S. Millman, P. Kusch, Phys. Rev. 58, 438 (1940); 
60, 91 (1941). 


An external field H (taken along the z axis) 
may be described by the vector potential 


A=}(Hxr] (1) 


and this induces a diamagnetic current density in 
the atomic electrons of* 


S=eAp(r)/mc, (2) 


where —e is the charge on the electron and p(r) is 
the charge density at r. The induced field is then 
given by 


A'(1) = (e/2me?) f dr’ (3) 


When p(r) is spherically symmetrical, this may 


4 This is most easily derived by considering the induced 
current in a conducting ring of radius r in the xy plane. 
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be written as 


+f 

r’>r (4) 

whence the induced magnetic field along the z 
direction is 


= (2? — 4r?) dr’r’?p(r’) 


r'<r 


+ (ell /3mc*) dr'p(r’)/r’. (5) 
r’>r 
The induced field at the center of the atom, i.e., 
that acting on the nuclear spin is just 


(0), = (eH /3mc?) dr’ p(r’)/r’ 


=(el1/3mc*)v(0) (6) 


and depends only on the electrostatic potential 
v(0) produced at the nucleus by the electrons. 
One has 0(0) = —Ze(1/r), and using an effective 
screening radius r~do/Z!, one finds 


H'(0).= = —ya?Z*"H, (7) 


where y is a number of order unity, do is the Bohr 
radius, and a is the fine structure constant. The 
value of y may be estimated by use of the Fermi- 
Thomas atom model. According to this5 


V(r) = (Ze/bx) o(x) (8) 


is the total electrostatic potential at a distance 
from the nucleus r = bx, where b =0.885a0/Z!, and 


(x) =1—1.588x%+4/3x!+--- (9). 


for small x. The electronic contribution to the 
potential at the nucleus is then 


v(0) = —1.588Ze/b, 


whence y has the value (1.588/0.885) /3 =0.598 


and 
H’'(0),= (10) 


5 E. Condon and G. Shortley, Theory of Atomic Spectra 


(Cambridge, 1935), p. 336; L. Brillouin, L’Atome de 


Thomas-Fermi (Hermann, Paris, 1934). 


The correction in magnetic field acting on the 
nucleus must also be applied to the field produced 
by an orbital electron when one calculates 
hyperfine structure. It is easy to show that this 
case may be combined with that of an external 
field by keeping the applied field constant, and 
changing the nuclear g factor to g(1—0.319 
xX 10-4Z*/*), except for some negligible errors due 
to the field produced by the outermost electrons. 

These simple formulae can be improved in any 
particular case by use of the Hartree model. For 
a number of atoms, values of v(0) have been given 
explicitly by Hartree and others. Thus for® 
wK, 20Ca, e6Fe, 37Rb, ssCs, soHg, one 
finds that the number 0.319 in Eq. (10) should be 
replaced by the numbers 0.259, 0.259, 0.263, 
0.268, 0.270, 0.274, 0.277, 0.280, respectively. 
Thus it is seen that for these elements, the Fermi- 
Thomas Z*/8 law is checked very well, although 
the numerical coefficient is too high by an 
amount varying smoothly from 19 to 12 percent. 

In the radiofrequency experiments with atoms,’ 
one has beside the nuclear spin, an orbital electron 
usually in a *S, state. We shall now discuss the 
effect of the induced field on this outer electron. 
It is seen from Eq. (5) that the induced field H’, 
at any distance r may be divided into two parts. 
One is due to electrons inside the sphere of radius 
r, the other due to electrons outside of the sphere. 
The first contribution varies like the second 
Legendre polynomial P»2 (cos@) with angle, and 
averages to zero in its effect on an s electron. The 
second contribution does not vary with angle, 
but it is small for an outer electron which does not 
penetrate much into the core. We can estimate 
the effective induced field for such an electron as 
follows. We must average 


H" dr'p(r’)/r’ (11) 


r’>r 


over the orbit of the outer electron. With Poisson’s 
equation 


4rp(r) (12) 
r 
dr? dr 


® The references used are given in F. Seitz, Modern 
Theory of Solids (McGraw-Hill, 1940), p. 251, except for 
the case of Fe which was obtained from M. Manning and 
L. Goldberg, Phys. Rev. 53, 662 (1938). 
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the integral, which is just the potential due to the 
electrons of the core beyond 7, can be expressed as 


f =V(r)+r(dV/ar), (13) 


where V(r) is the potential acting on the orbital 
electron. If we write V(r) =eZ(r)/r, this becomes 
edZ(r)/dr, so that the effective magnetic field is 


(e?/3me?) (dZ(r) /dr) w. (14) 


By applying the virial theorem to a potential of 
the above form, one finds 


W—T, (15) 


where W is the binding energy of the orbital 
electron and T is its average kinetic energy. The 
two quantities are of course equal for a pure 
coulomb field, but the average kinetic energy is 
larger than the binding energy for an atomic field 
because of the penetration of the orbital electron 
into regions of the core where the effective Z is 
larger than unity. A very simple, although rough 
way to estimate the difference T—W is the 
following: For orbits with small penetration, one 
can expand Z(r) as a power series in inverse 
powers of r. For large r only, one could write 


Z(r) =1+¢(ao/r), where ¢ is a constant. 


Then 
T — (16) 


Now the binding energy of the electron is equal to 
the hydrogenic value Wo=e?/(2n%ao) increased by 


Cae?(1/r?) ay. 


Hence T—W=—W-— Wo, and we can express our 


results simply in terms of the excess binding 
energy of the electron over the hydrogenic value, 
so that 


w= — (H/3mc*)(W— W). (17) 


This induced field is always of the order 10° of 
the applied field, and hence does not affect the 
discussion of the radiofrequency experiments in 
their present state of precision, although this can 
probably be increased. 

It will be noticed that the atomic electron 
behaves as if its g value were reduced from the 
value 2 to the value 


g=2(1—(T—W)/3me). (18) 


This is of the same order, but distinct from the 
relativistic change in the g factor discussed by 
Margenau,’ who found 


g=2(1—3(T/mc*) (19) 


His effect arises because of the change in the 
magnetic moment of an electron due to its motion 
in a fixed potential field, while we have con- 
sidered the fact that the potential field itself is 
changed by the presence of an external magnetic 
field, at least in the present case where this field 
is produced by electronic charges. Both cor- 
rections must be considered in order to obtain the 
effective g value of an s electron. 

These calculations were suggested to me by 
members of the molecular beam laboratory of 
Columbia University. I am indebted to them for 
helpful discussions. 


7H. Margenau, Phys. Rev. 57, 383 (1940). It has been 
shown by M. Phillips, Phys. Rev. 60, 100 (1941) that 
interconfiguration interactions do not affect the g value 
of an alkali atom in the ground state. 
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As is well known, Peierls and Debye found that, if 7>>@ (T is the temperature, @—the Debye 
temperature), the thermal conductivity « of dielectrics varies inversely to the temperature. 
However, this turns out to be true only for dielectrics obeying the definite law of sound velocity 
dispersion and the definite law of the dependence of sound velocity on the direction. For all other 
dielectrics at T>>@, x is inversely proportional to 7. This results from investigation of the free 


path of long elastic waves. 


EIERLS'! and Debye? have found the depend- 
ence of thermal conductivity « of dielectrics 
on the temperature for the cases in which T is 
very much greater than @ the Debye temperature. 
In order to clarify the problem we shall consider 
the situation which arises in the body with no 
dispersion of sound and with isotropic sound 
velocity. 

As is known,’ the absorption coefficient of the 
transversal sound (i.e., of transversal phonons) 
varies as w, w being a circular frequency. It is 
easily found that at 7>@ the phonon free path / 
is of the form: 


l~a(Ms?/KT)(K0/hw), (T>8). (1) 


Here M is the mass of the elementary cell, K 
the Boltzmann constant, s the velocity of sound, 
a the lattice constant. Absorption coefficient of 
longitudinal sound depends on the existence of 
the possibility of the absorption of a long wave- 
length longitudinal phonon by short wave-length 
phonons. From the conservation laws it follows, 


that 
f+F=F’; (2) 


Here f stands for the wave vector of a long 
wave-length longitudinal phonon, F for wave 
vector of a Debye short wave-length phonon. 
We get: 


Ow 
w(f) =f—-; Siong f =fVr. (3) 
oF 


1R. Peierls, Ann. d. Physik 3, 1055 (1929). 

2 P. Debye, Vortrage uber die Kinetische Theorie etc., 
s. 43 (Teubner, 1914). 

3L. Landau and G. Rumer, Sowjetunion Physik. Zeits. 


11, 18 (1937). 


Siong IS a longitudinal sound velocity, Il’; is the 
group velocity of a Debye phonon. 

If, as for the body under consideration, 
Siong> Vp, the process (2) is impossible. For 
bodies having such a dispersion law and an 
anisotropic sound velocity, the absorbtion of a 
long wave-length longitudinal phonon is only 
possible if it splits into two small transversal 
phonons or into one small transversal and one 
small longitudinal phonon. The absorption coeffi- 
cient turns out to be proportional to w‘ (in 
distinction to w*, as was obtained by Slonimsky‘). 
The free path /'!! may be written, if the order of 
the value is only concerned, as follows: 


(4) 


Moreover, the absorption of longitudinal 
phonons due to collisions involving four phonons 
must be taken into account as well. The corre- 
sponding free path is of the form: 


I'V = (5) 


Equations (4) and (5) give equal quantities at 
the frequency wo, at which 


hwo=K0(KT/Ms*)', (T>). (6) 


It is evident that in these conditions the heat is 
transferred by small longitudinal phonons be- 
cause they have enormous free paths. 

The thermal conductivity « equals: 


f Ce dO=sinddddyg. (7) 


~ 4Slonimsky, J. Exper. and Theor. Phys. (Russian) 7, 
1457 (1937). 
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Here C, is the specific heat of a given phonon. 
Because small frequencies alone are of any 
importance, C, may be put equal to K. Intro- 
ducing Eq. (4) into (7) we get after integration 
up to wo, this limit resulting from Eq. (6): 
—(— KO 


hs hwy 


Ms? 1 
KT hs 
Ms?\} 1 /K@a\? 
KT / hs 
Let us designate 1/a* by n, which is the quantity 


of the same order as the total number of atoms 
(or molecules) per cc, and use the relation: 


Ké~ (hs /a)(2n*)}. 


«= aKs— 


KT 


Then 
(T>@6). 


Thus for the bodies, for which Siax2 Vr, Smax 
being the velocity of the oscillations having the 
greatest velocity, the thermal conductivity at 
high temperatures varies as 7~!, in distinction to 
the results of Peierls' and Debye.* 

If Smax< Vy, the difference with respect to the 
absorption between ‘‘longitudinal’’ and ‘‘trans- 
versal’”’ sound vanishes and the result obtained 
by Peierls! and Debye? remains to be true. 

At low temperatures (7<6), according to 
detailed calculations, only the admixtures, the 
reflection from boundaries of small crystals and 
the collisions of triples of phonons are essential 
for the evaluation of «x. In distinction to the case 
treated by Peierls,’ here the finite value of 


5 R. Peierls, Ann. d. Physik 3, 1097 (1929). 


thermal conductivity cannot be obtained if we 
consider only the scattering of phonon on the 
atoms of admixture and the collision of triples of 
phonons, because in both effects the free path 
varies as w~*. The finite value of thermal con- 
ductivity can be obtained only if the reflection 
from small crystal boundaries is taken into 
account. The resulting expression for « turns out 
to be inverse to L?, L being the linear dimension 
of the crystal. Therefore, at low temperatures the 
thermal conductivity, as found experimentally, 
depends (although not strong) on the shape of 
the crystal. This holds even if the phonon free 
path for the frequency of the order KT /h is 
much less than the crystal dimensions. If the 
temperature is decreased and becomes less than 
the Debye temperature, thermal conductivity 
at first increases inversely to 7?, but finally it does 
not depend on temperature at low temperatures 
for which the waves with w= KT7/h have the free 
path of the same order, as the crystal dimensions. 
Further decrease of temperature as is known, 
leads to variation of x proportional to 7*.*7 

The temperature range in which x does not 
change at all is especially wide for the diamond 
because of its high value of Debye temperature. 
This result explains the experimental data con- 
cerning the thermal conductivity of diamond,* 
according to which data «x reveals no noticeable 
change in the temperature range from 24 to 343°. 

Detailed calculations have been published in 
Journal of Physics [4, 259 (1941) ]. In the same 
publication the calculations involving the region 
of low temperatures are to appear shortly. 


6 W. J. de Haas and T. Biermasz, Physica 5, 619 (1938). 
7H. B. G. Casimir, Physica 5, 495 (1938). 

* Eucken, Physik. Zeits. 12, 1055 (1911). 

® W. J. de Haas and T. Biermasz, Physica 5, 47 (1938). 
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The distribution of characteristic vibrations of the KCI crystal has been calculated following 
Born, v. Karman and Blackman under the assumption of both electrostatic and repulsive forces 
between the ions, the latter calculated for the 18 nearest neighbors from the elastic constants 
and the residual ray frequency. Application of this distribution has been made to the evaluation 


of the specific heat at low temperatures. 


INTRODUCTION 


give detailed explanation of many proper- 

ties of crystals, e.g., temperature de- 
pendence of the specific heat and the diffuse 
scattering of x-rays, it is desirable to know the 
equation of motion and the distribution of 
characteristic vibrations. These can be calculated 
from the lattice configuration and the forces 
between the atoms. The latter have often been 
calculated by means of the elastic constants, but 
this method allows the determination of only 
a few force constants and is therefore only useful 
as long as short range forces are considered. This 
method is hence applicable to atomic crystals, 
such as the tungsten crystal, which has been 
treated in detail by Fine.' In the case of ion 
crystals it is necessary to notice that the electro- 
static forces act also between ions of greater 
separation and have to be considered in a dif- 
ferent way. The method of considering the elec- 
trostatic forces separately has been used by 
Lyddane and Herzfeld? for calculation on NaCl, 
more recently by Kellermann® for the same 
crystal. Ljungquist* and Foldy® use this method 
on the KCI crystal. Lyddane, Herzfeld and Foldy 
use theoretical statements for the short range 
forces. Kellermann calculated the force constants 
between 6, Ljungquist between 18 nearest 
neighbors using the experimental values of the 
elastic constants. 


* Now at the University of Chicago. 

1P. C. Fine, Phys. Rev. 56, 355 (1939). 

*R. H. Lyddane and K. F. Herzfeld, Phys. Rev. 54, 
846 (1938). 

3E. W. Kellermann, Phil. Trans. Roy. Soc. 238, 513 
(1940). 
*S. Ljungquist, Uppsalg. To be published. 
5L. L. Foldy, Phys. Rev. 60, 64 (1941). 


THE FREQUENCY SPECTRUMt 


To calculate the distribution function for the 
simplest cubic crystal, the KCl crystal, the 
masses of the ions may be assumed to be equal 
and, under the assumption of equal forces 
(except for the difference in sign of the charge) 
between the ions, the crystal may be considered 
as a simple cubic lattice. In order to find the 
number of characteristic vibrations in given 
frequency intervals it is necessary to consider 
the equation of motion of the ions, which may be 
taken to be oscillating as a whole. One obtains 
these equations from the general equation of 
motion in the following way: the gradient of the 
potential energy of the ions, expanded in a power 
series of the displacement of the particles, is 
substituted for the force and the motion is taken 
to be oscillatory. If one assumes a plane wave 
solution, the displacement of the ions, whose 
rest position is characterized by the vector l 
with integral components, may be described as 
follows: 


u, = exp[i(gl) (1) 


U indicates the amplitude of the oscillations with 
frequency w=2z», ¢ is a vector in the direction 
of wave propagation of the magnitude 27a/n, a 
is the distance between nearest lattice points and 
the wave-length. The vector ¢, whose com- 
ponents g; are usually called phase components, 
characterizes thus both wave-length and direc- 
tion of the waves. To be free from the boundary 
conditions of the crystal it is convenient, follow- 
ing Born, to consider a “‘cyclic lattice;’’ that is 
an infinite lattice which repeats the finite lattice 
inclusive of its motion periodically in each direc- 


t For the formulae used see reference 11. 
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tion. If the considered finite lattice has n ele- 
mentary cells in each direction, each cell con- 
taining one particle, the condition of periodicity 
of motion makes ng; an integral multiple of 27. 
Since the wave function has physical significance 
only at the lattice points the g; may be restricted 
to The values ¢; (gi=27p;/n) build 
a cubic lattice with the lattice constant 2x/n 
and the edge length 27; this is the so-called phase 
cube. The number of lattice points in a volume 
element (dg), or the number of characteristic 
vibrations which is proportional to this number 
is given by 


dn =3(n/2m)3(d¢g) (2) 


since to each point of the phase cube three 
characteristic vibrations correspond. The inte- 
gration of expression (2) gives 3n*, the number of 
degrees of freedom. In order to find from the 
number of vibrations in the volume element of 
the phase cube the distribution of lattice vibra- 
tions along the frequency scale the relation 
between frequency and phase components must 
be found. One obtains this relation from the 


2 2 


(Nuttin) 


condition that the equations of motion give the 
ratios between the amplitude components, so 
that the determinant of the equations of motion 
vanishes. 

With the assumed oscillation (1) the equations 
of motion turn out to be: 


(Aut+mw*) +A 2U2+A13U3= 
U2+A23U3=0, 
mw?) U3=0. 


The coefficients Ax(¢) are given from the 
derivatives of the potential energy of the ions in 
direction of the i and k coordinates by 


Au(¢)= 2’ (4) 
These coefficients can be divided into two parts, 


one due to the electrostatic forces (A in) and the 


other to the repulsive forces (A ak, 

The method given by Ewald® allows the cal- 
culation of the electrostatic parts. Following 
Waller’* one obtains: 


m=1 


m 


2 


In these equations }°y is a sum over all triplets 
of odd numbers, and >>’; a sum over all triplets 
of integral numbers except (0, 0, 0), G(x) is the 
error function 


2 z 
G(x) = f exp(—a?)da 


and G’’(x) its second derivative. The length of 
the vector J is / and yw;=¢;/r. By assuming a 
proper value for the constant 7 it is possible to 
find sufficiently convergent series of both types. 
The calculations have been carried out with n=1 
which give the higher number of terms in the 
more convenient second sum. The coefficients 
have been checked by means of the Laplace 
equation or by using different values for n. To 


check the method the electrostatic parts of the 
elastic constants have been calculated and com- 
pared with the sums given by Born. 

The coefficients A’. have been calculated by 
following the method of Born, v. Karman,’ and 
Blackman.’ For this purpose the force constants 
P and Q have been used which Ljungquist* ob- 
tained from the experimental values of the 


6 P. Ewald, Ann. d. Physik 64, 253 (1921). 

7T am much indebted to Professor I. Waller for giving 
me these unpublished formulae. 

* At the same time Kellermann calculated in England 
the coefficients for the NaCl crystal. The coefficients 
given above may be obtained from his formulae as the 
special case for equal ion masses. 

8M. Born, Th. v. K4rman, Physik. Zeits. 13, 297 
(1912); 14, 15, 65 (1913). 

9M. Blackman, Proc. Roy. Soc. A148, 365, 384 (1935); 
A159, 416 (1937). 
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elastic constants (¢;; and ,) given by Forsterling,!® 

from the residual ray frequency (Ao=70.7 X 10-4 

cm) and from the equilibrium condition after 

splitting these values in electrostatic parts! and 
the parts due to the repulsive forces. 

Ajs($) = Pi2)(1 —cosg;) 

2P12(2 — cosg; 
P11)(2 —cosy; cosy, 
—COS¢; cosy) 
—4P11(1—cosy, cosg;), (6) 

Aix($) = —4Qna? sing; sing. i+k+1+i. 

The numerical values of the constants are: 
P\,= —0.0469 x 10°; —1.992 
Qiu= —0.0315 10%; Qie=2.5218X 10". 

TABLE I. Values of electrostatic parts of coefficients Aix. 


— 282 — 282 +563 0 0 0 
— 978 -—978 +1955 0 0 0 
—16.97 -16.97 +3393 0 0 0 
—20.01 —20.01 +4003 0 0 0 
612 +306 +306 0 0 — 269 
-1434 335 +1768 0 0 
22.97 - 9.98 +3295 0 0 —299 
—26.66 -12.83 +3048 0 0 0 
—25.95 +1298 +1298 0 0 — 6.0 
—38.51 +803 +3049 0 0 4.32 
+588 +3814 0 0 0 
—55.83 +2792 +2792 0 0 343 
—63.59 +2681. +3679 0 0 0 
-7244 +3622 +3622 0 0 0 


— 335 — 335 — 3.35 
— 7.70 — 7.70 +1540 — 5.10 — 5.08 — 5.08 
—15.89 —15.89 +3178 — 7.17 — 387 — 3.87 


—1945 —-1945 +3890 — 8.14 0 0 
—19.07 +954 + 9.54 -— 811 811 — 7.88 
—31.86 + 3.06 +2880 -—1190 — 649 — 6.14 
—37.66 + 0.16 +37.50 —13.76 0 0 
—50.95 +2548 +2548 973 9.73 — 4.94 
—60.07 +23.98 +3609 —11.84 0 0 
—71.03 +35.52 +35.52 0 0 0 
0 0 0 —13.52 -—13.52 —13.52 
—11.97 —11.97 +23.96 —21.54 —11.53 —11.53 
—17.98 —17.98 +35.96 —25.88 0 0 
—34.49 +17.25 +17.25 —20.64 -—20.64 —10.89 
—47.87 +1346 +3442 —26.42 0 0 
—67.58 +33.79 +33.79 0 0 0 
0 0 0 —25.82 -—25.82 —25.82 
—16.36 -16.36 +32.74 —41.23 0 0 
—64.09 +32.05 +32.05 0 0 0 
—62.66 +31.33 +31.33 0 0 0 (x,1,1)* 
1 1 15.67 —15.67 +3133 —46.99 0 0 (x,x,1) 
0 0 0 —31.33 —31.33 —31.33(x,x,x) 


* Indicates the direction in which this point is approached. 


” K. Forsterling, Zeits. f. Physik 2, 172 (1920). 

"M. Born and M. Géppert-Mayer, Handbuch der 
Physik Vol. 24-2 (1934); M. Born, Enzyklopédie d. Math. 
Wiss. Vol. 5-3 (1923). 


The coefficients A, are uniquely determined 
at all points of the phase cube except at the 
point (x, x, 7). Terms of the magnitude of the 
reciprocal vacuum wave-length have been ne- 
glected in comparison to terms of at least the 
magnitude of the reciprocal lattice constant, 
which means we have neglected the retardation. 
At the point (z, z, 7) these latter terms vanish 


and the coefficients A depend on the direction 
from which this point is approached. At this 
point only the retardation cannot be disregarded 
but the distribution of vibrations is not changed 
by the value at this cornerpoint. Even at this 
point the solution of the secular equation gives 
a distinct set of frequencies. ft 

Because the coefficients depend only on the 
magnitude of the phase components it is sufficient 
to carry out the calculations for the positive 
octant of the phase cube. As a permutation of 
the g; gives the same frequency triplets it is suf- 
ficient to calculate the frequencies at 35 points 
in order to know them at all points differing by 
tm in the phase coordinates. Because of the 
symmetry of the coefficients all but one of the 
secular equations split up into equations of lower 
degree. 

The numerical values of the coefficients and 
the frequencies are collected in Tables I-III. 
Table I gives the electrostatic parts of the coef- 
ficients calculated with the factor e?/a*=7.479 


X10*. Table II contains the coefficients Ax and 
Table III gives the frequencies. 

Born, v. Karman, and Blackman obtain with 
coefficients analogous to (6) at the point (7, z, 7) 
of the phase cube three equal frequencies cor- 
responding to the residual ray frequency. Con- 
sidering the long range Coulomb forces one 
obtains two vibrations with residual ray fre- 
quency and one of almost double this frequency 
(Ao = 39.9 10-4 cm). 

The three frequencies at each point of the 
phase cube have been arranged in three branches 
of vibrations, one containing the highest, one the 
medium and one the lowest frequencies at each 
point. It was possible in this way to combine all 
frequencies by three smooth surfaces of the 


t For more detailed discussion of the effect of retardation 
see reference 3, and J. H. Thompson, Proc. Roy. Soc. 
A149, 487 (1935). 


i 
t 
P 
Oo 
E 
it 
d 
fe) 
tl 
in 
of 
fr 
cc 
tic 


LATTICE VIBRATIONS OF KCI CRYSTAL 825 


TABLE II. Values of repulsive parts of coefficients Aix. 


TABLE III. Values of the frequencies. 


An he Aw in 
x10-3 x13 XK 103 


“2 

+166 +164 12.53 0 0 0 
0 0 4 +560 +560 42.76 0 0 0 

0 0 2 +956 +956 —73.00 0 0 0 
0 +1120 41120 85.53 0 0 0 
& +326 -1101 0 0 +0.62 
+718 -—734 —41.54 0 0 +0.88 
 & +1110 -—368 —72.08 0 0 +0.62 
#1 +1273 —216 84.72 0 0 0 
+1101 38.59 0 0 +1.24 
0 & +1484 35.64 69.84 0 0 +0.88 
0 #1 +1642  -82.78 0 0 0 
& +1857 -67.60 -—67.60 0 0 +0.62 
2 #1 +2011 66.68  —80.84 0 0 0 
1 +2164 80.04  —80.04 0 0 0 

¢ —951 -—951 -—951 40.62 +4062 
—226 -—226 -T7116 +4062 +062 +0.62 
1 —07%6 +0.62 0 0 

%—235 -3743 +4088 4088 +1.24 
+4118 -6897 +4088 +4062 +0.88 
1 +264 +0.88 0 0 
4 2 & +462 -66.77 -66.77 +4062 +062 +0.62 
4 2 1 +603 -80.13 +0.62 0 0 

+745 -7935 79.35 0 0 0 

-31.79 66.86 41.24 4088  +0.88 
1 -80.23 +41.24 0 0 

} 2 2 64.76 -64.76 +4088 +4088  +0.62 
2 1 63.89 -7842 +0.88 0 0 

4 1 1 -77.67 —T77.67 0 0 0 

2 62.75 -62.75 +0.62 +0.62 
1 -76.70 +0.62 0 0 

1 1 -7599 —75.99 0 0 0 

1 1 1 —75.30 —75.30 0 0 0 


¢i¢g2¢3w space. At points where two frequencies 
are equal the frequency branches change from 
one smooth surface to the other. The branch of 
highest frequency contains mainly longitudinal, 
the two other mainly transverse vibrations. 
This combination of vibrations enables one to 
interpolate frequency values at other points of 
the phase cube. The frequencies have been inter- 
polated at points in the center between each pair 
of points of distance jr. To find, following 
Blackman, the number of vibrations in frequency 
intervals the whole frequency range has been 
divided into 24 intervals and three sets of curves 
of constant frequency, corresponding to the 
three frequency branches, have been determined 
in the 9 planes g3=0, jw, By means of 
graphical integration the areas between curves 
of constant frequency have been evaluated and 
from these the volumes between surfaces of 
constant frequency for all three branches have 
been determined by repeated graphical integra- 
tion. These volumes are proportional to the 


0 0 } 1.06 0.44 0.44 
0 0 3 1,94 0.82 0.82 
0 0 i 2.51 1.09 1,09 
0 0 1 2.71 1.19 1,19 
0 3 i 1.27 0.97 0.68 
0 } 3 1,99 1.27 1.08 
2.52 1.47 1.39 
0 } 1 2.70 1.56 1.50 
0 3 3 2.21 1.84 1.55 
0 3 ; 2.55 2.07 1.96 
0 3 1 2.69 2.15 2.11 
0 3 2? 2.62 2.45 2.44 
0 } 1 2.67 2.65 2.54 
0 1 1 2.86 2.66 2.66 
1.55 1.05 1.05 
} 2.14 1.54 1.21 
3 2.56 1.94 1.36 
i 3 1 2.70 2.12 1.43 
3 2.55 1.85 1.60 
3 2.82 2.32 1.80 
3 3 1 2.76 2.68 1.85 
3.08 2.44 2.31 
i i 1 3.13 2.67 2.38 
} 1 1 3.20 2.66 2.66 
3.10 1.90 1.90 
5 3 i 3.41 2.38 1.95 
3 3 1 3.44 2.67 1.97 
3 3 3 3.79 2.45 2.29 
3 3 1 3.85 2.67 2.36 
3 1 1 3.90 2.66 2.66 
i i i 4.27 2.46 2.46 
1 4.38 2.66 2.47 
3 1 1 4.50 2.66 2.66 
1 1 1 4.72 2.67 2.67 


number of characteristic vibrations. The smooth 
curve in each frequency interval containing the 
number of vibrations which has been calculated 
in this way was assumed to be the distribution 
function. 

The two branches of lower frequency show 
maxima of vibration density near the residual 
ray frequency. The lowest branch has this fre- 
quency as maximum frequency, the medium 
branch has a small number of vibrations with 
frequencies higher than the residual ray fre- 
quency. Both have second density maxima at 
lower frequencies, thus causing a second maxi- 
mum of the distribution of all three branches 
together. The branch of highest frequency has 
its maximum also near the residual ray frequency 
but has a high number of vibrations between this 
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and almost the double frequency. Figure 1 gives 
the sum of the distribution functions of all three 


branches. 
Speciric HEAT 


The specific heat of a system of N oscillators 
of frequency v is given by 


Cah (hv/kT)* exp(hv/kT) 
(exp[hv/kT]—1)? dv 


In the case of a quadratic distribution function 
dn/dv=const. Xv, with a maximum frequency, 
which is determined by the total number of 
degrees of freedom, a characteristic temperature 
6 can be defined, which is related to the maximum 
frequency by 0=Aimax/k. 

For the lower frequencies the distribution 
function of the KCI crystal is in agreement with 
the Debye quadratic distribution of the con- 
tinuum with the characteristic temperature 
6= 230°. This @ value is in agreement with the 
value obtained from specific heat measurements 
because this parabola is a good approximation for 
the distribution function even at higher fre- 
quencies because it continues to frequencies 
higher than the residual ray frequency. However, 
in temperature intervals where the first maximum 
of the distribution function gives the main 
contribution one may expect a minimum of the 
corresponding 6 values of the continuum which 
would give the same specific heat at this tem- 
perature. Blackman’s qualitative considerations 
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Fic. 1. Density distribution of the lattice vibrations of 
the KCI crystal. 
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Fic. 2. Specific heat. @—Keesom and Clark. +— 
Nernst and Lindemann. g—Southard and Nelson. 
—calculated 


and Kellermann’s” calculation of the NaCl crys- 
tal show the same feature. It is hardly possible to 
give a quantitative discussion of this feature 
without calculating the distribution with elastic 
constants at lower temperatures. But under the 
assumption that the distribution function does 
not change with temperature, although the 
elastic constants do," the dependence of the 
characteristic temperature on the absolute tem- 
perature can be calculated. This leads to the 
function given in Fig. 2. In this diagram the 
experimental values given by different authors 
are also plotted. Keesom and Clark" find a 
probable increase of @ in a range between 2°K 
and 4°K but it is not possible to make any de- 
cision about the reality of such an increase from 
the present calculations because the frequency 
intervals used are too large and the number of 
calculated frequencies too small. Above these 
temperatures they find a decrease to about 17°k, 
but the given values may be interpreted as a 
slight increase for @ at the end of this interval. 
The measurements are extended to higher tem- 
peratures by Nernst and Lindemann,’ and 
Southard and Nelson.'® The combined values 
seem to indicate a minimum near 15°K. 

I am much indebted to Professor I. Waller for 
suggesting this subject and for much advice 
during the work and to Internationella Student- 
hjalpen (I.S.S.) Uppsala for making it possible for 
me to do this work. 


122 E, W. Kellermann, Proc. Roy. Soc. A178, 17 (1941). 

13M. A. Durand, Ph 7 Rev. 50, 449 (1936). 

4 W. H. Keesom and W. Clark, Physica 2, 698 (1935). 

6 W, Nernst and T. Lindemann, Zeits. f. elec. Chemie 
17, 817 (1911). 

16 J. C. Southard and R. A. Nelson, J. Am. Chem. Soc, 
55, 4865 (1933). 


\ 
I] 
r 
n 
fe 
P 
n 
al 
of 
vi 
pe 
cc 
th 
fir 
en 
pa 
ef 
to 
by 
eff 
stc 
th 
= 
Re 
anc 
K. 
(19 
D. 
2 


DECEMBER 1, 1941 


PHYSICAL REVIEW 


VOLUME 60 


Range-Velocity Relation for Fission Fragments in Helium 


J. K. 
Institute of Theoretical Physics, University of Copenhagen Copenhagen, Denmark 
(Received September 29, 1941) 


The range-velocity relation of fission fragment tracks in helium has been studied by means of 
individual close collisions as well as by statistical considerations of the branch frequency. The 
comparison with previous investigations of fragment tracks in argon indicates that the shape of 
the range-velocity curves for the two gases differs appreciably in the last part of the range just 
as was to be expected from theoretical considerations concerning the relative contributions to 
the stopping power by electron encounters and nuclear collisions in light and heavy gases. 


S reported in previous papers,' cloud-cham- 
ber tracks of fission fragments in argon as 
well as in helium have been obtained in this 
Institute. While the tracks in argon offered 
sufficient material for a detailed study of the 
range-velocity relation, the tracks in helium 
were so far mainly used for range measurements. 
The ratio between the range of the fission frag- 
ments and the range of a-rays was in helium 
found to be about 20 percent larger than the 
corresponding ratio in argon. As pointed out by 
Professor Bohr, this difference is just in agree- 
ment with the different shape of the range- 
velocity curves in light and in heavy gases to be 
expected from theoretical considerations.” 

In fact, the range-velocity curve consists of an 
almost linear slope extending over the larger part 
of the range, and an end part, where the rate of 
velocity loss to begin with decreases and, after 
passing through a minimum, rapidly increases 
corresponding to a steep descent of the curve at 
the very end of the range. The stopping in the 
first part is due almost entirely to electronic 
encounters, whereas the stopping in the later 
part is mainly caused by nuclear collisions, the 
effect of which becomes preponderant when the 
total charge of the fragment is nearly neutralized 
by electron capture. While the electronic stopping 
effect in different gases is proportional to the 
stopping power for a-rays of the same velocity, 
the nuclear stopping effect increases far more 


1K. J. Brostrgm, a gos Lauritsen, Phys. 
Rev. sd 651 (1940) ; As ild, K. J. Brostrgm, 
and T. Lauritsen, Phys. Rev. 58, 350 1940) ; J. K. Boggild, 
K. J. Brostrgm, and T. Lauritsen, Phys. Rev. 59, 275 
(1941); J. K. Bgggild, K. J. Brostrém, and T. Lauritsen, 
D. Kgl. Danske Vid. Selsk. Math.-fys. Medd. (Math.- 
phys. Comm., Acad. Sci. apenas 18, (1940). 

2 N. Bohr, Phys. Rev. 59, 270 (1941). 


rapidly with the atomic number of the gas. This 
means that the end part of the range should be 
relatively longer in a light gas than in a heavy 
one, and in helium it was estimated to be three 
times as long as in argon. Since in this gas the end 
part is roughly 10 percent of the range, it should 
be about 30 percent in helium, and the whole 
range in helium relative to a-rays ought therefore 
to be just about 20 percent longer than in argon. 

In view of the theoretical interest of a closer 
examination of this question the experiments 
were continued, and a large number of cloud- 
chamber pictures of fission fragment tracks in 
helium have recently been obtained which allow 
a direct investigation of the range-velocity curve 
in this gas. A few typical examples of such 
pictures are given in Fig. 1 where, for comparison, 
one of the older photographs in argon is also 
shown. As seen, the tracks in helium are less 
curved than in argon, as was to be expected from 
theory,’ since the mean deviation over a part of 
the range with a given velocity loss should be 
proportional to the square root of the atomic 
number; furthermore, the number of branches of 
a given length is considerably greater in helium 
due to the larger range of recoil particles of the 
same energy in helium than in argon. 

The previous papers contain a discussion of 
two different ways of velocity determination, viz., 
the usual direct method by means of individual 
close collisions, and a statistical method based on 
an examination of the branch frequency along 
the range. The first method has the advantage of 
giving absolute values but the reliability of the 
velocity determination depends upon the cer- 
tainty of the identification of the particles 


3N. Bohr, Phys. Rev. 58, 654 (1940). 
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Fic. 1. A. Track in argon showing fortuitous bending and several 
branches all more intense than the protons and a-ray tracks appearing 
as background. Fine parallel lines are wires fastened to the underside 
of the glass cover to provide conducting surface for sweep field. 
(A+H:0+C:2H,O, total pressure 10 cm). B. Track in helium showing 
a close collision with a helium nucleus rather early in the range and 
numerous collisions on the last part. C. Track with several long helium 
branches suitable for velocity determination especially the two indi- 
cated by arrows. D. Track in helium with at least 30 collisions giving 
branches of more than 2-mm range. The forked track in the center of 
the picture is not a branch but a helium recoil nucleus produced by 
neutron impact and colliding after a path of few mm in the chamber 
with another helium nucleus. (B —-D : He +H.0, total pressure 34 cm.) 


involved in the collision and of the knowledge of 
the range-velocity relation for these particles. 
The branch statistics gives information about the 
range-velocity relation by simply counting in 
various sections of the range the number of 
branches with lengths within given limits. This 
number is, in fact, inversely proportional to the 
square of the average velocity of the fragment in 
the respective range section. This method has the 
advantage that it is independent of other infor- 
mation, but the velocity determinations are only 
relative. 

In the investigation of the range-velocity 
relation in argon the statistical method was 
mainly used, since the mean error of the direct 
velocity determination even under the most 
favorable conditions was about 20-30 percent. In 
helium, however, the direct velocity determi- 
nation is practicable with greater exactness due 
to the fact that the collisions with helium nuclei 
are rather easy to distinguish from collisions with 
heavier and lighter atoms present in the cloud 
chamber. In the present investigation, both 
methods of velocity determination could there- 
fore be used and gave results supporting each 
other. The experimental arrangement was the 
same as that described in the previous papers 
with the only difference that the water-alcohol 
mixtures used at that time were now replaced by 
pure water in order to reduce the stopping power 
of the vapor in the cloud chamber. The uranium 
was placed as a “thick” layer on a copper strip 
near the wall of the chamber. 
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Fic. 2. Range-velocity curve in helium derived from 
individual close collisions with helium nuclei (full circles) 
and with heavier nuclei (open circles). 
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The determination of the velocity of the 
fragments is given by the formula 


i+k)? 7 


4 cos*6 


where & is the ratio of the masses of the gas 
nucleus and the fission fragment, u and v are the 
velocities of these particles after the collision, 
and @ is the angle between the direction of the 
stem of the track before the collision and the 
direction of the branch. In collisions with helium 
nuclei the velocity u was yielded by means of the 
energy-range relation for a-particles,‘ the range 
of the helium branches being converted into 
normal air values by use of the stopping power 
coefficient of the gas. This stopping power was 
measured directly by means of a beam of a-rays 
with a residual range of 12.5-mm normal air, 
admitted through a window in the side of the 
chamber. Since the stopping power of the gas 
turned out to be somewhat larger than it would 
corresponding to pure helium and water vapor, 
the helium gas was examined spectroscopically 
and was found to contain some nitrogen. There- 
fore, in the collisions with heavier nuclei, the 
mass was estimated to be 15, and the velocity u 
was determined by a curve halfway between the 
range-velocity curves for oxygen and nitrogen of 
Blackett and Lees® with the values of the velocity 
increased by 5 percent.® 
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Fic. 3. Range-velocity relation in helium derived from 
statistic of branch distribution. 


4 Cornell University curves (1938). 

’P. M. S. Blackett and D. S. Lees, Proc. Roy. Soc. 
A134, 660 (1931). 

® M. S. Livingston and H. A. Bethe, Rev. Mod. Phys. 9, 
§95 (1937). 
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Fic. 4. Range-velocity relation in argon (curve I) and in 
helium (curve IT). 


The present investigation includes 34 fragment 
tracks at about 25-cm gas pressure and 83 
fragment tracks at about 34-cm pressure. Among 
these, 49 collisions with helium nuclei and 15 
with heavier nuclei are suitable for velocity 
determination; the results of these measurements 
are reproduced in Fig. 2, where the velocity is 
given in relation to the residual fragment range 
converted into normal air. The small full circles 
refer to the collisions with helium nuclei and the 
open circles refer to the collisions with heavier 
nuclei, and it will be seen that the two kinds of 
collisions within the experimental errors give 
consistent velocity determinations. The curve 
drawn on Fig. 2 to fit all the points as closely as 
possible must be taken as an average curve for 
the range-velocity relations for the lighter and 
heavier groups of fission fragments, and the 
relatively small spread of the points shows that 
the velocity curves of the two fragments differ 
only little within the last 10 mm of the range. 

The experimental material was also subjected 
to a statistical examination on the basis of 83 
tracks at 34-cm gas pressure, including in all 
about 1100 branches of length between 2.3 and 
10 mm, i.e., 0.26-1.14 mm normal air. The 
relative velocities thus calculated are reproduced 
in Fig. 3 in such a scale that for the larger 
velocities agreement is obtained with the range- 
velocity curve derived from the direct velocity 
determinations which is here reproduced from 
Fig. 2. The circles with the vertical lines repre- 
sent mean values and probable errors of the 
relative velocities. Although the points fall 
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rather close to the curve, the deviations for the 
smaller velocities are larger than the estimated 
experimental errors indicating some minor de- 
fects in the one or the other of the two methods 
of velocity determination. In the statistical 
method such a defect might perhaps be due to a 
failure to count all branches in the last, more 
branching part of the range; in the direct method 
a defect might be due to a possible inaccuracy in 
the assumed range-velocity relation of the 
a-particles in helium for the smallest velocities. 

A comparison of the range-velocity relation in 
argon and in helium is attempted in Fig. 4. Here, 
curve I is an average range-velocity curve for the 
two groups of fission fragments in argon obtained 
in the previous paper. In curve II, the full- 
drawn part corresponds to the information con- 
tained in Fig. 2, while the dotted part is drawn in 
a way that the total average range is 4-mm 
normal air greater than in argon, in accordance 
with the earlier measurements. Although the 
information on which the upper parts of curves I 
and II are based is rather scanty, it seems that 


there can be no doubt about the general run of 
the curves. In particular, the point where the 
rate of velocity loss is a minimum is obviously 
much farther removed from the end of the range 
in helium than in argon. The shape of the curves 
seems altogether to fit the theoretical expecta- 
tions, and the curves can, at least within the 
experimental errors, be drawn in such a way that 
the rate of velocity loss is alike in the two gases in 
the first part of the range where the stopping is 
mainly due to electron encounters, whereas the 
end part of the range-velocity curve in helium is 
considerably longer than in argon, corresponding 
to the effect of the smaller energy loss by nuclear 
collisions in helium than in argon. 

The author wishes to express his heartiest 
thanks to Professor Niels Bohr for his constant 
interest in this work and for many helpful dis- 
cussions. I am also much obliged to Dr. E. 
Rasmussen for the spectroscopical analysis of the 
helium and to Mr. K. Lindberg-Nielsen and Mr. 
H. Arrge whose assistance in the experiments has 
been of great value. 
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On the Pseudoscalar Mesotron Theory of 3-Decay 


E. C. NELSON 
Department of Physics, University of California, Berkeley, California 
(Received September 30, 1941) 


The calculations of Sakata which indicate that the pseudoscalar mesotron theory can suc- 
cessfully account for both mesotron decay and $-decay are shown to be in error. While the 
pseudoscalar mesotron theory gives the correct spin and parity selection rules for 8-decay, it also 
gives B-decay lifetimes too long by a factor of 10° and having an inverse seventh power depend- 
ence on the upper limit of the 8-spectrum. As these same conclusions can be made when the 
mesotron field is strongly coupled to nuclear particles, it would seem that they are sufficient to 
show that pseudoscalar mesotrons are not responsible for B-decay. 


UKAWA’S! suggestion that 6-decay is to be 
interpreted as the emission of a mesotron by 

a neutron (or proton) and its subsequent decay 
into an electron (or positron) and a neutrino has 
received support from the observed instability? 
and apparent §-radioactivity of cosmic-ray meso- 
1H. Yukawa, Proc. Phys. Math. Soc. Japan 17, 48 (1935). 
2H. V. Neher and H. G. Stever, Phys. Rev. 58, 766 
(1940) ; B. Rossi and D. B. Hall, 7 Rev. 59, 223 (1941); 
W. M. Nielsen, C. M. Ryerson, L. W. Nordheim, and K. Z. 
Morgan, Phys. Rev. 59, 547 (1941); F. Rasetti, Phys. Rev. 


60, 198 (1941); E. J. Williams and G. E. Robert, Nature 
145, 102 (1940). 


trons. However, before it can be regarded as a 
correct description of nuclear 8-processes, it 
must give rise to a theory of 8-decay possessing 
the following features: (1) Fermi energy distri- 
bution* with its consequent inverse fifth power 
dependence of the lifetime on the upper limit of 
the spectrum,‘ (2) Gamow-Teller selection rules,* § 


3 J. L. Lawson, Phys. Rev. 57, 982 (1940). 

4M. G. White, E. C. Creutz, L. A. Delsasso, and R. R. 
Wilson, Phys. Rev. 59, 63 (1941). 

5 T. Bjerge and K. J. Brostrém, Kgl. Danske Vid. Sels. 
Math.-fys. Medd. 16, No. 8 (1938). 
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(3) calculated lifetimes of 8-radioactive nuclei, 
with the mesotron-electron-neutrino coupling 
constant evaluated from the disintegration of the 
free mesotron and the mesotron-nuclear particle 
coupling constant evaluated from nuclear forces, 
in agreement with experiment. The scalar meso- 
tron theory meets only the first point; the vector 
theory but the first two. Both these theories give 
8-decay that is too slow.** Its characteristic 
time, that factor in the lifetime independent of 
the energy and matrix element of the transition, 
is of the order of a hundred times greater than 
that evaluated from experiments on light nuclei. 

Christy and Kusaka’s® analysis of cosmic-ray 
bursts together with the spin dependence of 
nuclear forces indicates that mesotrons are best 
described by a pseudoscalar field.'°" Recently, 
the assertion has been made that pseudoscalar 
mesotrons provide a theory of 8-decay in agree- 
ment with experiment. Calculations by Sakata® 
yielded the necessary ‘‘allowed’’ 8-decay and 
slow disintegration of the free mesotron. How- 
ever, his formula for B-decay turns out to be 
incorrect; proper treatment of the problem 
shows that the lifetime discrepancy remains and 
is, in fact, increased from a factor of 10° to a 
factor of 10°. 

Yukawa treats the connection between nuclear 
forces and $-decay with an interaction of the 
mesotron field with nuclear particles and with 
electrons and neutrinos symmetrical in all re- 
spects except in strength. In the pseudoscalar 
theory this is expressed in the following Lorentz- 


®H. Yukawa, S. Sakata, M. Kobayasi, and M. Taketani, 
Proc. Phys. Math. Soc. Japan 20, 720 (1938). 

a id 0) A. Bethe and L. W. Nordheim, Phys. Rev. 57, 998 

940). 

( s fy Sakata, Proc. Phys. Math. Soc. Japan 23, 291 
1941). 

® 5 F. Christy and S. Kusaka, Phys. Rev. 59, 405, 414 
(1941). 

10 J. R. Oppenheimer, Phys. Rev. 59, 462 (1941). 

1 J. R. Oppenheimer and J. S. Schwinger, Phys. Rev. 60, 
150 (1941). 

12 Sakata has apparently neglected the singularity arising 
from the derivatives of the Yukawa potential in his 
Eq. (79d). As has been pointed out by Serber, reference 7, 
(exp which also occurs 
in the Bethe-Nordheim, reference 7, treatment of the vector 
theory, contains Upon integration 
over the electron-neutrino coordinates in the matrix 
element for 8-decay, this delta function term cancels in the 
usual approximation for ‘‘allowed’’ transitions, wave 
functions of electron and neutrino treated as constant over 
nuclear dimensions, the — (4%/3)@1-@2 from the other terms. 
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invariant terms in the Lagrangian: 


2M stir, 


2M 
(t2+1Ty) 
+g 


(t2—tTy) 


(1) 


U is the pseudoscalar field and U* its complex 
conjugate. y and ¢ are the neutron-proton and 
the electron-neutrino wave fields, y'=iy*8 and 
¢gt=i¢g*8 their respective adjoints. y' and 6 are 
Dirac matrices, rx and ry are isotopic spin 
matrices, /=proton mass, and n= the mesotron 
mass. The mesotron-nuclear particle coupling 
constant g is normalized so that the potential in 
the singlet state of the deuteron is g(e-*’/r), 
c=pc/h. 

It might be expected that interaction terms 
containing the first derivatives of U as well as 
ones containing U should occur in the Lagrangian. 
However, such terms do not really represent a 
different type of interaction, for in virtue of the 
equations of motion of y, they may be derived 
from the terms in (1) by partial integration. 


ha, 


ty 
2 2Mc 2 


The neglected terms are of the order g’, are 
bilinear in the U fields, and give no contribution 
to B-decay. They also give no contribution to 
nuclear forces in the usual approximations. Thus 
there is no need to introduce two independent 
constants, f and g, in the pseudoscalar theory as 
there is in the vector theory. 

Since (1) does not contain the time derivative 
of U, i.e., the canonically conjugate field mo- 
mentum, the interaction terms in the Hamiltonian 
are just the negative of the ones in L. 


int (2) 


In this formulation there are no terms directly 
coupling neutrons and protons to electrons and 
neutrinos in either the Hamiltonian or the 
Lagrangian. Such terms do appear in the 
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Hamiltonian in the equivalent formulation in 
which derivatives of U occur in the interaction 
Lagrangian. There the inclusion of such terms in 
the Hamiltonian is essential to preserve the 
covariance of the theory and to give the same 
matrix elements as (2). Any attempt to eliminate 
them by adding invariant direct interaction 
terms to the Lagrangian succeeds only in re- 
placing them by larger ones. 

The matrix element, /7;;, for B-decay follows 
from (2) by a straightforward calculation using 
second-order perturbation theory, and plane 
Dirac waves for the electron and neutrino. 


= (3) 


2Mi 


fa 


XW, exp[—ik-r.],  (3’) 
W; and W; are the initial and final wave functions 
of the nucleus. The summation s is carried out 
over all the particles in the nucleus. “, and u, are 
the spin functions of the electron and the 
neutrino. 


W = upper limit of the 8-spectrum, 
E,,=energy of the mesotron, 
hk =its momentum. 


Since by the equations of motion of y in non- 
relativistic approximation 


(tr, +17,) h 
(4) 
2 21Mc 2 


hk 


Cues 


hk 
XV; exp{ —ik-r, ]=—-B, (4’) 
Ch 


@ represents the Pauli spin matrices. The transi- 
tion probability for B-decay is: 


1 27 
me 


P.2dP, P, dP, 
dW 


E.=energy of the electron, 
P,=momentum of the electron, 
E,=energy of the neutrino, 
P,=momentum of the neutrino, 


P.+P,=hk, 
(5) 
E.E, 
1 11 
-=— -|B,|*J (wo), (6) 
T 183 


32 we? 
ta 


m=mass of the electron, 


I(wo) = fe = 


Xw(w?— 1)*(wo—w)? 


13wo'* 2 


= 1) 
70 140 105 


In[ wot (8) 


w=Ee/mc*, wo=W/mce*. 


It is instructive to compare B, the nuclear 
matrix element, and B’, the electron-neutrino 
matrix element, (3) in the special case of the 
decay of the neutron. Here all the particles may 
be appropriately described by using plane waves. 
The contribution of B’ is expressed in (5). Since 
both the neutron and proton have non-vanishing 
rest mass, the contribution of B differs from that 
of B’ by mass terms. 


2M\? Py -Ppc?+ My M pct 


spins of P and N Ex Ep 


(Pp)*c? 


while (5) is of the order unity. 
The presence of @ in B (4’) ensures Gamow- 
Teller selection rules. Although By* is a pseudo- 
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scalar, there is no parity change‘in the nuclear 
state in the most “allowed” transition. This 
is so because B vanishes in view of (4) when 
exp [—ik-r] is replaced by the first term, 1, in 
its expansion and the most “allowed” transition 
comes from the second term in the expansion, 
—ik-r, which when multiplied by By® does not 
change sign under reflection of the space coordi- 
nates. Since B is proportional to the momentum 
of the mesotron, the lifetime 7 has an inverse 
seventh power dependence on the upper limit of 
the 6-spectrum in disagreement with the inverse 
fifth power dependence shown by the series of 
light positron emitters.‘ 

Only slow mesotrons are involved in B-decay, 
hk<yc. Thus the Ak/yc in B contributes, upon 
averaging over the angles of the electron and 
neutrino, effectively (m/u)* to the transition 
probability, slowing up §-decay in that pro- 
portion. 

The probability per unit time that a free 
pseudoscalar mesotron decay, obtained after a 
similar but simpler calculation, is 


1/7, =(uc*/h)(g’*/he) 


from (7) and (9) 
g? “) 


with 7,=2X10-* sec., g?/hc=0.1, »=200 m, 

= 2.4(10)'° sec. while from (6) and experi- 
ments on z°*P=6.7(10)* sec. 
~10°, a discrepancy greater by (u/m)? than that 
in the scalar or in the vector theory. 

In deriving these results perturbation-theoretic 
methods have been used. The large mesotron- 
nuclear particle coupling renders these methods 
unsuited to many problems. However, recent 
investigation of this interaction in the limit of 
strong coupling, with extended sources used to 
ensure convergence, indicates that the normal 
state of the field of pseudoscalar mesotrons 
“bound” to nuclear particles differs from the 
classical static field only by a factor of 1/./6." 
Since the classical field gives for 8B-decay the 
same results as perturbation theory, no new 
feature, except the § which only increases the 
lifetime discrepancy, is introduced into this 


(9) 


3 The author is indebted to Professor J. R. Oppenheimer 
for information on this point. 
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problem by strong coupling. For allowed transi- 
tions B-decay does not depend on the details of 
the source. As preliminary results of strong 
coupling calculations indicate that g*/hc is not a 
large number," it is unlikely that a detailed 
treatment of the problem would produce results 
altered enough to give a correct description of 
B-decay. Therefore, it would seem that the 
magnitude of the lifetime discrepancy and the 
seventh power energy dependence are sufficient 
to show that pseudoscalar mesotrons are not 
responsible for 8-decay. 

Two alternative interpretations of 8-decay 
have been suggested: (1) There are two kinds of 
mesotrons—long-lived pseudoscalar mesotrons 
observed near sea level in cosmic rays and fast 
decaying mesotrons of another type to give 
B-decay.'* Mgller and Rosenfeld'® suggested, for 
other reasons, a theory giving a combination of 
pseudoscalar and vector mesotron fields in nuclei. 
Rozental has shown that the constants in this 
theory can be chosen to give correctly the 
observed mesotron decay and #-decay.'® (2) 
B-decay occurs as in the original Fermi theory 
by means of a direct interaction between the 
nuclear particles and electrons and neutrinos. 
Since the constant characterizing this interaction 
is dimensionally different from that which couples 
mesotrons to electrons and neutrinos, it does not 
seem possible to relate them in any unique way 
and establish thereby a simple connection be- 
tween 8-decay and the decay of a free mesotron. 
The suggestion’? that mesotron decay occurs 
through the virtual emission of a neutron-proton 
pair is dependent on a highly speculative appli- 
cation of a not yet substanstiated relativistic 
theory of neutrons and protons. 

The author wishes to thank Professor J. R. 
Oppenheimer for suggesting this problem and for 
his continued interest in this work and Dr. J. S. 
Schwinger and Dr. E. P. Cooper for their helpful 
discussions. 


“4 R. F. Christy and J. R. Oppenheimer, Phys. Rev. 60, 
566 (1941). 

4% C, Mgller and L. Rosenfeld, Kgl. Danske Vid. Sels. 
Math.-fys. Medd. 17, No. 8 (1940); C. M@ller, Kgl. Danske 
Vid. Sels. Math.-fys. Medd. 18, No. 6 (1941). 

16S. Rozental, Kgl. Danske Vid. Sels. Math. -fys. Medd. 
18, No. 7 (1941). S. Rozental, Phys. Rev. 60, 612 (1941). 
(194t} Sakata, Proc. Phys. Math. Soc. Japan 23, 2 
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Letters to the Editor 


ROMPT publication of brief reports of important dis- 
coveries in physics may be secured by addressing them 
to this department. Closing dates for this department are, for 
the first issue of the month, the eighteenth of the preceding 
month, for the second issue, the third of the month. Because 
of the late closing dates for the section no proof can be shown 
to authors. The Board of Editors does not hold itself responsible 
for the opinions expressed by the correspondents. 
Communications should not in general exceed 600 words 
in length. 


Slowing Down of Low Energy Neutrons 
in Water 


M. GOLDHABER AND R. D. O'NEAL 
Department of Physics, University of Illinois, Urbana, Illinois 
October 15, 1941 


HE slowing down in water of the fast neutrons 

emitted by a Rn-a-Be source has been studied ex- 
tensively by many investigators, particularly by Amaldi 
and Fermi.! Figure 1a, taken from their work, shows how 
the intensity inside water of neutrons of thermal energy 
(C neutrons) and neutrons of an energy of a few ev (D, 
A+B, I, in order of increasing energies) varies with the 
distance from the source. The striking feature, that at 
large distances from the source the curves for C and reso- 
nance neutrons are rather similar, was explained quali- 
tatively as due to the fact that the average mean free path 
for scattering of the fast primary neutrons which have 
energies up to 12 Mev is of the order of 5 cm, whereas that 
of the resonance neutrons is less than 1 cm, and the diffu- 
sion length of the C neutrons is 2.1 cm. At some distance 
from the source, equilibrium conditions are approached: 
The character of the curves is now governed in the main 
by the largest of these three values, i.e., the mean free 
path of the primary neutrons, and hence the intensities of 
resonance and C neutrons decrease in a similar manner. 
This implies that the intensity of the resonance neutrons 
does at no point in the water become negligible compared 
with that of the C neutrons. This important result has 
often been considered as quite generally true, independent 
of the neutron source used. 

However, a different situation might be expected when 
the neutron source emits neutrons of low, and preferably 
homogeneous, energy. For neutrons of a few hundred kev 
energy, the mean free path for scattering in water is about 
1 cm, only a little larger than for resonance neutrons, but 
well below the diffusion length of € neutrons. For such 
primary neutrons, one might therefore expect to find much 
less similarity between the Amaldi-Fermi curves for reso- 
nance and C neutrons. This is borne out by the following 
experiment. 

As a neutron source we used the photo-neutrons pro- 
duced in beryllium by the hard y-rays of radio-yttrium 
(100d). The RdY, of an initial intensity of about 30 
millicuries Ra equivalent, was produced in the 60-inch 
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Berkeley cyclotron by the reaction Sr(d,2n)Y.2. The 
RdY+Be neutrons have an energy of 220+20 kev, as 
determined from the energy of the recoil protons by 
Scharff-Goldhaber.’ The C and D neutrons were measured 
with an indium foil detector (54 min.) and the I neutrons 
with an iodine detector (PbI,). The results are shown in 
Fig. 1b. It is seen that at a large distance from the neutron 
source, the intensity of the resonance neutrons decreases 
much faster than that of the C neutrons. In this way one 
can obtain C neutrons practically free from resonance 
neutrons. 

We have observed similar effects with other photo- 
neutron sources. A detailed report of this work will appear 
in due course. 


I 
A+B 
D 
Cc 
(a) c 
te] 5 10 15 20 cm 


Fic. 1. Ordinates =activation Xr?, Abscissae =r =distance of de- 
tector from center of source. The units of the ordinates are chosen 
so that each curve includes the same area. (a) Taken from Amaldi 
and Fermi, reference 1, p. 918, Fig. 7. Rn-a-Be neutrons; (b) RdY +Be 
neutrons. 


The radio-yttrium source was produced in Berkeley by 
the late Dr. C. Pecher. We wish to thank Professor E. O. 
Lawrence for generously permitting us to use this valuable 
source indefinitely. 

1 E, Amaldi and E. Fermi, Phys. Rev. 50, 899 (1936). 

2C. Pecher, Phys. Rev. 58, 843 (1940). 

3G. Scharff-Goldhaber, Phys. Rev. 59, 937A (1941); see also J. R 


Richardson, Phys. Rev. 60, 188 (1941); J. R. Downing, M. Deutsch, 
and A. Roberts, Phys. Rev. 60, 470 (1941). 
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The Change in the Electric Resistance of a Single 
Crystal of Magnetite by Magnetic Field 
at Low Temperatures 


H. MASUMOTO AND Y, SHIRAKAWA 


Research Institute for Iron, Steel and Other Metals, Téhoku Imperial 
University, Sendai, Japan 


September 5, 1941 


HE change in electric resistance of a single crystal of 
magnetite at low temperatures, 0° (ice water), 
—95° (freezing point of toluol) and —195° (liquid nitro- 
gen) was measured up to 1200 oersteds in longitudinal 
magnetic field, and the intensity of magnetization of the 
specinien was also measured under the same condition. 
The specimen is a rectangular crystal rod, 3.43 mm? in 
sectional area and 6.0 mm in length, which was cut off 
from a block of perfect natural crystal from Kamaishi, 
Japan; the axis of the rod coincides with the direction of a 
principal axis [111]. Each end of the specimen was soldered 
to a copper lead wire by aid of platinum chloride melted 
in lavender oil. The resistance was measured with a poten- 
tiometer. The amount of current through the specimen at 
0°, —95°, and —195° was 0.02, 0.01, and 0.001 ampere, 
respectively. The intensity of magnetization was measured 
by means of the ballistic method; the demagnetizing factor 
used was that given by R. Mann. 

The results of measurement are shown in Fig. 1, in which 
the lower half is the relation between the change of AR/R 
and the longitudinal magnetic field H, the magneto- 
resistance curve, and the upper half is the relation between 
the intensity of magnetization J and H, the magnetization 
curve. From this figure, we see that the change of resistance 
is negative for all field strengths at temperatures below 
zero. This negative magneto-resistance effect has been 
observed only in Fe-Si alloys,! but not in the alloys of 
Ni-Cu,? Ni-Co, Fe-Ni and Fe-Co.* The change of resistance 
at 0° and —95° decreases at first rapidly and then slowly 


Fic. 1. Magneto-resistance and magnetization curves at low tempera- 
tures for magnetite having principal axis [111]. 
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with the increase of the field, while at — 195°, it decreases 
linearly. The different form in these curves is to be attrib- 
uted to the difference of the states before and after the 
allotropic transformation of magnetite in — 158°~— 166°.* 
By the same cause, the magnetization curves at 0° and 
—95° are different in form from that at —195°. In the 
transversal field the change of resistance shows the same 
behavior as shown in the ferromagnetic metals. 

The decrease of resistance in the ferromagnetic com- 
pound is very interesting and will be studied in the other 
principal axes of magnetite in the near future. 

1Y. Shirakawa, Nature 142, 395 (1938); Sci. Rep. Téhoku Imp. 
Univ. 27, 255 (1938). 

2H. Masumoto, Y. Shirakawa, Sci. Rep. Téhoku Imp. Univ. 25, 104 
ay "Shirakawa, Sci. Rep. Téhoku Imp. Univ., Honda Anniv. Vol., 


362 (1936); 27, 485, 532 (1938). 
4T. Okamura, Sci. Rep. Téhoku Imp. Univ. 21, 231 (1932). 


Meson Theory and Nuclear Forces 


HARALD WERGELAND 
Blindern Universitet, Oslo, Norway 
October 9, 1941 

S is well known, the nuclear interaction predicted by 
meson theory contains terms so strongly singular 
(~-—1/r’) that some kind of cutting off procedure must 
be adopted when applying it as potential energy in the 

Schrédinger equation. 

However, there is a second possibility discussed by 
Bethe,! but finally discarded by him, namely, that the 
diverging terms may have a significance analogous to the 
magnetic interaction between electrons. In that case the 
situation would be entirely different, because, as Breit has 
shown,? such terms cannot simply be considered as part 
of an energy operator in configuration space. 

In fact, there are reasons for believing that the meson 
interaction and Breit’s operator have a similar origin. To 
demonstrate this, we modify the derivation of nuclear 
forces from field theory in the following way: 

If the homogeneous equations governing the meson field 
in vacuum are 

one may ask what are the inhomogeneous equations for 
the field generated by the heavy particles. 

Clearly the right-hand side must be a four vector. Now 
the simplest four vector which can be associated with the 
motion of the heavy particles is the Dirac current 


ay; 
According to Gordon,’ ¢ can be resolved into two parts, 
viz., the convection current 


oy oy* 
ve 


(1,2 Dirac matrices, a4= 3). 


and the polarization current 


each of them separately satisfying the continuity equation 
0j,/dx, 
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A general wave equation would therefore be, 


Oxyox" 


with +; and +2 as independent constants. 

This together with the wave equation for the free heavy 
particles defines practically uniquely a Lagrangian and 
Hamiltonian in which the interaction heavy particle-meson 
field is represented by: 


+ 


Proceeding to the non-relativistic limit for the heavy par- 
ticles, neglecting terms ~v/c and leaving out 6-functions, 
one gets in the usual way the interaction potential: 


With yi=g1, v2=(2M/u)ge, (u=meson mass) this is the 
form given by Youkawa et al. 

If we take y: = y2=e and let ,—0, this expression changes 
continually into its electrostatic counterpart 


r? 


that is, the Coulomb interaction+the spin-spin inter- 
action which follows from Breit'’s operator in non-rela- 
tivistic approximation. 

From the point of view proposed above one would 
expect 7:= 2 also in meson theory, thus leaving only one 
constant free for each sort of field. Only tnen will the 
difference between electrodynamics and meson theory 
merely come from the finite rest mass of the mesons. On 
account of the conservation of the currents it is of course 
formally possible to choose y:#y2, but this is a separate 
assumption which cannot be inferred from the analogy 
with electrodynamics. 

A. Bethe, Phys. Rev. 57, 260 (1940). 


2G. Breit, Phys. Rev. 39, 616 (1932). 
3 W. Gordon, Zeits. f. Physik 50, 630 (1928). 


On the Gyromagnetic Effects in 
Ferromagnetic Substances 


C. J. GoRTER 
Zeeman-Laboratorium der Universiteit van Amsterdam, Holland 
September 16, 1941 


HE refined researches of Barnett and others! on the 
gyromagnetic effects in ferromagnetic substances 
have revealed positive deviations of the ratio p of mechani- 
cal to magnetic moment from the spin-only value mc/e. 
This obviously means that, besides the spins, the orbital 
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motion of electrons also contributes to the mechanical 
and magnetic moments. But further understanding was 
blocked by the apparent influence of minor structural 
differences in the samples upon the magnitude of the 
deviation. The most striking example was that the devia- 
tion mentioned for Yensen electrolytic iron, as observed 
in the Barnett effect, was three times larger than that for 
soft Armco or Norway iron. Even more incomprehensible 
was the difference between the gyromagnetic ratio p for 
Yensen iron and also for an iron-cobalt alloy, observed in 
the Barnett effect, and the lower gyromagnetic ratio ob- 
served in the Einstein-de Haas effect for the same 
substances. 

Now Barnett recently in a footnote? sheds doubt upon 
his own former results on the Barnett effect of Yensen 
iron and the iron-cobalt alloy, mentioning that a redetermi- 
nation appears to give much lower values of p, in agreement 
with those obtained in the Einstein - de Haas effect. It is 
the purpose of this note to point out that this considerably 
simplifies the situation. Considering Barnett’s results and 
omitting the doubtful data, one sees that the p-values 
seem to split into four different groups, namely, nickel and 
the alloys containing mainly nickel, cobalt and the alloys 
consisting primarily of cobalt, iron and the alloys con- 
taining mainly iron, and finally Heusler’s alloy. The 
approximate values of pe/mc for each group will be found 
in Table I. 


TABLE I 
pe/me pe ‘mc 
SUBSTANCES EXP CALC. 10N 
Heusler's alloy 1.01 1.00 Mn** 
Iron and alloys 1.03 1.07-1.11 Fe** 
Cobalt and alloys 1.08 1.12-1.27 Cot* 
Nickel and alloys 1.05 1.13-1.18 Nitt 


On the other hand, recently the gyromagnetic ratios of 
bivalent paramagnetic salts have been calculated* from the 
deviations of their susceptibilities from the spin-only 
values. In the right-hand side of the table these calculated 
values of pe/me are given. 

Apparently there is a parallelism between the two series 
of numbers, though the theoretical paramagnetic deviations 
from unity are about three times larger than the experi- 
mental ferromagnetic deviations. Does this mean that the 
ions responsible for ferromagnetism are doubly ionized and 
subject to an interaction with their surroundings which is 
similar to (though a few times larger than) that charac- 
teristic of the bivalent ions in paramagnetic salts? 


1S. J. Barnett, Rev. Mod. Phys. 7, 129 (1935); D. P. Ray-Chaudhury, 
Ind. J. Phys. 9, 383 (1935); I. K. Kokoin and S. V. Goobar, C. R. 
Acad. Sci., Moscow 19, 249 (1938). 

2S. J. Barnett, Proc. Am. Acad. 73, 401 (1940). 

3C. J. Gorter and B. Kahn, Physica 7, 753 (1940). 
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